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FOREWORD 


TEris  report  was  prepared  by  the  Douglas  Aircraft  Company,  Aircraft 
Division,  Long  Beach,  California,  for  the  Aerospace  Dynamics  Branch,  Vehicle 
Dynamics  Division,  Air  Force  Flight  Dynamics  Laboratory,  Wrigr.t-Patterson 
Air  Force  Base,  Ohio  under  Contract  F33615-70-C-1167.  This  research  was 
conducted  under  Project  137°,  "Dynamic  Problems  in  Military  Flight  Vehicles", 
and  Task  137003,  "Prevention  of  Dynamic  Aeroelastic  Instabilities  in  Advanced 
Military  Aircraft,"  Mr.  S.  J.  Pol’ock  of  the  Aerospace  Dynamics  Branch  was 
Task  Engineer. 

This  report  consists  of  two  parts  with  two  volumes  for  ev'.ch  part.  This 
volume.  Volume  II  of  Part  II  is  the  Computer  Program  N5KA.  Volume  I  of 
Part  II  contains  a  method  which  uses  an  image  system  and  an  axial  singularity 
system  to  account  for  the  effects  of  the  bodies.  Volume  I  of  Part  I  contains 
the  method  of  direct  application  of  nonplanar  lifting  surface  elements,  and 
Volume  II  of  Part  I  is  the  Computer  Program  H7WC. 

The  work  reported  herein  was  conducted  during  the  period  of  December 
1969  to  August  i971. 

The  Principal  Investigator  was  Joseph  P.  Giesing.  Mrs.  T.  P.  Kalman  was 
responsible  for  the  computer  programming  and  Dr.  W.  P.  Rodden  was  a  McDonnell 
Douglas  Company  Consultant.  Others  have  made  significant  contributions  to  this 
project  including  Messrs,  D.  K,  Larson,  D,  S.  Warren,  and  W,  E.  Henry. 

The  contractor's  designation  of  this  report  is  MDC-J0944.  The  report 
was  released  by  the  authors  in  August  1971  for  publication  as  an  AFFDL 
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This  technical  report  has  been  reviewed  and  is  approved, 

WALTER/ J 

Asst,  for  Research  &  Technology 

Vehicle  Dynamics  Division 


ABSTRACT 


A  technique  for  predicting  steady  and  oscillatory  aerodynamic  loads  on 
general  configurations  has  been  developed  which  is  based  on  the  Doublet-Lattice 
Method  and  the  method  of  images.  Chord-  and  spanwise  loading  on  lifting 
surfaces  and  longitudinal  body  load  distributions  are  determined.  Configura¬ 
tions  may  be  composed  of  an  assemblage  of  bodies  (elliptic  cross  sections  and 
a  distribution  of  width  or  radius)  and  lifting  surfaces  (arbitrary  planform 
and  dihedral,  with  or  without  control  surfaces).  Loadings  predicted  by  this 
method  are  required  for  flutter,  gust,  frequency  response  and  static  aero- 
elastic  analyses  and  may  be  used  to  determine  static  and  dynamic  stability 
derivatives.  Volume  I  presents  the  theory  and  calculated  results  while 
Volume  II  presents  the  details  of  the  computer  program  used  to  implement 
the  theory. 
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NOMENCLATURE 


Reference  total  area 
Average  body  width 
Local  body  width 
Radius  of  curvature 
Average  body  height 
Local  body  height 

Position  vector  to  center  of  c  ir'ature 

Rolling  moment  coefficient  (moiiprst  qA  )  (■•  right  wing  down) 

Pitching  moment  coefficient  (momont/qA-)  . f  nose  up) 

Yawing  moment  coefficient  (moment/qr.)  (+  nose  right) 

Pressure  coefficient 

Side  force  coefficient  (Force/qA)  (+  cut  light  wing) 

Vertical  force  coefficient  (Foroe/oA)  (+  ip) 

Local  chord  length 
Reference  chord  length 
Local  pitching  moment  coefficient 
Loval  normal  force  coefficient 
Center  of  pressure 

Matrix  relating  normal wash  to  lifting  pressures  for  lifting 
surface  elements 

Matrix  relating  normalwash  to  lifting  pressures  for  image 
elements 

Matrix  relating  normalwash  to  lifting  pressures  for  elements 
and  all  their  images 

Matrix  relating  normalwash  to  lifting  pressures  for  elements 
plus  their  images  plus  the  contributions  due  to  symmetry  and 
ground  effect 


Dj  Partitioned  matrix  [D  |  E],  relating  normalwash  to  lifting 

pressures  and  doublet  strengths 

D  Matrix  relating  the  flow  normal  to  a  body  surface  (at  the 

y  meridian  angle  e)  to  the  lifting  pressure  for  elements  and 

their  images 

D(y),D(z)  Matrix  relating  the  average  side-  or  upwash  at  a  body  due  to 
lifting  surface  elements 

D2D  Matrix  relating  the  doublet  strength  to  the  local  up-  or  side- 

wash  using  quasi -steady,  two-dimensional  slender-body  theory 

d  Spacing  of  doublets  or  vortices  within  slender  bodies 

(simulation  of  body  aspect  ration  (b/a)) 

E  Matrix  relating  normalwash  to  axial  doublet  strengths 

E  Matrix  relating  normalwash  to  axial  doublet  strengths  with  the 

effects  of  symmetry  and  ground  effect  included 

E^y),E(z)  Matrix  relating  the  normalwash  to  y-  or  z-oriented  axial  doublets 

e  Lifting  surface  element  semi-width;  also  cross-sectional 

element  semi -width 

F<z>,F(y)  Total  force  on  a  body  due  to  a  point  pressure  doublet.  Subscript 
p(z)  F(y)  indicates  direction  of  force;  superscript  indicates  direction  of 
y  ,hy  pressure  doublet 

f  Nondimensional  deflection.  Also  function  involving  Hankel  functions 

(2) 

'  Hankel  function  of  the  second  kind  of  order  v 

V 

h  Deflections  normal  to  a  lifting  surface 

hy,h2  Deflections  of  a  body  in  y-  and  z-directions,  respectively 

i,T,k  Unit  vectors  in  x-,  y-  and  z-directions,  respectively 

ip  Unit  vector  in  the  direction  of  the  body  force 

K  Velocity  kernel  function;  the  normalwash  due  to  a  point  pressure 

doublet;  also  (a£  -  b£)/4 

K,  Potential  kernel  function;  the  potential  due  to  a  point  pressure 

♦  doublet 


kr  Reduced  frequency  (wc/EU^) 

k  urM/U^ 
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The  normal wash  due  to  a  potential  doublet 
The  potential  due  to  a  potential  doublet 

Mach  number;  also  normalwash  due  to  a  point  source;  also  moment 
Orientation  of  pressure  doublet 
Outward  normal  and  tangent  vectors 
Function  involving  Hankel  functions 
Generalized  force;  also  modified  acceleration  potential 
Dynamic  pressure 
Generalized  modal  coordinate 
fix  -  c)2  +  e2r2 
[( y  -  n)2  +  (z  -  c)2 
(a  +  b)/2 
semi -span 

Freestream  velocity 

Normalwash  boundary  values 

Normalwash  due  to  image  lifting  surface  elements 

Normalwash  due  to  body  interference  doublet  distribution 

Ws  +  W, 

Normalwash  due  to  lifting  surface  elements 
w  -  AW 

Normalwash  in  the  circle  plane 
Coordinates  of  a  receiving  point 
Coordinate  about  which  moments  are  taken 
angle  of  attack 
/(  -  M2 

Dihedral  angle:  yr,  receiving  point,  y$,  sending  point 
Vortex  strength 
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Lifting  pressure 

Modified  acceleration  potential  jump 
Normalwash  due  to  slender  body  elements 
Longitudinal  length  of  lifting  surface  box 
Longitudinal  length  of  axial  element 
Potential  jump 

Symmetry  plane  indication  (1  symmetry,  0  no  symmetry,  -1  anti¬ 
symmetry);  also  a  delta  function;  also  a  virtual  displacement 

Elemental  area 

Ground-effect  indication  (-1  ground  effect,  0  no  ground  effect, 

1  anti ground  effect) 

z-coordinate  of  sending  point 

y-coordinate  of  sending  point 

Lateral  coordinates  in  the  plane  of  the  lifting  surface 

Meridian  angle  for  a  body  of  circular  cross  section 

Sweep  of  1/4-chord  of  lifting  surface  element;  also  inclination 
angle  in  z-y-plane  of  a  cross-sectional  surface  element 

Quadrupole  strength 

Doublet  strength  of  interference-body  elements 

Doublet  strength  of  slender-body  elements 

Multipole  strength  in  circle  plane;  v  gives  order  of  pole 

Doublet  strength  of  modified  acceleration  potential  distribution 
in  y-  and  z-di recti ons;  also  reduction  factors  for  image  doublets 

x-coordinate  of  sending  point 

Distance  from  center  curvature  to  external  singularity 

Source  strength 

Velocity  potential 

Acceleration  potential 

Frequency 


x 


5  Center  of  axial -body  element 

v» 

?1  Leading  edge  of  body  element 

^2  Trailing  edge  of  body  element 

Subscripts  and  Superscripts 
a  Body  axis 

b  Body 

I  Image 

LL  Lower  left-hand  quadrant 

LR  Lower  right-hand  quadrant 

n  Residual  or  interference  flow 

r,s  Receiving  and  sending  points,  respectively 

UL  Upper  left-hand  quadrant 

UR  Upper  right-hand  quadrant 

s  Steady 

y,q  y-  and  z-di recti ons 

e  On  the  body  surface 

1,2  Planar  and  nonplanar  parts,  respectively 

1/4  Quarter  chord  of  element 
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1 .0  INTRODUCTION 


Program  N5KA  is  the  result  of  implementing  the  equations  of  Section  2.0, 
Part  II,  Vol .  I  (Reference  1)  for  the  computer.  The  organization  of  these 
computations  is  outlined  as  follows: 

1)  All  of  the  data  required  (except  modal  data)  are  generated  from 
the  input  data  in  program  Segment  No.  2. 

2)  The  influence  coefficient  matrix  [Dy]  is  generated  in  program 
Segment  No.  3.  This  matrix  relates  the  normalwash,  upwash  and 
sidewash  to  the  lifting  surface  pressures  and  body  axial  doublet 
strengths  in  the  z-  and  y-directions. 

3)  The  normalwash,  sidewash  and  upwash  flow  fields,  aw,  caused  by  the 
slender  body  elements  is  generated  in  program  Segment  No.  5  for 
all  modes.  Modal  data  is  read  in  and  organized  in  program  Segment 
No.  4.  Currently  the  source  distributions  for  steady  flow  are 
not  included. 

4)  The  final  normalwash  (also  side  and  upwash)  boundary  condition,  Wy, 
caused  by  the  motions  of  the  lifting  surfaces  and  the  slender  body 
flow  field,  aw,  is  generated  in  program  Segment  No.  6. 

5)  The  augmented  matrix  [Dy!  Wy]  is  formed  in  program  Segment  1  and 
solved  in  program  Segment  7  for  the  lifting  surface  pressures  and 
body  doublet  strengths  (both  in  z-  and  y-directions). 

6)  The  solution  obtained  in  program  Segment  7  is  used  to  calculate  the 
body  forces  and  moments  in  either  program  Segment  8  or  program 
Segment  9  depending  on  the  method  of  calculations  desired. 

7)  The  lifting  surface  pressures  and  body  axial  force  and  moment 
distributions  are  integrated  to  form  aerodynamic  coefficients  in 
program  Segment  No.  10.  Also  in  this  segment  the  body  forces  are 
redistributed. 

8)  The  lifting  surface  pressures  and  body  axial  force  and  moment 
distributions  are  integrated  to  form  the  generalized  forces  in 
program  Segment  11. 

Configurations  are  built  up  of  lifting  surface  panels  and  bodies  to  any 
degree  of  complexity  desired  (see  Sketch  1.0-1).  A  lifting  surface  panel  is  a 
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SKETCH  1.0-1 
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trapezoid  (two  edges  parallel  to  the  x-axis)  defined  by  the  coordinates  of 
its  corner  points.  The  panel  is  divided  arbitrarily  chordwise  and  spanwise 
to  produce  a  surface  of  elements  or  boxes.  Hinge  lines,  fold  lines,  and  lines 
of  intersection  (of  two  or  more  panels)  must  lie  on  box  boundaries.  Lifting 
surface  panels  that  lie  one  in  back  of  another  must  be  aligned  so  that  those 

2 

box  edges  lying  in  the  same  plane  coincide.  Part  I,  Volume  I  of  this  report 
provides  additional  information  on  the  distribution  of  lifting  surface  boxes  ( 
pages  49-50).  Bodies  are  defined  in  two  ways:  1)  a  tube  of  constant  cross 
section,  of  aspect  ratio  b/a,  divided  into  interference  elements  and  2)  an 
equivalent  body  of  elliptic  cross  section  of  varying  radius  divided  up  into 
slender  body  elements.  The  interference  elements  are  provided  to  complete  the 
body  interference  and  are  concentrated  in  regions  of  maximum  interference  while 
the  slender  body  elements  and  radius  distribution  are  used  to  obtain  the 
slender  body  flow  field  of  the  body.  It  is  important  that  the  radius  distri¬ 
bution  be  accurate  because  a  numerical  derivative  of  the  radius  must  be  taken 
with  respect  to  x.  Also  provided,  for  bodies,  are  one  or  two  distributions 
of  pickup  points.  Pickup  points  are  the  points  on  the  body  surface  where  the 
normalwash  is  to  be  determined  in  order  to  calculate  an  average  up-  or  side 

wash.  The  distribution  of  pickup  points  is  given  in  terms  of  e,  which  for  a 

•1  az 

circle  is  the  meridian  angle.  For  an  ellipse  e  is  defined  as  tan  (^-). 

There  are  two  possible  e  distributions;  one  for  regions  where  many  (e.g.  8) 
points  are  required  and  regions  where  few  (e.g.,  4)  are  required. 

When  the  actual  lifting-surface/body  geometry  is  built  up  the  constant 
section  tube  is  used  to  represent  the  body.  It  is  important  that  the  lifting 
surfaces  adjacent  to  bodies  be  attached  to  the  body  surfaces  without  any  gaps. 
Even  small  gaps  cause  a  reduction  in  load  near  the  gap.  The  relationship 
between  lifting  surface  panels  and  bodies  may  be  incorrect  because  the 
average  tube  is  used  instead  of  the  actual  body  shape.  In  order  to  correct 
this  difficulty  and  obtain  the  proper  flow  field  at  the  lifting  surface  panel 
due  to  the  slender  body  elements,  a  shifting  of  the  panel  is  provided  for. 

The  shift  of  a  panel  may  be  different  for  different  bodies.  Thus  the  quantities 
An^p^,  a?^  are  input  into  the  program,  where  £  indicates  panel  and  b_  indicates 
body. 

To  increase  the  efficiency  of  the  method  a  further  input  is  required 
for  each  panel.  This  input  identifies  the  bodies  "associated"  with  that  panel. 
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Associated  bodies  are  those  that  are  required  to  possess  an  image  of  the  panel 
in  question.  It  is  important  to  have  only  as  many  images  as  is  necessary  for 
accuracy  since  each  image  doubles  the  number  of  kernel  functions  that  are  to 
be  calculated.  It  may  even  be  desirable  to  break  up  lifting  surfaces  into 
more  panels  than  necessary  so  that  panels,  located  a  considerable  distance 
from  the  body,  may  not  be  required  to  possess  an  image  within  that  body. 

Symmetry  planes  and  ground  effect  are  included  by  the  use  of  the  input 
quantities  6  and  e.  Symmetry,  no  symmetry  and  antisyrrmetry  are  activated  by 
setting  6  equal  to  1,  0  and  -1  respectively.  The  plane  of  symmetry  is  the 
y  =  0  plane.  Ground  effect,  no  ground  effect  and  antiground  effect  are 
activated  by  setting  e  equal  to  -1,  0  and  1  respectively.  The  ground  effect 
plane  is  the  z  =  0  plane.  Thus  if  ground  effect  is  desired  the  configuration 
must  be  placed  a  distance  above  the  z  =  0  plane. 

Polynomial  modes  are  input  for  panels  (motion  normal  to  the  surface)  and 
for  bodies  (motions  in  the  z-  and  y-di recti ons ) .  The  coefficients  of  the 
polynomials  are  the  input  quantities  to  the  program. 

The  following  list  gives  the  program  limits: 

1)  The  maximum  number  of  unknowns,  i.e.  the  total  number  of  all  the 
lifting  surface  elements  plus  the  interference  body  elements,  is  500. 

NWORK 

2)  The  maximum  number  of  modes  is  either  50  or  whichever  is 

smaller,  where  NWORK  =  10000  and  NT0T  is  the  total  number  of  unknowns 
for  the  case. 

3)  The  maximum  number  of  panels  is  99  while  the  maximum  number  of  bodies 
is  10. 

4)  The  maximum  number  of  spanwise  strips  per  panel  is  50  while  the 
maximum  number  of  chordwise  boxes  per  strip  is  50. 

5)  The  maximum  number  of  body  interference  elements  (for  all  bodies) 
is  100  while  the  maximum  number  of  slender  body  elements  is  200. 

The  work  area  dimension,  NWORK,  imposes  the  following  overall  restri¬ 
ction:  4  (NSTRIP  +  NBZ  +  NBY  +  NTZS  +  NTYS)  must  be  <  NWORK,  where 
NSTRIP  is  the  total  number  of  strips,  NBZ  is  the  number  of  z-oriented 
bodies,  NBY  is  the  number  of  y-oriented  oodies,  and  NTZS,  NTYS  are 
the  total  number  of  slender  body  elements  with  z-  and  y-orientations 
respecti vely. 
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6)  The  maximum  number  of  modal  coefficients  for  panels  is  150.  The 
maximum  number  of  modal  coefficients  for  bodies  is  150  for  z-motions 
and  150  for  body  y-motions. 

7)  The  maximum  number  of  reduced  frequencies  is  6. 

It  is  important  that  none  of  the  above  maximum  values  is  exceeded.  Not 
all  maximums  can  be  utilized  at  the  same  time  without  violating  others.  For 
example  if  the  maximum  number  of  spanwise  strips  and  chordwise  boxes  is  input 
then  the  total  number  of  lifting  surface  boxes  will  exceed  the  maximum  number 
of  unknowns.  The  maximums  outlined  above  are  tailored  to  allow  the  computer 
program  to  fit  into  a  core  (360/65)  of  260  k  bytes.  If  more  core  is  available, 
the  user  may  wish  to  increase  the  dimension  (NW0RK)  of  the  work  array  (WORK, 
see  Section  5.5.1)  in  order  to  accommodate  larger  cases  in  the  program. 
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2.0  INPUT  PROCEDURE  AND  EXAMPLE  CASE 


2.1  Input  Sheets 

The  input  sheets  for  program  N5KA  are  shown  on  the  next  two  pages.  The 
first  three  cards  represent  general  data  that  is  input  once  per  case.  The  next 
four*  cards  {sequence  numbers  4,  5,  6  and  7)  represent  panel  data  that  is 
repeated  per  panel . 

If  the  data  includes  bodies  also,  the  next  input  card  (#8)  contains 
general  information  for  the  first  body  of  the  case.  The  subsequent  two**  cards 
(sequence  numbers  9  and  10)  are  interference  body  element  data,  and  the  next 
two**  cards  (#11  and  #12)  represent  slender  body  element  data.  Next,  the 
0i-  and  02-arrays  are  input  (cards  #13  and  14)  which  describe  the  angular 
distribution  of  points  on  the  surface  of  interference  body  elements.  Card 
15  identifies  the  sections  of  the  interference  body  for  which  01-distribution 
is  specified.  The  next  card  (#16)  contains  the  y-  and  z-shift  information 
for  all  panels  if  this  is  also  desired  (see  control  items  in  card  8). 

Cards  8  through  16  represent  all  the  body  data  required  and  are  repeated 
per  body.  Note,  that  bodies  oscillating  in  the  z-direction  are  input  first, 
then  bodies  that  can  oscillate  both  in  the  z-  and  y-directions,  and  finally, 
bodies  that  oscillate  in  the  y-direction  only. 

4» 

The  last  four  cards  represent  the  polynomial  mode  information  for  the 
case.  Card  17  contains  the  '-odal  data  for  z-oriented  bodies,  card  18  contains 
the  same  for  y-oriented  bodio,  and  card  19  represents  modal  data  for  panels. 
Modal  data  may  be  input  in  any  order  with  a  maximum  of  three  sets  of  data  per 
card,  a  total  of  100  sets  of  data  per  case,  and  only  the  nonzero  modal  coeffi¬ 
cients  need  to  be  input.  Card  20  terminates  the  reading  of  modal  data;  always 
input  card  20  as  the  last  card  following  all  modal  information.  If  the  case 
specifies  more  than  one  reduced  frequency,  all  modal  data  cards  are  to  be 
repeated  for  each  reduced  frequency  (see  card  #3). 

A  detailed  description  of  all  data  items  is  given  following  the  input 
sheets  for  program  N5KA. 

*  There  may  be  more  than  four  cards  in  order  to  present  all  0  and  t  values. 

**  There  may  be  more  than  two  cards  in  order  to  present  oil  ?I,  RI  elements 

for  interference  bodies,  and  ?S,  RS  elements  for  slender  bodies, 
t  There  may  be  more  than  four  cards  needed  to  enter  all  modal  information. 
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PROGRAM  NO. 


PROCRAM  NO. 


2.2  Description  of  Input  Data 


CARD 

NO. 

ITEM 

NO. 

MNEMONIC 

SYMBOL 

— 

FIELD 

SOURCE 

DESCRIPTION 

0 

HDR( 15) 

1-60 

15A4 

Header  information 

1 

1 

IBFS 

61-70 

no 

Body  force  calculation 
method  flag, 
alternate  #1  IBFS  =  1; 

alternate  §2  IBFS  =  0. 

* 

2 

FMACH 

M 

1-10 

Mach  number,  usual 

definition 

3 

REFA 

A 

11-20 

Reference  area;  usually 
total  area  of  both  wings 

4 

9EFS 

S 

21-30 

5F10.0 

Reference  semi span 

5 

REFC 

c 

31-40 

Reference  chord; 
usually  average  chord 
of  wing 

2 

6 

XM 

XM 

41-50 

MAIN 

Moment  axis 

7 

ND 

6 

51-52 

12 

Symmetry  flag  (y  =  0 
plane) 

6  =  1  for  symmetry 

6  =  -1  for  antisymmetry 
5=0  for  no  symmetry 

8 

NE 

z 

53-54 

12 

Second  symmetry  flag 
(z  =  0  plane) 
e  =  1  for  biplane 
effect  (symmetry) 

*Use  #1  for  circular  and  #2  for  elliptic  cross-sections 
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■ 

ITEM 

NO. 

MNEMONIC 

SYMBOL 

mm 

■ 

FIELD 

SOURCE 

DESCRIPTION 

e  =  -1  for  ground 
effect  (antisymmetry) 
e  =  0  no  symmetry 

9 

NP, 

55-56 

12 

Total  number  of  panels 

NOPAN 

t 

on  all  lifting  surfaces 

10 

NB 

57-58 

12 

i 

Total  number  of  bodies 

11 

NK 

59-60 

12 

i 

Total  number  of  reduced 

frequencies;  max.  6 

per  case 

12 

MK1 

61-63 

13 

Sequence  number  of 

first  box  on  first 

panel  representing  a 
body  surface,  whenever 
this  body  is  at  zero 
incidence;  otherwise 

MK1  =  0 

13 

MK2 

64-66 

13 

Sequence  number  of  last 
box  on  last  panel 
representing  a  body 
surface,  whenever  this 
body  is  at  zero  inci¬ 
dence;  otherwise 

MK2  =  0. 

Note  that  panels  on 
body  surfaces  need  not 
be  input  last 

2 

14 

Nl, 

67 

11 

MAIN 

Print  flag  for  solu- 

NPR1 

tions;  Nl  =  1  means 

all  solutions  are 
printed,  Nl  =  0  means 
no  print.  Usual 
setting  is  Nl  =0 
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FIELD  SOURCE 


DESCRIPTION 


Control  flag  for 
pressures  and  general¬ 
ized  forces.  N2  =  0 
means  no  pressures  and 
no  generalized  forces; 
N2  =  1  means  pressures 
printed,  and  general¬ 
ized  forces  computed 
according  to  AGARD 
definition;  N2  =  2 
means  pressures  and 
conventional  general¬ 
ized  forces  -  see 
Sec.  5.11.1 

Data  flag;  N3  =  1  means 
DT  matrix  print;  N3  =  0 
means  no  print. 

Usually  N3  =  0 

Detail  print  flag  for 
subroutines  RDM0DE,  SB, 
WANDWT  and  BFM; 

N4  =  0  means  no  print; 
N4  =  1  means  detail 
print  in  subroutine 
BFM  only; 

N4  =  2  means  detail 
print  in  all  four  sub¬ 
routines.  Usually  N4  = 


Array  of  reduced 
frequencies 

Kr  2F 


ITEM  MNEMONIC  SYMBOL  CARD  FIELD  SOURCE 
NO.  COLUMN 


6F10.0 


DESCRIPTION 


Panel  edge  coordinates 

( X-| , Y-j  ,Z-j )  -  inboard 
leading  edge 
( X2 ’ Y-j  >Z-j)  -  inboard 
trailing  edge 
(X3,Y2,Z2)  -  outboard 
leading  edge 
(X4,Y2,Z2)  -  outboard 
trailing  edge 


2F10.0 


Y1  Zl' 


27  NC 


28  NS 


29  NAB(IO) 


31  TAU( 50) 


nc  21-30 


ns  31-40 


41-60 


ei  1-60 


Number  of 
chordwise 
Number  of 
spanwise 


divisions 

for 

panel 


Associated  bodies;  a 
max.  of  six  per  panel 


Fractional  chordwise 
divisions  for  panel. 
Usually  varies  from  0 
at  leading  edge  to  1.0 
at  trailing  edge 


Fractional  spanwise 
divisions  for  panel. 
Usually  varies  from  0 
at  inboard  edge  to  1.0 
at  outboard  edge  ( 

Repeat  Items  #19  through 
31  for  all  panels 
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CARD 

ITEM 

MNEMONIC 

32 

ZC 

33 

YC 

34 

RAD 

35 

AR 

8 

36 

NBE 

37 

NSBE 

38 

NZY 

39 

NRI 

40 

NRS 

FIELD  SOURCE 


DESCRIPTION 


C 


1-10 

11-20 

21-30 

31-40 


41-43 

44-46 


47-48  12 


51-52 


4F10.0 


49-50  12 


z-coordinate 

y-coordinate 


of  body 
axis 


Average  characteristic 
semi -width  of  body 

Cross-sectional  aspect 
ratio  of  body 


Number  of 

interference  body 
Number  of  elements 
slender 


NZY=1  -  z-oriented  body 
NZY=2  -  z-and  y-oriented 
NZY=3  -  y-oriented  body 
Input  bodies  in  this 
order,  i.e.  z-oriented 
bodies  first,  then  z- 
and  y-,  then  y-  bodies 

Interference  'radius1 
flag;  NRI  =  1  -  RI  - 
array  is  input  (see 
below);  NRI  =  0  - 
RI -array  is  not  input, 
but  rather  is  taken 
as  'a',  i  .e.  RI^=a 
for  all  i  =  1,  NBE 

Slender  body  'radius' 
flag;  NRS=1  -  RS-array 
is  input  (see  below) 
NRS=0  -  RS-array  is  not 
input;  instead,  RS..=a, 
all  i  =  1,  NSBE 
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■ 

ITEM 

NO. 

MNEMONIC 

SYMBOL 

— 

FIELD 

SOURCE 

DESCRIPTION 

41 

NSH 

53-54 

12 

DATA 

Number  of  An -As  pairs 
for  body  -  see  Items 

#51  -  53 

8 

42 

NT! 

Nei 

55-56 

12 

Number  of  elements  in 

the  ei  -array 

-  see 

item  #48 

1 

43 

NT2 

N02 

57-58 

12 

Number  of  elements  in 

the  02  -array 
1! 

-  see 

item  #49.  Note  that, 

if  NT2=0,  the  02  - 
v 

array  is  not  input 

9 

44 

1-60 

6F10.0 

x-coordi nates 

of 

inter- 

ference* 

10 

45 

RI(IOO) 

Rli 

1-60 

6F10.0 

Average  char- 

body 

element 

acteristic 

endpoint; 

semi -wi dths 

of 

i=l, 

(NBE+1 ) 

Omit  RI  if  NRI  =  0 

11 

46 

XIS(IOO) 

?Si 

1-60 

6F10.0 

x-coordi nates 

Slender 

body 

element 

endpoints; 

1+1, 

(NSBE+1 ) 

12 

47 

RS(IOO) 

RSi 

1-60 

6F10.0 

Oi 

Average  char¬ 
acteristic 

semi -widths 

of 

Omit  RS  if  NRS  =  0 

*  Omit  Items  #44  and  45  if  NBE  =  0 
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CARO  ITEM  MNEMONIC  SYMBOL  C 
NO.  NO.  C 


48  THl  (24) 


49  TH2 (24) 


LI 

L2 

50  L3 


FIELD 

SOURCE 

6F10.0 

6F10.0 

DATA 

6110 

L4=18  L3=13 


l2=io 


DESCRIPTION 


Angular  orientation  of 
the  points  V  on 
interference  body 
surfaces  -  first 
array,  y  =  1,  NT1 


Second  array  of  eu's 
for  interference 
bodies;  p  =  1 ,  NT2 
Omit  this  item  if 
NT2  =  0 


First,  and 

Last  elements  for  inter¬ 
ference  body  with 
e  ^-distribution; 
a  max.  of  three 
pairs  per  body 


•Interference  Elements 

Ly 


20^19-18  •u.iyg  12  11  JO- 
62  61  02 


4  3  2-"  1  i 


Distributions 


62  Distribution 


Distribution 


Items  #48  and  #49 
are  input  in  degrees 

Omit  Items  #48-50  if 
NBE  =  0 
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ITEM  MNEMONIC  SYMBOL  C 


51  NCD1 


52  CD2 


53  CD3 


An.j  11-20 


An^  41-50 
A?i+1  51-60 


FIELD  SOURCE 


1-10  no 


2F10.0 


A?.  21-30 

31-40  no 


2F10.0 


DESCRIPTION 


Panel  member  for  An> 


Ac  pair 
y-shift 


DATA  z  shi  ft 


of  panel,  first 
set 


Another  set  of  the 
above  three  items. 

Repeat  for  all  i  =  1,  NSH 


ta 


a  _ 


u 


Omit  Items  #51-53  if 
NSH  =  0 

Repeat  Items  #32  through 
53  for  all  bodies;  omit 
same  if  there  are  no 
bodies  for  case. 
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FIELD  SOURCE 


DESCRIPTION 


Punch  * AZ '  in  cc  1-2  of 
all  modal  data  cards  for 
z-oriented  bodies 


Mode  number 
for  the 


first 

set 

second 

of 

third 

modal 

data 

first 

second 

value 

thi  rd  I 

Power  of  x/c  first 

in  the  mode  second  value  I 

polynomial  third 

Coefficient  first 
of  (x/c)n  in  second  value 
the  mode  third 
polynomial 


Punch  'AY'  in  cc  1-2  of 
all  modal  data  cards  for 
y-oriented  bodies 


As  Items  #55-58, 

but  now  for  y-oriented 

bodies 


Omit  cards  #17-18  if  NB=0 
Omit  card  #17  if  NBZ=0 
Omit  card  #18  if  NBY=0 

Punch  'A'  in  cc  1  of  all 
modal  data  cards  for  panels 


ITEM 

NO. 

MNEMONIC 

SYMBOL 

pffjM 

FIELD 

SOURCE 

DESCRIPTION 

3-4 

first 

set 

65 

NQ 

q 

23-24 

12 

Mode  number 

second 

of 

43-44 

for  the 

third 

modal 

j 

data 

5-7 

66 

NRP 

P 

25-27 

12 

Panel  number  -  3  sets 

45-47 

8 

Power  of  y/c  in  the 

67 

M 

m 

28 

11 

mode  polynomial,  where 

48 

y  is  a  'spanwise* 

coordinate  along  the 

lifting  surface,  e. 

9*. 

for  fins  y  = 

=  z 

68 

N 

n 

9 

11 

RDMODE 

Power  of  x/c  in  the 

29 

mode  polynomial 

49 

10 

Flag  that  sets  the 

69 

N8 

N8 

30 

11 

origin  of  the  spanwise 

50 

coordinate  y. 

N8=0  means  origin  of 

coordinates; 

N8=l  means  inboard 

edge  of  panel 

See  Sec.  5.4.1  (Sub- 

routine  RDMODE )for 

i 

details. 

11-20 

1 

Coefficient  of 

70 

31-40 

F10.0! 

(x/c)n(y/c)m  in  the 

a(p) 

qmn 

51-60 

mode  polynomial 
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■ 

ITEM 

NO. 

MNEMONIC 

SYMBOL 

IBB 

FIELD 

SOURCE 

DESCRIPTION 

20 

71 

-1 

-1 

1-2 

12 

RDMODE 

Punch  '-V  in  cc  1-2  of 

card  following  last 

modal  data  card  for  case 

Note  that  only  the  nonzero  coefficients  need  to  be  input;  each  of  the  cards 
#17,  18  and  19  accommodate  3  sets  of  modal  data. 

If  NK>1,  the  modal  data  cards  #17-20  have  to  be  repeated  for  each  reduced 
frequency,  i.e.,  input  these  NK  times. 

2.3  Example  Case  Input  Sheets 

The  following  case  is  to  be  viewed  only  as  an  example  of  the  proper 
input  procedure  and  not  as  an  optimum  idealization.  The  configuration  is 
shown  in  Sketch  2.3-1. 
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The  configuration  consists  of  3  panels  and  2  bodies.  The  panels  are 
given  as  follows: 


O) 

Xi 

3 

2.0  X2 

= 

4.0 

x3 

3 

2.5 

X4 

3 

4.0 

Yi 

= 

0.86603  Y2 

= 

2.0 

2i 

= 

0.5 

lz 

= 

0.5 

0 

= 

0.0,  0.5,  1.0 

and 

t 

3 

0.0, 

0.5,  1.0 

Associated  body  is  body  #1. 

(2) 

Xi 

= 

2.5  X2 

a 

4.0 

x3 

3 

3.0 

x4 

3 

4.0 

Yi 

3 

2.0  Y2 

a 

3.0 

Zi 

3 

.05 

Z2 

3 

1.0 

e 

3 

0.0,  0.5,  1.0 

and 

T 

= 

0.0, 

0.5,  1.0 

Panel 

#2  has  no  associated  bodies. 

(3) 

Xi 

= 

2.5  X2 

a 

4.0 

X3 

3 

2.5 

X4 

= 

4.0 

Yi 

3 

2.0  Y2 

at 

2.0 

Zi 

3 

0.5 

z2 

= 

0.0 

6 

3 

0.0,  0.5,  1.0 

and 

T 

3 

0.0, 

1.0 

Associated  body  is  body  #2. 

The  body  data  is  as  follows: 

Body  #1  is  z-oriented  with  Zc  3  0.0,  Yc  =  0.0  and  is  divided  into  3 
interference  elements  and  5  slender  body  elements  with  end  points 


ei 

li 

o 

o 

«# 

ro 

• 

o 

• 

o 

* 

6.0 

es 

3  0.0, 

* 

o 

« 

C\J 

o 

• 

4.0, 

The  interference  body  radius  is  a  3  1.0,  and  the  array  of  the  slender  body 
element  radii  is 

RS  3  0.0,  0.5,  1.0,  1.0,  0.5,  0.0. 

The  second  element  of  the  interference  body  has  0i  distribution  on  its  sur¬ 
face,  while  the  other  elements  have  e2  distribution: 


Qi 

=  0, 

60, 

120, 

180,  240,  300 

e2 

3  0, 

90, 

180, 

270. 
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Body  #2  is  both  z-  and  y-oriented  with  Zc  =  -0.5  and  Yc  =  2.0.  The 
element  endpoints  are  defined  by 

51  =  2.0,  2.5,  3.25,  4.0,  4.5 

£S  =  2.0,  2.25,  2.5,  3.25,  4.0 

The  interference  body  radius  is  0.5  and  the  slender  body  radii  are  given  by 
the  array 

RS  =  0.0,  0.25,  0.5,  0.0 

All  interference  elements  of  body  § 2  have  0!  distribution: 

01  =  45,  135,  225,  315. 


The  configuration  is  given  three  modes  of  motion: 


Mode  #1  is  a  plunging  motion.  The  coefficients  of  motion  are  given 
by  aij^  =  cos  y  for  all  panels  p,  where  y  =  dihedral  angle  of  panel,  and 

az^b'  =  1.0  for  bodies  with  z-orientation. 


Mode  #2  is  a  pitching  motion.  The  coefficients  are  a^  =  cos  y  for 


,b? 

all  panels  p,  and  az}  ;  =  1.0  for  the  z-bodies. 


Mode  #3  is  a  rolling  motion.  The  coefficients  are  a^  =  (Y-j  cos  y 

+  Z,  sin  y)/c  and  =  1.0  for  all  panels  p.  Body  #1  has  no  motion  in  this 

u  (2)  -  (2) 
mode  while  body  #2  has  two  modal  coefficients  aZg  '  =  2.0/c  and  ayQ  '  =  0.5/c. 


The  input  sheets  are  shown  on  the  next  several  pages,  while  the  out 
put  for  the  case  is  found  on  the  subsequent  pages. 
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PROGRAM  NO. 


PROGRAM  NO, 


PROGRAM  NO. 


2.4  Example  Case  Output 


N5<A  EXAMPLE  CASE 


•  •  ARRAY  Oh  fttOIJCEO  FREQllfcNC  1 E  S  *  * 


0.50000^ 


pErfftCsci  chord  » 
affCa^NCC  SFMl-SPAN  • 
REFERENCE  Aft  t  A  » 

HACK  NU^HFR  « 
h)MfNl  AXIS  ■ 

SYMMtTKY  FLAGS  «*  OK 
S’JHftEft  'F  PANELS  •  3 

MM'JFft  Of  ODOIFS  •  2 


l .Sf^oo 
j .o^ooc 
o  ,4009C 
0.  *15000 
2.00000 

A  >  1  CPSIl’lN  *  0 


••  PANEL  V'.  1  INPUT  VALUES 


xl  • 

2.0LOOOO 

<2  • 

4. CLOD DO 

Y 1  * 

0.  56'»v3P 

n  - 

O.SOOOoO 

x3  » 

2.500000 

X4  * 

A .000000 

Y2  « 

2.000)00 

1 2  ■ 

0.500000 

NC  * 

2  NS  ■  2 

ASSOC! 

AT  FI)  bODIFS 

»  1 

0  0 

0 

0  ) 

J  OrikOMSC  DIVISIONS  F  )«  PANFL  l 


0.0 

0 • SO^OOOOOF  00  0.10000000*  01 

1  SPANWISF  DIVISIONS  FOP  p AUEl 

l 

0.0 

v. sotooopof  o')  ). loojooeot  u 

«•  PAN«l.  '.3.  Z  INPUT  VAIUIS  •• 

XI  » 

2.4A0COO 

*2  .  A.OOO’PO  y  |  .  ?.  000003  Zl  « 

0.500000 

X)  • 

3.0P0CPO 

XA  «  A.OOOOCO  VZ  .  3. 000000  ZP  « 

l. 000000 

NC  * 

y  NS  «  2 

ASSOCIATED  60*)!  fS 

■  0  0  0  0  0  0 

3  CHOP  Owl  SC  DIVISIONS  F'JP  PANEl 

2 

0.0 

O.SOrtCOOOOE  OF-  O.IOoOOOOOC  01 

3  SPANWISF  01 VI  SI  INS  F  )P  PANtl 

2 

r'.o  O.S'-'O  >OOCOF  r0  0.nOC99''OE  01 

«*  PANEL  M*.  1  INPUT  VAL  IF  S  ♦* 


XI  * 

2.500*0 

X2  * 

4.000 OCO  Y l  » 

2.000000  Zl  * 

0 • 5 ''0000 

X)  - 

2.50000) 

X4  * 

4. 0000.00  Y2  « 

2.0CO00O  Z?  ■ 

o.p 

NC  * 

7  MS  •  I 

ASSOC  1 

iatfo  onoirs 

«  2 

0  0  0 

0  ) 

3  CH0ft04»ISf 

DIVISIONS  f-  PANEL 

3 

O.o 

9.5C9COOCOF  r 

0  O.IDOC-OOE  PI 

7  v>AMwI$f  9! VISIONS  PANEL  3 


0.  l.>orOOAOC  01 
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**  summary  of  pancl  data  ** 


PANfl  NC  NS  N8-ARRAY  0  IMF  ORAL  NO.  Of  ASSOC. 

anuu  o»io  i  e  s 


LIST  OF  ASSOCIATED  RUDIES 


4 

8 

10 


O.o 

26.56505 

270.00000 


•*  GEOMETRY  ARRAYS  FOR  ALL  PANFLS  ♦* 


PANCI 

NO. 

•STRIP 

NO. 

BOX 

NO. 

3/4  Chord 

X 

1/4  CHORD  X-CDOROINATFS 

IN HOARD  CENTFR  OUTBOARD 

BOX  CHORO 
QFLTA-X 

1/4  CHORO 

sweep  anglf 

1 

1 

1 

2.B2813 

2.25000 

2.35938 

2.46875 

0.93760 

0.36022 

1 

I 

2 

3.76563 

3.25000 

3.296d8 

3.34375 

0.93750 

0.16307 

I 

2 

3 

2.98438 

2.46875 

2.57813 

2.68750 

0.81250 

0.36822 

l 

2 

4 

3.79688 

3.34375 

3.  39063 

3.43750 

0.8125C 

C. 16387 

2 

3 

5 

3.14063 

2.68750 

2.79688 

2.90625 

0.68750 

0.37299 

? 

3 

6 

3.82813 

3.43750 

3.48433 

3.53125 

0 .68750 

0. 16616 

2 

4 

7 

3. 2968 H 

2.90626 

3.01563 

3.12500 

0.56250 

0.37299 

2 

4 

8 

3.85930 

3.53125 

3.5781 3 

3.62500 

0.56250 

0.16616 

3 

5 

9 

3.0625C 

2.68750 

2.68759 

2.68750 

0.75000 

0.0 

3 

5 

10 

3.81250 

3.43750 

3.43750 

3.43750 

0.750PP 

0.0 

NO. 

Y 

l 

OCITA-Y 

OfLTA-2 

F 

CHORD 

X-L.F. 

1 

1.14952 

0.50000 

0.56698 

0.0 

0.28349 

1.87500 

2.12609 

2 

1.71651 

0.50000 

0.56698 

0.0 

0.28349 

1.62500 

2.37500 

3 

2,25000 

0.62500 

0.60000 

0.25000 

0.27961 

1.37500 

2.62600 

4 

2. 75000 

0.87500 

0.50000 

0.25000 

0.27951 

1.12500 

2.87609 

5 

2.00000 

0.25000 

0.0 

-0.50C00 

0.25000 

1.50000 

2.50C00 
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**  BUOY  NO 


1  INPUT  VALUES 


CENTER  Of  BODY  COORDINATES  v  *  0.0  l  *  0.0 

AVERAGE  HALF-WIDTH  OF  BODY  *  1.000000 

CROSS-S tCT l ON A l  ASPECT  PATIO  *  l .000000 

NUMBER  OF  INTERFERENCE  ELEMENTS  ON  30DY  »  3 

NUM3FR  OF  SlFNOFR  BUOY  ELEMENTS  "  5 

2-Y  ORIENTATION  FLAG  »  l 

RI -FLAG  *  0  «-S  CL AG  •  1 

number  or  oclta-eta  OFLTA-ZFTA  'AIRS  *  0 

4  X  I- 1  ELEMENTS  FOR  BODY  1 

0.0  0. 20000000C  01  0.40000000E  Cl  0.60000000E  01 

4  R-l  ELEMENTS  FOR  BODY  1 

0.  lOOCOOOOF.  01  0. 10000000E  01  ()•  lOOOOOOOf  01  0 •  10000000E  01 

6  Xl-S  ELEMENTS  FOR  BODY  1 

c,o  o.icooooooF  ci  o. aooocoooe  01  0.40000000E  01  o.sdoooccoe  01  c.<.coooo.>oe  "i 

6  R-S  ELEMENTS  FOR  BODY  l 

0.0  O.SOOOOOOOL  00  O.lOOODOOOf  01  0 •  10300000E  0)  0.tf»00CCft«C  00  f>. 0 

4  THETA- I  ELEMENTS  for  BODY  l 

C.o  0. 60000000E  02  0.12000000r  03  O.HOOOOOOE  03  C.24000COCP  01  0.  .'OOOOOOOt  03 

4  THETA-2  ELEMFNTS  FOR  BOOV  l 

0,0  0.90000000E  02  0. 18000000F  03  0.270000006  03 

THE  FIRST  ANO  LAST  BODY  ELEMENTS  FOR  THETA- 1  (IN  BODY  1 

2  2  0  0  0  0 
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**  BODY  NU.  2  INPUT  VALUES  *• 

CENTER  Of  BODY  COORDINATES  Y  *  2.00C0O0  l  *  -0. *>00000 

average  half-width  OF  BODY  *  O.SOOOOC 

CKOSS-SCCriONAl  ASPbCT  RATIO  »  l. 000000 

NUMRfcR  or  INTERFERENCE  ELEMENTS  ON  BODY  a  4 

NUMBER  OF  SLENDER  BODY  ELEMENTS  «  4 

/-Y  ORIENTATION  HAG  -  2 

P I -FLAG  *  0  R-S  FLAG  •  I 

number  OF  OFLIA-ETA  OCL T A— / C T A  PAIRS  «  0 

5  X  I - 1  ELEMENTS  FOR  BODY  ? 

0.  200000C0C  01  0. 2SBCOOrtOF  01  0. 325COOOOE  01  0.40000000E  01  0.4SOO''OOOF  01 

5  ft- I  ELEMENTS  FOR  BUOY  2 

0.50OO0000F  00  O.5O000DU0E  CO  O.SOOCOOOOF  00  0.50000000C  00  O.SOOOOOOOE  00 

5  Xl-S  ELEMENTS  FOR  BODY  2 

0.20000POOF  01  0 •  225000C0E  Cl  0.250000001  01  0.325000006  01  C .4r>OC(/OOOE  Cl 

S  ft-S  ELEMENTS  FOR  BODY  2 

0.0  0.250000006  00  0.50000000E  00  0.500G00G0E  OC  O.C 

4  THETA-1  ELEMENTS  FOR  300Y  2 

0. 45000000E  02  0.13500000F  03  0.22500000F  03  O.31SO0000F  03 

THE  FIRST  AND  LAST  BODY  CLEMENT S  FOR  THETA- l  ON  BODY  2 

1  4  0  0  0  0 
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***** 

EXECUTION 

CPU 

I/O 

TOTAL 

***** 


EXECUTION 

CPU 

I/O 

TOTAL 

***** 


EXECUTION 

CPU 

I/O 

TOTAL 


**  RECEIVING  POINT  ARRAYS  X,  Y,  l,  GAMMA  ** 


0.202812506 

01 

0.376562506 

01 

0.298437506 

01 

0. 379687506 

91 

0.314062506 

Cl 

0. 382812  50F 

01 

0. 329687506 

01 

0.385937506 

ot 

0.306253006 

01 

0.301 25000C 

01 

0. 190000006 

01 

0 • 3000Q09CF 

01 

0.500000006 

Cl 

0.225000006 

01 

3.2875COOC6 

01 

0.362500006 

o  l 

0.4  25COOOOE 

01 

0. 1 1496? 106 

01 

0.171650/06 

Cl 

0.225000006 

<H 

0.27500000F 

*1 

O.2CCC0cnftg 

01 

0.0 

0.200000006 

01 

0. 5  OOCOOOQF 

00 

0. 50000000E 

00 

0.6250f000fc 

00 

O.075OOOOOF 

OC 

C.25CC00096 

oc 

0.0 

-0.500000006 

00 

0.0 

0.0 

0.46364/436 

00 

0.46364743F 

00 

0.47 1 21 d  7| 6 

01 

0.0 

0.0 

**  LENDING  POINT 

ARRAYS  XII*  X12,  A0» 

AO-PRIME  ♦* 

0.0 

0.100000006 

01 

o.2ooooooor 

01 

C. 400000006 

01 

0 .590COOCL  6 

01 

0.2000C00CF 

Cl 

*.225000006 

01 

0. 2500000CF 

01 

0.325000006 

01 

0. 10009000E 

01 

0, 200000C0E 

01 

0, 4Q0000C0E 

01 

0.5000000CC 

01 

0.600CGOOGF 

01 

0. 225900006 

Cl 

0. 25000000E 

01 

0.325C00006 

01 

0.45000000F 

01 

0.25000000F 

00 

0.760000006 

CO 

0,  lOOOOOf'OE 

01 

0 • 750009001 

90 

9.2500*0*0t 

GO 

0. 12500000F 

00 

0.375000006 

oc 

0. 50000000C 

00 

0.250000006 

00 

Q.6COOOCOOE 

00 

0.5000^0006 

00 

0.0 

-0 .600009006 

CO 

-0. 5*0f  0*r0L 

90 

0. IC00C09CF 

01 

o. loooroooa 

01 

0.0 

-0.39999998F 

00 

TIME(MINUTES) 

0.0067 

0.0186 

0.0253 


T I M6 I M I NOTES ) 
0.0 
0.0 
0.0 


THEIMINUTCSI 
0.3162 
0.0351 
0. 1512 
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INPUT  MOOAL  OAT  A  FOR  —  A  — 

1001000  1.0000  100200C  O.Q660  2001010  1.0000  2002010  0.8660 
300IOCI  0.6774  3001101  1.0000  1002001  1.3214  3002101  1.0(00 
3003001  -0.3333  3003101  1.0000 


INPUT  MOOAl  OATA  FOR  —  A /  ~ 

10010OO  l .0000  1002000  1.0000  2C011C0  1.0000  2002100  1.0000 

3002000  1.6667 


INPUT  MODAL  DATA  FOR  --  AY 
30C2000  0.3333 


CXFCUT1CN  T iME(MlNUreS) 
CPU  0.0000 

I/O  0.0062 

TOTAL  0.0062 

***** 


EXECUTION  T  I  (MINUTES) 
CPU  0.0 

I/O  0.0 

total  0.0 


CALCULATION  "F  02  AND  OY.  N  -  2 1  NM  *  3  NWS  *  9  NTYS  «  4 

MAX.  COftf  N66DCO  ■  320  C0R6  AVAIL.  =  10O00 

ALL  CALCULATIONS  DONE  IN  CORC. 


MODE  -  1 
NT/  S  -  9 
NTYS  -  4 


U-Z,  U-Y, 

CP-Z*OFLTA-A 

,  CP-Y*06L  T A-A 

0.0 

0.78546  00 

o.o 

0.7*696 

01 

0.0 

0.12576 

02 

0.0 

0.7069E  01 

0.0 

0.78546 

00 

0.0 

0.1963F 

OC 

0.0 

0.17o7E  01 

0.0 

0.3142E 

01 

0.0 

0.78546 

OO 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

-0.6236E 

00  0. 3 1426 

01 

-0.47126 

01 

C.9426C 

01 

-0.1676E  02 

0.0 

-0.4  71 2f 

01  -0.94266 

01 

-0.52366 

00 

-0.31426 

01 

-0.3272001 

0.7854C  °C 

-0.29456 

00  0.2356F 

01 

-0* 15  7  IF 

01 

C.O 

-0.6545F  00 

-0.3 1421  01 

0.0 

0.0 

0.0 

c.o 

0.0 

0.0 

0.0 

0.0 

MODE  - 

2 

ntzs  - 

9 

NTYS  - 

4 

U-Z,  U-Y, 

CP-Z  *OEL  TA-A 

,  CP-Y*DFL T A-A 

0.7854P  00 

0.26I8E  00 

0.70696 

01  0.70696 

01 

0. 12576 

0? 

0.2513F 

02 

0.7(696  01 

0.2  1  2 IF  02 

0.78546 

OC  0.2880E 

01 

0.19636 

OC 

0.2782F 

OC 

0.1767E  01 

0.27986  01 

0.31426 

01  0.602  If 

01 

0.78546 

OC 

o.2C29r 

01 

0.0 

0.0 

0.0 

0.0 

o.o 

0.0 

0.0 

0.0 

0.2967E 

01  0.20946 

01 

0.47126 

01 

0.18856 

02 

-0. 3361 C  02 

0.3351F  C2 

-0.23566 

02  -0.18856 

02 

-0.50516 

01 

-0.10476 

02 

0. 7390F  00 

0.11786  01 

0.18906 

01  0.4320F 

01 

-0.301  IF 

J1 

0.31426 

01 

-0.4832F  01 

-0.6807F  01 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

9.0 

MODE  - 

3 

NTZS  - 

9 

NTYS  - 

4 

U-Z,  U-Y, 

CP-Z  *OELT  A-A 

,  CP-Y  *PEL  T  A- A 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

c.o 

0.0 

0.0 

0.32726 

00 

0.0 

0.29466  01 

0.0 

0.52366 

01 

0.0 

C. 13096 

01 

0.0 

C. 65456-01 

0.0 

0.5H90F 

OC 

0.0 

0.1047F 

01 

0.0 

0.26186  CO 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

-0.54546-01 

0.13096  01 

-0.49P9L 

00  0.39276 

01 

-0, 261  86 

01 

0.0 

-0. I0°lf  01 

-C.6236E  01 

-0 . 109  IF- 01  0.2618E 

oo 

-0.90176* 

■01 

0.78546 

00 

-0.5236F  00 

0.0 

-0.21821 

00  -0.1C47E 

01 
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COLUMN  l  l/f- 

1  -0.1676$? 

3  -0.170*8! 

8  -0.16'76E 
7  -c.m«*F 

o  o,? 14746 

U  C. 75136b- 
13  -0.  i«596F- 
15  -0.10968C 
P  -0. 1 70356 
19  -0.13416E-' 
21  -0.914786 


ObLTA-  W- 
*0  -0.2f>7?56  00 
CO  3.97032E-01 
00  0. 19870E-01 

00  -0.©3007F-02 
Cl  r.$65B4C-01 
*2  -0. 13743C-01 
3*  -0.381C8E-01 
00  -0.  106A3F  00 
00  -0. 1 8190F  00 
01  -0.  10852b  00 
02  -0.13222E  00 


2  -0.27663b  00  -0.40265C  00 
4  -0.227291  00  -0.17359b  00 
6  -0.»33C6f  00  -0.10716b  CO 
Q  -0.  *.96211  00  600316-01 


10  0.277M3V--01  0  752236-01 

12  -0.7/854E-01  -0. 406*96-0: 

19  -0.W996E  00  -O.P**'76  CO 
16  -0.17229b  00  -0. 117915  00 
18  -0.20927F-01  -0.12301b  00 

20  -0.997639-01  -0. 12259b  00 


COLUMN  2  Ob  OfclTA-  W- 


1  -0.5O579E  00 
?  -0.26295F  00 
5  -0.27637F  00 
7  -0.2231 56  00 
9  0.10460F  00 

It  -0.19271C-02 
13  -0.99692F-01 
15  -0. 30299T  CO 
17  -0.60529F  00 
19  -0.197306  00 
21  -0.  196646  C*) 


-0.29402F  00 
0.32520E  00 
0.25650E  00 
0.15608E  OC 
6.80B71E-01 
-0. 33497F-01 
-0. 81319E-C1 
-0 • 5757  IE- 01 
-0 •  16  2536  00 
-0.16584b  00 
-0.29960F  00 


2  -0.106816  01  -0.450376  00 
4  -0.69030F  00  -0.86960F-02 
6  -0.509486  00  0.58633F-01 
B  -0.338296  00  C.1183C6  00 
10  0.15612F  00  0.125436  00 
12  -O.IB143F  00  0. 563486-02 
14  -0.337 Jbf  00  0.27377b-02 
16  -0.47061 C  00  0.l4893b-01 
18  -0.18410E  00  -0.184976  00 
20  -0.24727E  00  -0.19016E  00 


COLUMN  3  OF  fill  TA-  H- 

1  -0.77452F-01 

C. 703026-01 

2  -0. 678306-01 

-0.5  7944E-01 

3  -0.920546-0. 

0.3217CF  00 

4  -0.64814F-01 

0.2 199JF-01 

5  -0.84912E-01 

0.23316E  00 

6  -0.42804F-01 

0. ©98846-01 

7  -0.52960F-01 

0.70497F-01 

8  -0. 377146-01 

0.29099E-CI 

9  0.81054E-02 

-0.13012F  00 

10  -0.1 5097E-01 

-O.U022E  00 

U  0.12523F-01 

-0 • 22905C-01 

12  -0.12142C  00 

-0. 676996-01 

13  -0.326606-02 

-0.63509E-CI 

14  0.12716E-01 

0.968486-02 

15  -0.29452E-G2 

0.189966-01 

ib  -O.17280E-01 

0*  10752E-01 

17  -0.19371E-0! 

0. 14849F-02 

18  0.132096-01 

0. 397836-02 

19  0.1121 lc-01 

21  0. 9  7084E-09 

0.64860E-02 

C. 174236-01 

20  0.13454E-Q1 

0. 1083 lE-Ol 

EXECUTION  r I MC  f  M I  NOTES) 

CPU  0.0754 

I/O  0.0404 

TOTAL  0.1158 

ABSOLUTE  COORDINATES  USFO  FOR  PANEL  l. 

ABSOLUTE  COORDINATES  USED  FOR  PANEL  2. 


W  FOR  MODE  I 


0.0 

0.1000F 

01 

0.0 

0.10006 

01 

0.0 

0.10006  0! 

0.0 

O.IOOOC 

01 

0.0 

0.86606 

00 

0.0 

0. 96©06  00 

0.0 

0.86606 

00 

0.0 

A. 86606 

00 

0.0 

0.0 

0.0 

0.0 

— Wl  —  FOR 

MODE  1 

0.16  766  OC 

0.1 2676 

01 

0.27666  00 

0. 1403F 

01 

0.1786b  00 

0.9530E  90 

0.2273S  ^0 

0.1  1  7  46* 

01 

0.162H6  CO 

0.84o5F 

00 

0. 183 IF  00 

0.97326  00 

0.125??  00 

0.87236 

00 

0.15426  OC 

0.9 1 2  IF 

00 

-0.21476-01 

-0.56886-01 

-0.27  78 6*. 01 

-0.76226 

-01 

-0.7514E-02 

0.1379E- 

‘01 

0.72856-01 

A.4C62F-91 

0.19606-02 

0.38  11C* 

-01 

C.1300F  00 

0.10106 

00 

0. 1 0976  00 

0.10646  00 

0.1723F  00 

0.1L74F 

CO 

0.1703F  00 

0. 18196 

00 

0. 2043F-01 

0.12306  00 

0.1 >9 '6-01 

0.I085F 

00 

0.4476E-01 

0.1225? 

00 

0.  <*1686- 02 

0.1322E  00 

ABECture  coordinates 

USEO  FOR  PANEL 

1. 

A3S0LUTC  COORDINATES 

US60  FOR  PANEL 

2. 
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W  FOR  MODE 


2 


O.lOOOf  01  0.1885E  01  t. 10006  01 

0.10006  01  0.25 3 IF  01  0.8660E  00 

0.86606  00  0. 19036  01  0.8660F  00 

0.0  0.0 

— wr—  FOR  MIJOF  2 

0.15966  01  0.2179C  01  0.20656  01 

0.1690F  01  0.25406  01  0.11426  01 

0.10896  01  0.1747F  01  0.12046  01 

-0. 1561 F  00  -0.1264F  00  0.1927F-02 

0.44646-01  *.01326-01  0.33746  0" 

0.4706F  00  -0.14896-01  0.60536  00 

0.14736  C<*  0.18586  00  0.24736  00 

P6UT1VF  COUROINATCS  USFO  FOR  PANEL 
RCLAI 1 VF  COORDINATES  USED  FOR  PANEL 
RFLATIVF  COORDINATES  USEO  FOR  PANEL 


0.25106  01  0.10006  01  0.1990E  01 
0.18136  01  0.8660F  00  0.22106  01 
0.22286  01  0.0  0.0 


0.2961F  01  0.1282C  01  0.1664E  01 

0.15576  01  0.13766  01  0.21516  01 

0.21106  01  -0.10466  00  -0. 80876-01 
0.3345F-01  0*  1 8146  00  -0.5635F-02 

-0.27380  02  *.30306  30  0.5757F-01 

0.1625F  OO  0.18416  )C  0.1850C  OC 
0.19026  00  0.1968F  K  0.29966  00 

1. 

2. 

3. 


W  FOR  MODE  3 


0.0 

0.7663E  00 

0.0 

0.0 

0.1144F  01 

0.0 

o.o 

0.18806  01 

0.0 

0.0 

-0.16676  00 

--WT--  FOR 

MODE  3 

0.7745F-01 

0.6880E  00 

0.67836-01 

0.64816-01 

*.11226  01 

0.849IE-01 

0.5296F-01 

0.1810F  01 

0. 37  7 1 C- 01 

0.15106-01 

-0.56446-01  -0.12526-01 

0.32666-02 

0.6351F-0 l  - 

•0.12726-01 

0.172RF-01 

-0.1075F-01 

0. 19376-01 

-0.11216-01 

-0.64866-02  -0.13456-01 

*  **«» 

EXECUTION  T 

IMFIMINUTFS) 

CPU 

0.0133 

I/O 

0.0106 

TOTAL 

*.0240 

***** 


EXECUTION  T 1  ME ( Ml  NOTES  I 
CPU  0.0111 

I/O  0.0133 

TOTAL  0.0244 


0.7663F 

00 

0.0 

0.1144E 

01 

0.15C86 

01 

0.0 

0.15086 

01 

0.1880E 

01 

0.0 

-0.16676 

on 

0.8243E  00  0.9205F-01  0.8226F  00 

0.1275E  01  0.42806-01  0.1438E  01 

0.1851E  01  -0. 8105F-02  -0.3655C-*l 
0.22916-01  0.1214F  OC  0.67706-01 

-0.9655E-02  0.2945C-02  -0.19006-01 

-0.  1485E-02  -0.  1 32  IE-01  -0.3978F-02 
-0. 1083F-0 1  -0.9708E-02  -0.17426-01 
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COLUMN  NU 


1  OF  CAHHAS  FOLLOWS 


0.196&56236  01 
-ft* l 7455091 6  01 
0.16948707F  01 
-ft  *76321  1856cftl 
-0.54306471F  Cft 

0.410557006  or 

C.  102442.466  01 


0.92423382*  01 
0*  4937H624F  01 
0.43l4R623f  31 
0.163301  UC  01 
-0.O650265KC  00 
-0,8697093  76  00 
0.947699596  00 


-0. 11037083b  01 
0. 1858071 3F  01 
-0. 1492551 86  01 
0. 573873501-01 
0.  302  7  84  74C  00 
-O.12UO062C  00 
0.  1373  3<*73l  Cl 


0. 462 1 1 bB06  01 
0.436878406  Cl 
0.10070200b  ftl 
0  •  1  9*.  78 1>4  7b  oft 
0.160268106  00 
-0.16611433-  w' 
0. I 5H  77 1 80t  ftl 


0.2603ft857t  01 
-ft.  l>m7rf6MF  ftl 
* .  1  JSO’i^SF  01 
r  .2??6f  00 

-r  .4*»4‘  9-47f  -ft  * 

ft.»4t<.2  33J»  OL 
*•143  1002  7f  01 


0.481331541  01 
0.3Q424414F  01 
0.1C761557C  01 
-r.3300047sr  oi 
-0, 866360471  OC 
0.40218151*  00 
0.00667246b  00 


COLUMN  NO.  2  OF  CAMMAS  FOLLOWS 

0.H3234043R  01  0.5546*4536  01 

0.24164743F  01  0.144840936  02 

0.790010236  01  0.739005456  01 

0.202006376  01  0.39445601*  01 

-0. 24903203b  01  -0.1O5297046  01 

-0.100113o3C  00  -C. 28204*326  01 

0.206001136  01  0.855753716  00 


ft.  364484126  01  0. 1  2030080b  02 

0.34354019b  01  C. 603129306  01 

-'..00620325b  00  0.724732881  01 

0.200872096  CO  0.114335096  Oft 

0.653377126  00  -0 . 1 4427 942fc  OC 

-ft. 31O25070F  00  -0.556054206  00 

0.470909606  Cl  C.10422626F  01 


0.10454  *34C  02  C. 6221  31366  01 

-0.347ftoS$8r  00  0. 10 1  ‘•2693b  0? 

7.3764W07HF  01  0.34082502b  OC 

- ). '44790906  01  -C. 751627161  01 

-C. 723036^96  ftO  -C. 1 22750446  01 

C» 109 1 2 321 C  01  0.290803406  00 

0 .42601 456 F  01  0 .8242081 46  OC 


COLUMN  NO.  3  OF  CAMMAS  FOLLOWS 

0.1 7746496F  01  0.27028K09F  01  -ft.958l627QC  00  0.409359076  01  0.217061606  01  0.417509986  01 

-0.211270006  01  0.474133406  01  0.160891366  01  0.62Po3R466  01  -0. 30868311b  01  0. 36716051b  01 

0.133191786  01  0.770206476  01  -0.30920459F  01  0.241830796  01  0.71831417F  00  -0.4301627^6-01 

0.102434706  00  0.883634006  00  0.856153966-01  0.100280096  00  0.153521026  01  -0.21076336F  01 

-0.386846846  00  -0.39633286F  00  -0.15489656F  00  -0.1939S346F  00  -0.25440876b  OC  -0.Q7 1 3050 7f  Of 

0. 69175780b— 01  -0.13080824C  01  -0.445067026  00  -0.805&20646  00  0.203273006  00  -0.157312616-02 

0.5 1055562F  Oft  0.  270523136  00  0.4711  71026  CO  0.485300426  00  0.414597076  00  0.195830446-01 

THE  21  X  21  MATRIX  WITH  3  RIGHT  SlUfcS  WAS  S0LV60  DIRECTLY  IN  0.005  MlNUTfcS. 

*«**« 

EXECUTION  timeiminutcsi 
CPU  0.0211 

I/O  0.0146 

TOTAL  0.0357 


LOADS  ON  SLENDER  BODY  ELEMENTS  DUE  TO  LIFTING  SURFACE  BOXES 


1 

o.O 

F2 

0.0 

0.0 

FY 

0.0 

0.0 

M2 

0.0 

0.0 

MY 

0.0 

2 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

3 

0 .5731 E 

00  O.7088E  01 

0.0 

0.0 

0.11426 

01  -0.9615C  OO 

0,0 

0." 

4 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

5 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

LOAOS 

ON  SLENDER 

BODY  ELEMENTS 

i  DUE 

TO  lifting 

SURFACF  BOXES  AND 

INTERFERENCE 

BODY 

ELEMENTS. 

FI 

FY 

M2 

MY 

l 

0.28606-01 

0.77396-01 

0.0 

o.o 

o.o 

0.0 

0.0 

0.0 

2 

0.20606-01 

0.7/396-01 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

3 

0.5476F  00 

O.6012E  01 

0.0 

0.0 

0.1246F  01 

-0.94816  09 

0.0 

A,  A 

4 

-0.4089F  00 

0.7298E  CO 

0.0 

o.o 

-0.2773E-01 

C.3652C-01 

c.o 

0.0 

5 

-0.29R0F  00 

0.  5H*BF  00 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

ELEMENT  LOADS 

FOR  SLFNOFR 

BODY 

NUMBER 

t.  MODE  NO. 

1 

F2 

TY 

M2 

MY 

1  -0. 23111 

CO  O.1N40F 

01 

0.0 

o.o 

0.0 

0.0 

0.9 

0.0 

2  -D. 23276 

01  0.47906 

01 

0.0 

o.o 

0.0 

0.0 

0.0 

0.0 

3  -0.8306F 

01  0.70736 

01 

0.0 

0.0 

0.1212F  01 

-0.2*>b  16  00 

O.C 

0.0 

4  -0.2B19F 

01  -0.4243F 

01 

0.0 

0.0 

-v,.  41326-01 

-0. 28706-01 

0.0 

O.o 

5  -0.55986 

00  -0.98706 

oc 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

LOADS 

ON  SLEN06R 

BODY  ELLMENTS  DUE 

TO  LIFTING 

SURFACF  BOXES. 

F2 

FY 

M2 

MY 

l 

0.0 

0.0  0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

2 

O.C 

0.0  0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

3 

0.91936  01 

0. 1 b 1 3F  02  0.0 

0.0 

0.9295F 

CO  -0.44796  01 

0.0 

o.o 

4 

0.0 

0.0  0.0 

0.0 

o.o 

0.0 

0.0 

0.0 

5 

0.0 

0.0  0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

LOADS  ON  SI LNUER 

BODY  ELEMENTS 

OUE 

TO  LIFTING 

F2 

FY 

I  0.I454F  00 

0.57176-01 

0.0 

0.0 

2  0.1454F  00 

0.57176-01 

0.0 

0.0 

3  0.8732F  01 

0*14356  0? 

0.0 

0.0 

4  -0.54086-01 

0. 1 86 IF  01 

O.o 

c.o 

5  -0.10466-01 

0.1391F  01 

0.0 

0.0 

SURFACE  BOXES  AND  INTERFERENCE  BODY  EL6MFNTS. 


m;  My 


0.0 

0.0 

0.0 

O.C 

0.0 

0,0 

0.0 

o.o 

0.1097F  01 

-0.4&02F 

01 

0.0 

c.c 

0. 10916-01 

<5.11  756 

00 

0.0 

0.0 

0.0 

0.0 

0.0 

O.C 

ELEMENT  LOAOS 

FOR  SLCN06R 

BODY 

NUMBER 

1*  MOOc  *M* 

2 

F2 

FY 

M2 

MY 

l  0.1629F 

01  0.11046 

01 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

2  0.25026 

01  0.94826 

01 

0.0 

0.0 

0.0 

0.0 

9.0 

9.0 

3  -0.87046 

01  0.32056 

02 

0.0 

0.0 

0.1774b 

01  -0.32166  01 

C.O 

0.0 

4  -0.12246 

02  -0.81296 

01 

0.0 

0.0 

-0.U12C 

00  -0.23796-01 

0.0 

O.o 

5  -0 .254  IF 

01  -0.3845F 

01 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 
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LOADS  ON  SLENDFR  BOOY  ELEMENTS  OU6  TO  LIFTING  SURFACE  80X6$ 


1 

0.0 

FZ 

0.0 

0.0 

FY 

0.0 

0.0 

MZ 

o.o 

°.o 

MY 

o.o 

2 

0.0 

0.0 

0.0 

o.o 

0.0 

o.c 

0.0 

0.0 

3 

0.1904E 

00  0.63566  01 

o.o 

o.o 

0.1214E 

01  -0.12006  01 

O.C 

o.c 

4 

0.0 

0.9 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

5 

0.0 

0.0 

0.0 

o.o 

0.0 

0.0 

0.0 

C.o 

LOAOS  ON  SL6N0CR 

K0«1Y  eiEMCNTS  0*16 

TO  LIFTING 

SURFACE  90X6$  ANO  INTERFERENCE 

BOOY 

CLEMENTS. 

f  L 

FY 

MZ 

MY 

l  0.42816-01 

0.9964E-01 

0.0 

0.0 

0.0  0.0 

0.0 

0,0 

2  0  *  42  8  IE -01 

0.99646-01 

0.0 

0.0 

0.0  0.0 

0.0 

0.0 

3  0.13356-02 

0,60226  01 

0.0 

0.0 

0.11816  01  -0.12776  01 

0.0 

0.0 

4  -0.58756  00 

0.51366  00 

0.0 

0.0 

-0.2674E-01  0.21446-01 

0.0 

0.0 

5  -0.48056  00 

0.4?78E  00 

0.0 

0.0 

0.0  0.0 

0.0 

0.0 

ELEMENT  LOAOS  FUR 

SL6N06R  e«)OY  NiHBFR 

l.  MODE  NO. 

3 

FZ 

ry 

MZ 

MY 

1  0.428  IE-01 

0.99646-01 

o.o 

0.0 

0.0 

0.0 

o.o 

0.0 

2  0.42816-01 

0.99646-01 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

3  -0.79136  00 

0,64566  01 

0.0 

0.0 

0.1125F 

01  -0.12176  00 

0.0 

0.0 

4  -0.67816  OD 

0,  78766-01 

0.0 

0.0 

-G.iMOF-Ol  -0.8727E-01 

0.0 

0.0 

5  -0.48C56  00 

0,42786  00 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

LOADS  ON  SLEN06R  BODY  ELEMENTS  DUE  TO  LIFTING  SURFACE  BUX6S. 


1  -0.8696E-03 

2  -0.17606-01 

3  -O.I2636  00 

4  0.98306-01 


l 

-0,18766-02 
-0.32096-01 
-0.26816  00 
-0.26496  00 


0.15076-01 

0.45336-01 

o.zonr  oo 

-0.55556-01 


0.32506-01 
0.92176-01 
0.30896  00 
0.3871F  00 


MZ 


MY 


0.54356-04 
0.99106-03 
-0. 28976-02 
0.22216-01 


0.11726 -03  -0 . 94 l 7C -03 
0.17716-02  -C.9499F-03 
0.53006-02  0.20626-01 

-0.67096-01  -0.21936-01 


-0 .203  lk-02 
-0.16996-02 
C.  1436F-01 
°  .92 1 2F-0 1 


LOAOS  ON  SLCNDER  »O0Y  ELEMENTS  OUc  TO  LIFTING  SURFACE  BOXFS  ANO  INTERFERENCE  OODY  ELEMENTS . 


F  2 


l  0.15196 

00 

0.32266-01 

0.18796 

00 

0.23586 

00 

0. 54356-04 

0.11726-03 

2  0.13526 

00 

0.52056-01 

0.21826 

00 

0.29546 

00 

0.99106-03 

0.177U-02 

3  -0.1019E 

00 

-0.53736  00 

0.46176 

00 

0.60716 

00 

-0.289 76-02 

0.53006-02 

4  -0.21406 

00 

0.20376  00 

0.71976- 

-01 

0. 104&F 

00 

0. 11866  00 

-0.18276  00 

0.2062F-01  0.1439E-01 

-0. 6069F-02  U.2040E  00 


ELEMENT  LOAOS  FOR  SL6N0FR 
FZ 

1  0.8094E-01  0.48116 

2  -0. 68  736-01  0.12496 

3  -0.1063E  01  -0.53736 

4  -0.77066  00  -0.15146 


800Y  NUMBER  2.  mode  NO. 


00 

01 

00 

01 


FY 


C. 15896  00 
0.18926  00 
0.3748F  00 
0.47616-01 


0.23586  00 
0.29546  00 
0.60716  00 
0,39356-01 


1 

I 

0.54356-04 
0.99106-03 
-0.28976-02 
0. 1 9076  00 


0.11726-03 
0.17716-02 
0.53006-02 
-0.13696  00 


-0.94176-03 
-0. 94996-03 
0.2C62E-01 
-0.15576-01 


-0.20316-02 
-0.16986-0? 
0.14386-01 
0.2285E  00 
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LOADS  ON  SLENOf*  BUOY  ELEMENTS  OU6  TO  LIFTING  SURFACE  OOXCS. 

F2  F  Y  MJ  MY 

1  -0.36806-0?  -0.24S2C-02  0.63776-01  0.42496-01  0.23C06-03  C. 15336-03  -C.39B6F-02  -0.265'F-O? 

2  -0.6892F-01  -0.4 1536-C1  0.1865F  00  0.1201E  00  0.38476-02  C. 22896-02  -C.3698F-02  -0.2194<„-02 

3  -0,53936  OO  -0.39386  00  D.7?60t  00  0.34896  00  0.83946-03  0.1785E-01  C.5?61F-01  -0.6045F-0? 

4  -0.11376  CO  -0.7R51F  00  0.4096P  00  9.10126  01  -0.34101-01  -0.19426  00  0.6107E-01  0.25576  00 


IUAUS  ON  SLENDER  BODY  ClEMtNTS  DUE  TO  LIFTING  SURFACE  BOXES  AND  INTERFERENCE  BODY  H6MFN1S. 

T2  FY  M2  MY 

1  9.3216F  AD  -0.79076-01  0.6104*  CO  0.18846  00  0.23006-03  0.I533E-03  -C.3986F-02  -0.26561-02 

2  0.25H4F  03  -0.1182P  00  0.7332E  00  0.266CE  00  0.18476-02  0.2289C-0?  -0.36RB6-02  -0.21946-02 

3  -0.74646  00  -0.10836  Cl  0.1640E  01  0.76296  00  0.8394C-03  0. 17896-01  0.526IF-0I  -0.6O45F-02 

4  -0.1315F  00  0.86136  00  0.80826  00  0.52646  00  0.2832E-01  -0.59896  OC  C.245HF  00  0.4771F  00 


FLCHfNT  LOADS  FOR  SLFNOFR  dOOY  NlJMftfft  2.  MODE  NO.  2 

F?  FY  M2  my 


1  0.5902P 

00 

0.62796 

CO 

0.55256 

00 

0.24646 

OC 

0.2  300F-03 

0.15336-03 

-0.3986F-02  -0.26566-0? 

2  0.1109F 

Cl 

C. 21846 

01 

0.6752F 

CO 

0.3240F 

on 

0.38476-02 

0.22396-02 

-0.36836-02  -0.21946-0? 

3  -0. 2601 F 

01 

C. 838 16 

00 

0. 1466C 

01 

0.91686 

DO 

0.83946-03 

0.  1785E-01 

0 . 526 1 r -01  -C. 60456-0? 

4  -D.3217F 

01 

-0.2814F 

01 

0.17126 

no 

0.96996 

00 

0.6 l44t-0 1 

-0.4256F  00 

0.2634F  CO  0.8695F  CD 

LOADS  ON  SLENDER  BODY  FLCMCNTS  IME  TO  LIFTING  SURFACC  BOXES . 

F2  FY  M2  MY 

1  -0.7847C-C3  -C. 11966-02  0.I36OF-OI  0.20716-01  0.4904E-04  C. 747OE-04  -0.R498E-O3  -9.1294F-0? 

2  -0.1440F-01  -0. 2697F-01  0. 39496-01  0.64086-01  0. 8021E-03  0.15366-02  -0.76886-01  -0.1473E-02 

3  -0.10226  00  -0.27  356  00  0.1241C  00  0.32306-01  -0.2875F-C2  0.8202F-02  0.1P42F-C1  0. 2076F-01 

4  0.1360f  00  -0.2694E  00  0.29666-01  0.23406  00  0. 3091E-01  -0.6906F-01  -0.3501F-02  0. 48116-01 


LOADS  UN  SLCNOER  BODY  ELEMENTS  OUF  TO  LITTING  SURFACE  30XcS  AND  1 NTFRFCRENCF  BODY  FL6MENTS . 

F2  FY  M?  MY 

1  -0.8021F-01  -0.95586-01  0.1164F  00  0.1939F-01  0.4904C-04  0.74796-04  -0.94986-03  -0.1294F-0? 

2  -0.9383f -01  -0.1214C  00  0.14236  00  O.6366E-0I  0.3021F-03  0.15366-02  -0.76B8F-03  -C. 14736-02 

3  -0.47616-01  -O.S62Ct  00  0.19096  00  0.1544F  00  -C. 23756-02  0.82026-02  0.18426-01  0.2076E-01 

4  -0.1990F  00  -0.32176-01  O.I218F  DO  -0.M83F-01  0.12106  00  -0.  14786  00  -0.2215F-01  0. 11836  00 


6LCM6NT  LOAOS  FOR  SL6NDFR  6DDY  NUW86R  2.  MflDE  NO.  3 

CZ  FY  M2  MY 

1  -0. I 079F  00  0.56916  00  0.1109F  00  0.1523F  00  0.49046-04  0.74706-04  -0.84986-03  -0.1294F-C2 

2  -0.3431F  00  9.1373C  01  0.92466-01  0.4625E  00  0.80216-03  O.1536E-02  -C.7688F-03  -0.14736-02 

3  -0.1 3776  01  -0.8620F  00  -0.7495F-01  0.15446  00  -0.29756-02  0.8202F-02  0.1842F-01  0.2076C-01 

4  -0.74296  00  -0.26916  01  C.11016-01  -0.5936F  00  0.12L0C  00  -0.1473F  00  -0.22156-01  0.11836  00 

***** 

cxecurioN  tihetminutfs) 

CPU  0.0300 

I/O  0.0319 

TOTAL  0.0619 
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MOO 6  NO  l 


STRIP  NO, 

Y 

7 

YDS 

UFT  COEFFICIENT 

MOMENT  C0ETFIC1FNT  CP- >  A 

nj  CP-  1 

l 

1*1495 

0.5000 

0.3832 

0.4J1427 

4.4J2753 

0.329856 

-0.40lo98  -0.61457 

0. 38574 

2 

1.7165 

0. 5000 

0.5727 

0.4/3788 

4.075587 

C. 495600 

-0.625017  -0.79604 

0. 37019 

3 

h* 2500 

0.62*0 

0.7530 

0,025103 

3.706612 

0.455104 

-0.297408  -17.87962 

0.31026 

4 

2. 7600 

0.8/50 

0.9167 

0.096159 

3.156390 

0.385158 

-0.104931  -J./6641 

0.2M324 

5 

2.000' 

0.7500 

3.6667 

0.MI614 

1.355033 

0.099202 

-0.239099  0.09526 

0.42645 

ROOT  CLfM. 

Y 

l 

m 

Uri  COEFFICIENT 

MOMENT  COEFFICIENT 

1 

0.0 

c.o 

O.OR33 

-0.213106 

1.648189 

0.0 

C.O 

1 

0.0 

0.0 

0. 003  J 

O.C 

0.0 

0.0 

0.0 

2 

0.0 

0.0 

0.2500 

-2. 327499 

4.789781 

0.9 

0.0 

2 

0.0 

0.0 

0.2500 

0.0 

0.0 

0.0 

0,0 

3 

0.0 

0.0 

0.5000 

-4.162803 

3.536185 

0. 605896 

-0.127619 

3 

0.0 

0.0 

0.5000 

0.0 

0.0 

0.0 

O.C 

4 

0.0 

0.0 

0.7500 

-2. « 19400 

-4.243449 

-0.041323 

-0.02870* 

4 

0.0 

0.0 

0.7500 

0.0 

0.0 

0.0 

0.0 

5 

0.0 

0.0 

0.9167 

-0.659796 

-0.907041 

0.0 

0.0 

5 

0.0 

0.0 

0.9167 

0.0 

0.0 

0.0 

0.0 

BODY  ELFM. 

V 

; 

*0l 

MODIFIED  LIFT  COEFFI C IENI Si  7-  ANO 

Y-01RCC1 IONS 

1 

0.0 

0.0 

0.0833 

1.437044 

1.568631 

0.0 

0.0 

2 

0.0 

0.0 

0.25O0 

-2.298344 

6.563907 

0.0 

0.0 

3 

0.0 

0.0 

0.5000 

-8.022783 

12.580199 

0.0 

0.0 

4 

0.0 

0.0 

0.7500 

-6.620901 

-4.239679 

0.0 

).0 

5 

0.0 

c.o 

0.9167 

-1.623099 

-2.1/5225 

0.0 

0.0 

BODY  NH» 

Y 

l 

X-CtMIER 

U1TAI  UFT  COFFF. 

TOTAL  M0MEN1  LOCFF. 

1 

0.0 

0.0 

3.0000 

-2.225857 

1.29J789 

4.735582 

-0.619459 

l 

O.o 

0.0 

3.0000 

0.0 

0.0 

0.0 

0.* 

BODY  FLFM. 

Y 

/ 

XOL 

UFT  COcFF  IC IFNT 

MOMENT  COFFficK  T 

6 

2.0000 

-o.rooo 

0.0500 

0.323774 

1.924397 

0.000217 

0.000469 

6 

2.000') 

-0. 500 J 

0,0500 

0.634755 

0.943069 

-0.003767 

-0.008 . 26 

7 

2.0000 

-0.5  00 

0.1600 

-C,«. 74999 

4.994  20  > 

0.003964 

0.007005 

7 

2.0000 

-0.6400 

0.1*00 

0. 754814 

1.101733 

-0.003799 

>0.00679) 

8 

2. 0*00 

-0.5000 

0.3500 

-l.-  6893 

-0. 7164*2 

-0.903862 

0. "07966 

8 

2.0000 

-0.5000 

0.3500 

0.49966 ) 

C.9J9400 

6.027492 

C. 919173 

9 

2.0000 

-o.soro 

0.7600 

-0.616455 

-1.211120 

0.080586 

-0.109485 

9 

2.0000 

-0.5)00 

0.7539 

-0.  O3-.O60 

0.031491 

-0*012456 

0.102784 

BODY  ELEM. 

Y 

L 

<0L 

MODIFIED  LIFT  CI'EFF IC1FNTS i  7-  AND 

Y-D1RECT IONS 

6 

7.0000 

-0.5000 

0.0500 

0. 740464 

2.146836 

0.705551 

0.56074  < 

7 

2.0000 

-0*5000 

9.16)0 

-0.314993 

3.342095 

C. 703525 

0.8331 J6 

8 

2.0000 

-0.5000 

0.3600 

-2.589956 

-0.5781*3 

0.96*40? 

1.410345 

9 

2.0000 

-0.5000 

0.7609 

-*.869009 

-3.230846 

-C. 112961 

0.158218 

BODY  NO. 

Y 

l 

X-CENTFK 

F3TAL  UFT  COEFF. 

TOTAL  MOMfNT  COfFF. 

2 

2.0000 

-0.5000 

3.2500 

-0.  204524 

-0.050244 

0.259284 

0.276117 

2 

2.0000 

-0.5000 

3.2500 

0. 105513 

0.184000 

-0.03399C 

-0.052717 

<*l  ■ 

-2.492867 

S* 429*83 

CY  «  0.0 

0.0 

CM  ■ 

2.574017 

-1.015308 

CM  »  0.0 

0.0 

CSf 

0.0 

0.0 
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MODE  NO  2 


STRIP  NJ. 

Y 

l 

YOS 

LIFT  COEFFICIENT 

MOMENT  COEFFICIENT 

CP-R  i 

VNO  CP-1 

1 

1.1405 

0.5000 

0.3832 

5.984121 

9.267561 

-0.163195 

-2.088826 

0.27727 

0.47539 

2 

1.7165 

0.5000 

0.5722 

6.435728 

10.352689 

C. 200348 

-2.326933 

0.21887 

0.47477 

3 

2.25CO 

0.6250 

0. 7500 

4.043768 

8.536967 

0.592437 

-1.468541 

0.19349 

0.42202 

4 

2.7500 

0.8750 

0.9167 

3.696903 

7.31 3690 

0.686186 

-0.897620 

C. 96377 

0.37273 

5 

2.0000 

0.2500 

0.6667 

’.892935 

2.142693 

-0.143624 

-0. 718103 

0.29966 

0.58523 

onnY  hem* 

Y 

l 

xut 

LIFT  COEFFICIENT 

MOMENT  COEFFICIENT 

i 

0.0 

0.0 

0.0833 

1.628965 

1 .1043FJ 

0.0 

0.0 

i 

0.0 

0.0 

0«O8’3 

.0 

0.0 

0.0 

9.0 

2 

0.0 

0.0 

0.2500 

e* 501630 

9.481941 

0.0 

0.0 

2 

0.0 

0.0 

0.2500 

0.9 

0.9 

0.0 

0.0 

3 

0.0 

0.0 

0.50'JO 

-4.35292  i 

16.024943 

0.887204 

-1.607999 

3 

0.0 

0.0 

0.5000 

O.C 

0.0 

0.0 

0.0 

4 

0.0 

0.0 

0.7500 

-12.276346 

-8.129966 

-0.IH233 

-0.023787 

4 

o.c 

0.0 

0.7509 

0.0 

0.0 

0.0 

0.0 

6 

o.c 

0.0 

0.9167 

-2.541182 

-3.845133 

0.0 

0.0 

4 

0.0 

0.0 

9. 9U7 

0.0 

C.O 

c.o 

0.9 

BODY  ELFM. 

Y 

l 

X1L 

MooiFieo 

LIFT  COEFFICIENTS.  Z-  AND 

y-direct ions 

1 

0.0 

0.0 

0.0833 

3.075882 

-1.427946 

0.0 

0.9 

2 

0.* 

9.0 

0.2500 

3.6601 36 

10.856790 

0.0 

0.0 

3 

>.o 

0.0 

0.5000 

0.109172 

46.399570 

0.0 

0.0 

4 

n  } 

0.0 

0.7509 

-22. 057449 

-2.866698 

0.0 

0.0 

5 

0.0 

0.0 

0.9167 

-7.123833 

-6.340250 

0.0 

o.o 

BODY  NO. 

Y 

7 

X-CENTER 

TOTAL  LIFT 

COEFF. 

TOTAl  MOMFNT  COFFF. 

1 

0.0 

0.0 

3.0000 

-3.023586 

4.790936 

9.60921C 

-5.578701 

1 

0.0 

c.o 

3.0000 

0.0 

0.0 

0.0 

0.0 

BODY  CLEM. 

Y 

L 

xnt 

LIFT  COEFFICIENT 

MOMENT  COEFFICIENT 

6 

2.0000 

-0.5000 

0.0500 

2. 360761 

2.51  15  39 

0.000920 

0.000613 

6 

2.0000 

-0.5000 

0.0500 

2.20980  9 

C. 985413 

-0.015943 

-0.010623 

7 

2.0000 

-0.5000 

0.1590 

4.437428 

8.736561 

0.015389 

0.009157 

7 

2.0000 

-0.5000 

0. 1500 

2.700893 

1.295815 

-0.014750 

-0.008777 

8 

2.00CC 

0.50CD 

0. 3500 

-3.468510 

1.117612 

0.001119 

0.023007 

8 

2.0CCC 

-0.5000 

0. 3500 

1.95464  f 

1.249996 

0.070149 

-0.008061 

9 

2.0000 

-0. 5000 

0.7500 

-2.5/3796 

-2.250  )78 

0.949156 

-3. 340607 

9 

2.0000 

-0.5000 

0.7500 

3.  136960 

C. 455902 

0.210737 

0.456594 

BODY  HEM. 

Y 

Z 

XflL 

MODIFIED 

LIFT  coefficients,  z-  and 

y-dikcctions 

6 

2.0000 

-0.5000 

0.0500 

3.423564 

2.662498 

1.930327 

0. 205229 

7 

2.0000 

-0.5000 

0.1500 

2.955782 

6.027204 

2.287436 

0.71861 1 

8 

2.0090 

-0.5000 

0. 3500 

-5.456857 

3.  166936 

3.402420 

2.092484 

9 

2. 00r0 

-0.5000 

0. 7500 

-7. 399488 

-5.741970 

0.469216 

1.369649 

BODY  NO. 

Y 

l 

X-CCNTER 

TOTAL  LIFT 

COEFF. 

TOTAL  MOMENT  COEFF. 

? 

2.00^0 

-0.5C0C 

3.2600 

-0.643606 

0.130693 

0.821367 

0.337447 

2 

2.0000 

-0.5000 

3.2500 

0.  44  7f  5 

0.32453 2 

-0. 1 68504 

-0.154367 

C7  • 

1.013682 

14.231640  ri 

*  0. 

0 

0.0 

CM  * 

4.635420 

-5.763388  CN 

-  9. 

,0 

0.0 

CSI  • 

0.0 

0.0 
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MOOF  NO  3 


SHI*  N*. 

Y 

2 

YOS 

LIFT  COEFFICIENT 

MONFNT 

COEFF ICIFNT  CM-*  ANO 

i  CP-I 

1 

1.14** 

O.5CO0 

0. 3*1> 

0.40824 ) 

3.398235 

0.290571 

-0. 590ol8  -0.46176 

0.42616 

2 

I,  fU* 

0.50*0 

0.5722 

0.03)454 

4.458212 

0.532)58 

-0.629056  -15.66328 

0.39089 

1 

*.***• 

0.42*0 

3. 25*0 

-0.603054 

4.978993 

0.604063 

-0.295527  1.23704 

0.30935 

4 

2. 1500 

0.3750 

3.0147 

-0.53046) 

5.105230 

0.66)15) 

0.033554  1.00319 

0.24343 

5 

2*0990 

0.2*00 

0.666? 

0.455)74 

0.420  309 

0.00881 3 

-0.168)70  0.2)065 

0.65059 

8"0y  IUm. 

Y 

l 

<01 

LIFT  COEFFICIENT 

MOMENT 

COEFFICIENT 

1 

0.0 

0.0 

3.33)3 

0.042304 

0.099540 

0.0 

0.0 

l 

0.0 

0.0 

0.04)3 

0.0 

0.0 

0.0 

0.0 

2 

0.0 

0.253* 

0.0428*3 

0.099540 

0.0 

0.0 

7 

0.0 

3.C 

0. >5)0 

0.0 

0.0 

0.0 

O.C 

\ 

0.? 

0.0 

0.6000 

-0. 19565  * 

3.223248 

0.562)77 

-0.060)6) 

) 

0.0 

0.0 

0.5000 

0.0 

0.0 

0.0 

0.0 

4 

0.0 

0.0 

3. 7500 

-0.678054 

0.073760 

-0.049384 

-0.087266 

4 

0.0 

0.0 

0. 7400 

O.C 

0.0 

0.0 

o.o 

9.0 

0.0 

0.0167 

-0.440516 

0.427825 

0.0 

0.0 

* 

3.0 

0.0 

0.4167 

0.0 

0.0 

0.0 

0.0 

?oot  eiM. 

Y 

l 

xot 

MODI  FIFO  LIFT  C0EFFI< 

Cl  ENTS.  Z-  ANO  Y-01FECT IONS 

i 

0.0 

0.0 

0.08)3 

0.361440 

0.043729 

0.0 

0.0 

2 

0.0 

0.0 

0.2500 

-0.6)0596 

-0.051170 

0.0 

0.0 

1 

o.c 

3.0 

0.5000 

1.2154  »i 

9. >*40580 

0.0 

0.0 

* 

0,9 

0.0 

0.7500 

-1.268764 

0.836*  63 

0.0 

0.0 

* 

0.0 

0.0 

0.0167 

-0.312666 

0.9061 7  3 

0.0 

0.0 

800 Y  NO. 

Y 

l 

X-CENTER 

TOTAL  LIFT 

COfFF. 

TOTAL  MO**FNT  COEFF. 

1 

0.0 

3.0 

3.0000 

-0.291291 

1.119119 

0.94)517 

-2.3422*6 

„ 

0.0 

0.0 

3.0000 

0.0 

0.0 

0.0 

0.0 

e*w  tifM. 

1 

l 

xni 

LIFT  COEFFICIENT 

U0MF  NT 

COEFFICIENT 

& 

2.000* 

-o.ooco 

0.0500 

-0.431625 

2.27*594 

0.000196 

0.000299 

5 

2.0000 

-0.5000 

0.0500 

0.44)492 

0.609)56 

-0.003399 

-0.005177 

1 

2.30OC 

-0.5000 

0.1500 

-1.372)94 

7.491276 

0.00)208 

0.006145 

1 

- 0 . 5SC0 

0. 1500 

1.84007) 

-0.00)075 

-r .004806 

8 

2.0000 

-0.5000 

0.3500 

-1.836075 

-0.749)41 

-0.0038)3 

0.010936 

8 

2.0500 

-0.5000 

1.  3500 

-0. 099934 

0.205834 

0.024563 

0.027676 

9 

2.000* 

-0.5000 

0.7500 

-0.594350 

-2 . 152955 

0.096763 

-0.U8245 

0 

2.0500 

-0.5000 

0.7500 

0.008805 

-0.474388 

-0.017719 

0.094607 

800Y  HE«. 

Y 

l 

xou 

MUOIFIEO 

nrl  COEFFICIENT*,  i-  AND  V-DIHEC f  1-INS 

6 

2.0000 

-0.5000 

0.0500 

0.451001 

2.983956 

0.572218 

0.627556 

7 

2.C0OO 

-0.4000 

0.1500 

-1.039745 

4.824802 

0.344584 

1.192759 

S 

2.o«ro 

-0.5000 

0.3500 

—  3 .49091  ) 

-0.419818 

-0.217383 

0.532570 

0 

2.3000 

-0.5000 

0.7500 

-1.962565 

-5. 7782C6 

-0.053222 

-1.257548 

800 Y  Nil  • 

V 

2 

X-CENTE* 

TOTAL  LIFT 

COFFF. 

TOTAL  MOMENT  COEFF. 

2 

2.0000 

-0.5000 

3.2500 

-0.401718 

-0.126736 

0.297814 

0.48188? 

2 

2.0000 

-0.5000 

3.2500 

0.021779 

r. 027438 

-*.0009)0 

0 .00601 1 

C2  • 

- 1 . 2533)1 

5.24512) 

CY 

•  0 

.0 

0.0 

CM  - 

1.816684 

-2.3825 

76 

CN 

•  0 

.0 

0.0 

CSl  * 

0.0 

0.0 
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M30E  NO  l 


PANFL  NO 

STRIP  NO 

BOX  NO 

xoc 

X 

Y 

7 

PRESSURES 

l 

1 

1 

0.12500 

2.35938 

1.14952 

0.50000 

1  .966562 

4.242338 

1 

1 

2 

0.62500 

3.29688 

1.14952 

0.60000 

-1.103708 

4.623160 

1 

; 

3 

0.12500 

2.57813 

1.71651 

0.50000 

2.693086 

4.813315 

1 

2 

4 

0.62500 

3.39063 

1.71651 

0.50000 

-1.745509 

4.937862 

2 

3 

5 

0.12500 

2.79688 

2.25000 

0.625C0 

1.858071 

4.368785 

2 

3 

6 

0.62500 

3.48438 

2.25000 

0.62500 

-1.807866 

3.042441 

2 

4 

7 

0.12500 

3.01563 

2.75000 

0.87500 

1  .684071 

4.314062 

2 

4 

a 

0.62500 

3.57813 

2.75000 

0.87500 

-1.492552 

1.997920 

3 

$ 

9 

0.12500 

2.68750 

2,00000 

0.25000 

1.359550 

1.076156 

3 

5 

10 

0.62500 

3.43750 

2.00000 

0.25000 

-0.076321 

1.633911 

THE 

2)  HJOIPIEO 

H-EIEMENTS  FOP 

MODE 

NO.  1 

1.59464 

1.59464 

1 • 3*202 

1. 38202 

C. 99851 

0.99851 

0.81696 

0.81696 

0.0 

0.0 

1.50000 

1.50000 

1.50000 

1.50000 

1.50000 

3.00000 

3.00000 

3.000C0 

3.90C00 

0.0 

0.0 

0.0 

0.0 

THE 

9  MODIFIED 

0H/DX-6LFMFNTS 

FOR 

MODE  NO.  1 

0.0 

0.0 

0.0 

0.0 

0.0 

o.c 

O.o 

0.0 

0.0 

OEFICCTIPN  moof 

PRESSURE  M'JOC 

PANELS 

O.ENERAL  I  260 
ONLY 

FORCES 

TOTAL 

1 

t 

0. JD132967E 

00 

0.41 844883F  01 

-0.249356946 

01 

0.53777828E 

01 

1 

2 

0.527738006 

01 

0. 90802803E  01 

0 .966585 16E 

00 

0.  141325986 

02 

l 

3 

-0.14895672E 

00 

0.431721 78F  01 

-0. 1 2436838F 

01 

0.51828642F 

01 

PANEL  NO 

STRIP  NO 

BOX  NO 

XOC 

NODE  NO  2 

X 

Y 

l 

PRESSURES 

t 

1 

1 

0.12500 

2.35938 

1.14952 

0.50000 

8.323404 

5.54604$ 

1 

1 

2 

0.62500 

3.29688 

1.14952 

O.SUOOO 

3.644041 

12.909089 

1 

2 

3 

C.  12590 

2.57813 

1.71651 

0.50000 

10.454984 

6.221313 

1 

2 

4 

C. 62500 

3.39963 

1.71651 

0.50000 

2.416474 

14.484083 

2 

3 

5 

0.12500 

2.79638 

2.25000 

0.62500 

8.435404 

6.931283 

2 

3 

6 

0.62500 

3.48438 

2.25000 

0.62500 

-0.347866 

10.142653 

2 

4 

7 

0.12500 

3.01563 

2. 75000 

0.87500 

7 . 990 103 

7.380054 

2 

4 

8 

0.62500 

3.57313 

2. 75000 

0.87500 

-0.996293 

7.247329 

3 

5 

9 

0.12500 

2.68750 

2.00000 

0.25OC0 

3.764908 

0. 340826 

J 

5 

10 

0.62500 

3.41750 

2.00000 

0.2500C 

2.02C964 

3.944560 

41 


2 


2.  *0824 
1.6424) 

3. C30O0 
7.78000 


I.0009O 
2  *00000 


lHf 

23  Mooineo 

H-EL6MLNTS  FOR 

Moot 

NO.  2 

3  *  4Q489 

2.  37538 

3.12395 

1.86180 

2.31945 

1.94*79 

0.0 

0.0 

0.50000 

1.50000 

4.50030 

5.50000 

4.250C0 

4.750^0 

5.  75000 

o.c 

0.0 

0.0 

c.o 

IMF 

9  MODIFIED 

OH/OX-FL6M6NTS 

FOR 

MODE  no.  2 

l.WVOO 

1.00000 

l.OCOOO 

l.COOOO 

2.00000 

2.01000 

2.00000 

OFFLFCHON  «00f  PrfFSSURr  MUPE 


GFNCRAL ! 7bD  FORCFS 
PANfL$  ONLY  TOTAL 


0.1 1 8030  7  IE  00 
0.95952950E  01 
■0.85063422E  00 


0.83292486F  Ot 
0.189728196  02 
0. 8 7O69708F  01 


-0.860979276  01 
-0.2q987698F  01 
-0. "31  874332E  01 


0 . 8040292 7f  01 
0.23680984F  02 
0.964641676  01 


MOOF  V)  3 


PANFL  NO 

STRIP  NO 

OCX  NO 

xoc 

X 

Y 

l 

PRESSURES 

1 

l 

1 

0.12500 

2.36938 

1.14952 

0.80000 

1.774650 

2. 702881 

1 

l 

2 

0.62500 

3.29688 

1.14952 

0.80000 

-0.958163 

4.093591 

1 

> 

3 

0.12500 

2.57813 

1.71651 

0.50000 

2.179616 

4.176091 

l 

2 

4 

0.62500 

3.39063 

1.71651 

0.50C00 

-2.112709 

4.741335 

2 

3 

5 

0.125C0 

2.79008 

2.25000 

0.62500 

1.698914 

6.286385 

2 

3 

6 

0.62500 

3.48438 

2.25000 

0.62500 

-3.086831 

3.671605 

2 

4 

7 

0.12500 

3.01563 

2.75000 

0.87800 

1.33l9ie 

7.792065 

2 

4 

ft 

0.62500 

3.57813 

2.75000 

0.87500 

-3.092846 

2.418398 

3 

9 

0.12500 

2.68750 

2.00000 

0.25000 

0.718314 

-0.043016 

3 

5 

10 

3.62500 

3.43750 

2.00009 

0.25000 

0.192435 

0.883634 

THC 

23  MUOIFIFO 

H-FL6MENT  S  FOR 

MODE 

NO.  3 

1.372  74 

1.37274 

1.97330 

1.97330 

1.84577 

1.84577 

2.03752 

2.03752 

-0.09175 

-0.09375 

C.O 

0.0 

0.0 

C.O 

0.0 

4,99999 

9.99999 

4.99999 

4,90999 

1.00900 

i.oooco 

1.00000 

1.00000 

THF 

13  MODIFIED 

om/dx-ele^ents 

FOR 

MOOF  NO.  3 

0.0 

O.C 

0.0 

O.C 

C.O 

0.0 

0.0 

C.O 

0.0 

0.0 

0.0 

0.0 

0.9 


06FLFCM'’ 

\ 

PRESSURE  MOOE 

PANELS 

generalized  i 

ONLY 

FORCES 

TOTAL 

3 

1 

0.356097886 

00 

0.60I03865E  01 

-0.521972786 

00 

0.5965S657E 

01 

3 

2 

0.733997346 

01 

0.  132518856  C2 

0.549304776 

01 

0.139038016 

0? 

> 

3 

-0.51878241C 

oc 

0.687810716  01 

-0. 1 843522  If 

01 

0.647394756 

01 

EXECUTION  HvCIHINUnS) 
CPU  '.0410 

I/O  9.037* 

TOTAL  0.<"  788 
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3.0  BLANK  AND  LABELED  COMMON  BLOCKS 
3.1  Blank  Common  Block 


ITEM 

NO. 

SYMBOL 

MNEMONICS 

DESCRIPTION 

1 

NP 

No.  of  panels 

Input 

9 

2 

NTP 

Total  no.  of  boxes  on  all  panels 

3 

NB 

No.  of  bodies 

Input 

10 

4 

NBZ 

No.  of  bodies  with  z-orientation 

5 

NBY 

No.  of  bodies  with  y-orientation 

6 

NTZ 

Total  no.  of  body  elements  with 

z-orientation 

7 

NTY 

Total  no.  of  body  elements  with 
y-orientation 

8 

FMACH 

Mach  Number 

Input 

Items 

2 

9 

A 

REFA 

Reference  area 

3 

10 

S 

REFS 

Reference  span 

4 

11 

c 

REFC 

Reference  chord 

5 

12 

6 

ND 

Symmetry  flag  for  y  -  0  plane 

7 

13 

e 

NE 

Symmetry  flag  for  z  =  0  plane 

8 

1.  No.  of  interference  elements 

(NBE)  for  all  bodies 

14 

NBEA( 10,2 

2.  Z-Y  flag  (NZY)  for  all  bodies 

15 

NSBEA(IO) 

No.  of  slender  body  elements 
(NSBE)  for  all  bodies 
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fill 

SYMBOL 

MNEMONICS 

DESCRIPTION 

15 

NSBEA(IO) 

No.  of  slender  body  elements 
(NSBE)  for  all  bodies 

16 

nc 

NCARAY(IOO) 

No.  of  chordwise  boxes  for 
all  panels 

17 

ns 

NSARAY(IOO) 

No.  of  s panwise  strips  for 
all  panels 

18 

nba 

NBARAY(IOO) 

nba.  =  e  nc-ns.,  where 
i  jsl  J  J 

i  =  1,  NP 

19 

zc 

ZB(10) 

z-coordinate  of  body  center 

• 

20 

yc 

YB(10) 

y-coordinate  of  body  center 

21 

xl.E. 

XLE(IO) 

x-coordinate  of  leading 
edge  of  body 

Eq. 

(5.2.1-17) 

22 

XT.E. 

XTE(IO) 

x-coordinate  of  trailing 
edge  of  body 

Eq. 

(5.2.1-18) 

23 

m 

ARB(IO) 

Cross  sectional  aspect 
ratio  of  body 

24 

a 

AVR(IO) 

Average  characteristic 
half-width  of  body 

Eq. 

25 

a0 

AO (2 00) 

Local  body  half-width 

(y-di recti  on;  radius 

for  circle)  Slender 

(5.2.1-22) 

26 

a0 

A0P(200) 

x-derivative  of  body 
half-width  arra*s 

Eq. 

(5.2.1-23) 

27 

csi 

XIS1 (200) 

x-coordinate  of  slender 
body  element  leading  edge 

Eq. 

(5.2.1-20) 

28 

5S2 

XIS2(200) 

x-coordinate  of  slender 

body  element  trailing  edge 

Eq. 

(5.2.1-21) 
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SYMBOL 

MNEMONICS 

XI  J(200) 

c 

CS  (2.10) 

yPs.nPs 

YS{210) 

zPs,£;Ps 

ZS(2;10) 

es 

EE(200) 

sin  Y$ 

SG(200) 

COS  Ys 

CG(200) 

xP»  xl 

X(500) 

CP 

XIC(500) 

k 

KR 

r 

XM 

XM 

AxP,AxI 

DELX(500) 

tan  a 

XLAM(500) 

RIA(IOO) 

DESCRIPTION 


x-coordinate  of  leading  edge 
of  strips  (on  all  panels) 

chord! ength  of  strips 

y-coordinate  of  strip  center- 
lines  and  body  center-lines 

z-coordinate  of  strip  center- 
lines  and  body  center-lines 

Half-width  of  strips 


of  the  dihedral 
cosine  angle  of  strip 

3/4  chord  x-coordi nates  of  all 
boxes  (receiving  points)  and 
interference  body  section 
midpoints 

1/4  chord  x-coordi nates  of  all 
boxes  (sending  points) 

Reduced  frequency 
Moment  axis 

Average  chord-lengths  of  all 
boxes,  and  interference 
body  element  lengths 

Tangent  of  the  sweep  angle 
of  the  1/4-chord  line  (bound 
vortex)  of  all  boxes 

Radii  of  all  body  inter¬ 
ference  elements 


Eq. 

(5. 2. 1-9) 

Eq. 

(5.2.1-10) 

Eq. 

(5.2.1-11) 

Eq. 

(5.2.1-14) 

Eq. 

(5.2.1-19 

Eq. 

(5.2. 1-1 
and  -2) 


Input  Item  6 
Eq. 

(5.2. 1-7 
and  -8) 


ITEM 

NO. 

SYMBOL 

MNEMONICS 

DESCRIPTION 

43 

44 

THIA(IOO) 

TH2A( 1 00 ) 

ei's  for  all  bodies 

02 's  for  all  bodies 

45 

NFL(IO) 

Number  of  body  sections  with 

ei -distribution  -  for  all  bodies 
(max  is  3  per  body) 

46 

I FLA( 30 ,2 ) 

1.  Sequence  number  of  the  first 
body  element  with  ei  distribution 

2.  Sequence  number  of  the  last 

body  element  with  ei -distribution  - 

for  all  bodies 

47 

NTT  2(10,2) 

1 .  Number  of  ei ’s 

2.  Number  of  02  's 

all 

bodies 

48 

NAS ( 1 00 ) 

Number  of  associated  bodies  for 

all  panels 

49 

NASB(200) 

Associated  bodies  - 

all  panels 

50 

YIN ( 1 00 ) 

y-coordinate  of  inboard  edge  of  panel 

51 

ZIN(IOO) 

z-coordinate  of  inboard  edge  of  panel  - 
all  panels 

52 

An 

DETA( 50) 

y-shift  elements 

all  panels 

53 

54 

55 

AC 

DZET(50) 

NOBODY (50) 

NPANEL(50) 

z-shift  elements  and  a11  bodles 

Body  number 

for  the  An,  Ac  pairs 

Panel  number 

56 

NTO 

NTP  +  BB  NBE. ,  where 
i=l  1 

NTP  =  total  number  of  boxes  on  all 
panels  and  NBE^  =  number  of  interference 
elements  on  body  i 

46 


Subroutines  using  the  Blank  Common  Block  are: 

MAIN,  DATA,  GEND,  SUBP,  SUBB,  SB,  BFM,  AERO  and  GENF. 


3.2  Labeled  Common  Blocks 

Common  DLM 


ITEM 

NO. 

SYMBOL 

MNEMONICS 

DESCRIPTION 

1 

K,<s> 

k2<s> 

KIO 

Planar  part  of 

the  steady 

contribution 

2 

K20 

Nonplanar  part  of 

to  the  kernel 

3 

ReK-| 

K1RT1 

Re(  aK-j  )T-j 

4 

ImK1 

K1IT1 

Unsteady  part 

of  modified 

5 

ReK2 

K2RT2P 

Re(AK2)T2* 

kernel 

6 

ImK2 

K2IT2P 

Im(AK2)(T2* 

7 

K10T1 

K1  T1 

8 

K20T2P 

K  (S)  T  * 

1^2  l2 

Subroutines  using  Common  DLM  are: 
FLLD,  INCRO  and  TKER. 
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Common  KDS 


ITEM 

NO. 

SYMBOL 

MNEMONICS 

DESCRIPTION 

1 

IND 

Flag,  0  or  1.  -  If  0,  the  total 
kernel  is  computed;  if  1,  the 
incremental  part  only. 

1 

Re(^) 

KD1R 

'X* 

Real  part  of 

% 

El 

Im(K1) 

KD1I 

Imaginary  part  of  K-j 

'V 

Re(K2) 

'X, 

KD2R 

Real  part  of  K2 

i 

Im(K2) 

KD2I 

Imaginary  part  of  K2  - 
and  K2  are  defined  in 

Sec.  5. 5. 1.4  Subroutine  FLLD. 

Subroutines  using  Common  KDS  are: 
FLLD,  INCRO  AND  TKER. 
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4.0  LOGICAL  TAPE  UNITS 


4.1  Tape  Numbers  and  Symbols 


TAPE  NUMBER 

SYMBOLS 

PRESENT 

ASSIGNMENT 

USER  SUBROUTINES 

NTAPE(l), 

1 

GEND,  SOLVIT 

ITP1 ,  NM, 
NDZ 

MAIN 

NTAPE(2), 

2 

GEND,  SOLVIT 

ITP2,  NO, 
MTAPE ,  NDY 

MAIN 

NTAPE(3) , 

3 

AERO,  BFM,  GEND,  GENF, 

ITP3,  NW, 
NPTAP,  NDW 

SB,  SOLVIT,  WANDWT,  MAIN,  RWREC 

NTAPE(4) , 

4 

MAIN,  GEND,  SB,  WANDWT 

ITP4,NW5,NWT 

NTAPE(8) , 

8 

MAIN,  GEND,  AERO,  GENF 

ITAPE, 

NEWBFM 

NTAPE(9), 

9 

BFM,  SB,  MAIN 

NPSTAP.NW2 

NTAPE(IO) , 

10 

MAIN,  SB,  SOLVIT 

NI,  NW3 

NTAPE(ll), 

11 

AERO,  BFM,  SB,  MAIN,  BFSMAT 

NBFM,NW4 

In  addition  to  the  above  tapes,  units  #5  and  6  are  used  throughout  the 
program  as  the  standard  ‘card-read1  and  'print'  units.  The  following  table 
explains  the  use  of  the  logical  tape  units  within  the  user  subroutines  in 
the  order  of  the  program  flow. 
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4.2  Table  of  Input  Output  Units 


SUBROUTINE 

NAME 


WANDWT 


SOLVIT 


BFSMAT 


MAIN 


TAPE  NUMBER 

INPUT 

CONTENTS 

OUTPUT 

CONTENTS 

SCRATCH 

I  TAPE 
NWT 


NPSTAP 

NN 


[DT  !  WT] 

i 


I  TAPE 


NW1=NDW 

NW2 

=  NPSTAP 


{WT} 

={W}+{AW} 


Solutions 
{P}  NO 


Bps** 


NEWBFM 


NBFM 


NEWBFM  BFM* 


*  BFM  represents  all  the  body  forces  and  moments  for  all  modes. 

**  BFS  represents  the  body  forces  FZ  and  FY  described  in  Sec.  5.8.1. 
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5.1  Segment  1 

5.1.1  MAIN 


5.0  DESCRIPTION  OF  SUBROUTINES 


Functional  Description 

The  MAIN  part  of  program  N5KA  reads  the  header  card,  the  reference 
variables,  control  and  print  flags,  and  the  array  of  reduced  frequencies  for 
the  case.  It  calls  subroutine  DATA  which  computes  the  basic  data  arrays  - 
these  are  saved  in  the  blank  common  block.  The  other  10  major  subroutines  are 
called  from  MAIN  in  an  overall  frequency  do-loop  for  each  element  of  the 
frequency  array. 

Input  Output  Variables 


MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

FREQ(IO) 

kr 

IN 

Card 

Reduced  frequency  array 

HDR( 15) 

IN 

Card 

Header  information,  one  card 

IA(2, 

150,3) 

IN 

AR6  of 
RDMODE 

See  Section  5.4.1 

NA(3) 

ARG  of 
RDMODE 

DT ( 500) 

[DT] 

IN 

I/O 

One  row  of  the  [DT]  matrix  gener- 

(I TAPE) 

ated  in  subroutine  GEND 

RHS(IOO) 

OUT 

One  row  of  right-hand-sides  for  one 
unknown  (panel  box,  or  interference 
body  element)  and  all  modes 

WORK 

(10000) 

[WT] 

IN 

I/O 

(NWT) 

Complete  mode  matrix  for  all  un¬ 
knowns  and  modes  generated  in  sub- 

routine  WANDWT 

RA  (2, 
10000) 

OUT 

A  two-dimensional  'real'  work  array 
equivalenced  with  the  complex  work 
array  WORK  (10000)  used  in 
subroutine  SOLVIT 

NTAPE (20) 

OUT 

Data 

initia¬ 

lization 

Variable  name  of  the  logical  tape 
unit  array 

I ERROR 

OUT 

MAIN 

Error  flag  initialized  to  0 

NWORK 

OUT 

MAIN 

Work  array  size 
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IN/OUT 


SOURCE 


DESCRIPTION 


IN  NTAPE 
(20) 


IN  Card 


See  the  argument  list  description  of 
the  user  subroutines 


Mach  Number 
Reference  area 
Reference  semi  span 
Reference  chord 
Moment  axis 
Symmetry  flag 
Second  sys.netry  flag 
No.  of  panels 
No.  of  bodies 

No.  of  reduced  frequencies 


Print  and  control  flags 

Option  flag  for  body  force  calculation 
method;  0  or  1 

Number  of  logical  tape  unit  containing 
body  forces 
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IN/OUT 

SOURCE 

DESCRIPTION 

EM 

6 

6  =\/T  -  M2 

■n 

Print  flag  for  S0LVIT 

IBB 

Control  flag  for  GENF 

NPRINT ,N3 

Print  flag  for  GEND 

NT0T 

Total  number  of  unknowns 

NTIBE 

Number  of  interference  body  elements 

NTSBE 

Number  of  slender  body  elements 

N 

OUT 

MAIN 

Number  of  panel  boxes  plus  number  of 
interference  body  elements 

NSTRIP 

Total  number  of  strips  on  all  panels 

See  NTOT  above 

Total  dimension  of  the  real  work-array  RA 

See  NTOT  above 

KR 

S 

j-th  element  of  reduced  frequency  array 

MD,NMT 

IN 

Arg  list 
of 

RDMODE 

Total  number  of  modes 

NM0DE,MD 

NIM 

OUT 

Total  number  of  mor'es 

NFC 

See  Section  5. 5. 8.1 

NWKR 

JX1 

Index  of  first-(JXl)  and  last  (0X2) 

JX2 

IN 

MAIN 

elements  in  the  modal  matrix  array  when 
read  from  logical  unit  number  NWT  one 

column  at  a  time 

JM 

Indices  used  in  the  do  loop  transposing 

JJ 

the  modal  matrix  array 

JX 

Subroutines  Called  DATA,  GEND,  RDMODE,  SB,  WANPWT,  RWREC,  SOLVIT,  BFM,  BFSMAT, 
AERO  and  GENF 


Common  Blocks  Blank  Common  Block 
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Begin  DO  170 


Read  one  row 
of  DT  matrix; 
,from  tape 
ITAPE 


y 

_ 

Augment  DT 
matrix  rows  by 
the  right  hand 
sides  WT 

_ 3 

r 

Write  augmented; 
matrix  rows 
k(DTjWT)  on 
tape  NI 


Ifi 


6.1.2  Subroutine  ATAN3  (Y,  X,  T) 

Functional  Description 

Subroutine  ATAN3  evaluates  T  =  atan  (Y,X)  by  considering  the  signs  of 
both  Y  and  X,  and  thereby  providing  a  result,  T,  that  lies  in  the  proper 
quadrant.  The  resulting  angle,  T,  is  returned  to  the  calling  program  in  radians 

Calling  Subroutine  DATA 

5.1.3  SUBROUTINE  DZY  (X,  Y,  Z,  SGR,  CGR,  XII,  XI2,  ETA,  ZETA,  AR,  AO,  KR, 

CBAR,  BETA,  FMACH ,  IDZDY,  DZDYR,  DZDYI) 

Functional  Description 

Subroutine  DYZ  calculates  the  effect  of  a  finite  length  doublet,  a  multi 
finite  length  doublet,  or  a  trapezoidal  vortex,  of  unit  strength,  at  a  field 
point.  Each  slender  body  is  composed  of  a  series  of  segments  or  elements. 

This  subroutine  calculates  the  effect  of  such  elements  on  field  points  (lift¬ 
ing  surface  boxes,  etc.).  If  the  cross  sectional  shape  is  circular  then  a 
doublet  is  used.  If  AR  >  1,  two  doublets  are  used.  If  AR  <  1  a  trapezoidal 
vortex  is  used.  This  subroutine  is  used  for  bodies  oriented  in  both  the  z- 
and  y-di recti cns. 

Input  Output  Variables 


MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

X 

X 

x-coordinate  of  the  receiving  point 

Y 

Y 

y-coordinate  of  the  receiving  point 

Z 

Z 

z-coordinate  of  the  receiving  point 

AO 

ao 

Radius  of  sending  body 

AR 

/R 

Aspect  ratio  of  sending  body 

KR 

kr 

IN 

ARG 

Reduced  frequency 

CGR 

COSYr 

Cosine  of  receiving  point  dihedral  angle 

ETA 

n 

y-coordinate  of  the  sending  strip 

SGR 

sinyr 

Sine  of  receiving  point  dihedral  angle 

XII 

U 

Leading  edge  of  sending  element 

XI2 

K2 

Trailing  edge  of  sending  element 
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MNEMONIC  SYMBOL 


IN/OUT  SOURCE 


DESCRIPTION 


BETA 

8 

/l  -  M  z 

CBAR 

c 

IN 

ARG 

Reference  chord  length 

ZETA 

c 

z-coordinate  of  the  sending  strip 

DZDYI 

HI 

OUT 

ARG 

Imaginary  part  of  Dz  or 

DZDYR 

H 

OUT 

ARG 

Real  part  of  D  or  D 
z  y 

M 

MACH  number 

ffifgEB 

IN 

ARG 

FLag  indicating  whether  Dz  or  Dy  is  to 
be  calculated 

Calling  Subroutines  ROWDYZ 

Called  Subroutines  FLLD,  TVOR 


Flow  Chart  -  Subroutine  DZY 
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5.1.4  SUBROUTINE  FLLD  (XOl  ,  X02,  YO,ZO,  SGR,  CGR,  SGS,  CGS,  KR,  CBAR, 
FMACH,  E,  L,  KD1R,  KD1I,  KD2R,  KD2I) 

Functional  Description 

This  subroutine  calculates  the  velocity  normal  to  a  surface  of  dihedral, 
due  to  a  finite  length  l_ine  doublet. 


Input  Output  Variables 


MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

E 

e 

Semi  width  of  singularity 

L 

Option  flag  for  use  in  subroutine  TKER 

KR 

kr 

Reduced  frequency 

YO 

y-n 

Difference  in  lateral  coordinates  of 

receiving  and  sending  points  respectively 

ZO 

2-? 

Difference  in  vertical  coordinates  of 
receiving  and  sending  points  respectively 

CGR 

cosYr 

CGS 

C0SYS 

IN 

ARG 

SGR 

sinyr 

SGS 

sinYs 

XOl 

X-?1 

Difference  in  longitudinal  coordinates 
of  receiving  and  leading  edge  of  sending 
points  respectively 

X02 

X-?2 

Same  as  above  except  42  1S  trie  trailing 
edge 

CBAR 

c" 

Reference  chord 

KOI  I 

Im(Kdl) 

See  Equation  (  5.1 .4-1) 

KD1R 

Re(Kdl) 

KD2I 

Im(Kd2) 

OUT 

ARG 

See  Equation  (  5.1 .4-2) 

KD2R 

Re(Kd2) 

FMACH 

M 

IN 

ARG 

Mach  number 

E2 

e2 

OUT 

DLM 
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MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

IND 

L 

OUT 

KOS 

KK1 1 

ym) 

KK1R 

Re(KKl) 

See  Section  3.2 

KK2I 

Ik(KK2) 

IN 

KDS 

KK2R 

Re(KK2) 

K10T1 

KZ0T2P 

IN 

DLM 

See  Section  3.2 

Calling  Subroutines  TVOR,  DZY 

Called  Subroutines  and  Common  Blocks 
TKER,  KDS,  DLM 

Equations 


Kdl  =  MCi)  eikrA5/c 

-  K1U2)  e"ikrA5/c  +  Kdlr  L 

(  5. 1.4-1) 

Kd2  =  K2(?i)  eikrA^c 

-  K2U2)  e"lkrAC/C  +  Kd2r  L 

(  5. 1.4-2) 

where 

A?  =  ?2  "  <ii 

K(c)  =  output  from  subroutine  TKER 
Kdir  =  "  (KioUi)  "  Ki o (€2 ) ) 

^d2r  =  “  ^2  (^20(51)  -  ^20(^2 ) ) 

T 1 ,  T2  ,  K10,  K2o  are  output  from  TKER 
L  is  an  option  flag  for  TKER  and  is  either  unity  or  zero. 
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5.1.5  SUBROUTINE  FWMW  (ND,  NE,  SGS,  CGS,  IRB,  AO,  ARB,  XBLE,  XBTE,  YB,  ZB,  XS,  YS. 
ZS,  NAS,  NASB,  KR,  BETA2,  CBAR,  FWZ,  FMY,  MWZ,  MWY,  IF1,  IPRNT,  IBFS) 

Functional  Description 

Given  a  unit  pressure  doublet  this  subroutine  calculates  the  effect  of 
this  doublet  plus  any  contributions  due  to  images,  symmetry  plane  and  ground 
effect  on  a  given  body. 

Input  Output  Variables 


MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

ND 

<5 

Symmetry  flag. 

NE 

e 

Ground  effects  flag. 

SGS 

sinYs 

Sine  of  sending  point  dihedral  angle 

CGS 

cosYs 

Cosine  of  sending  point  dihedral  angle 

IRB 

Number  of  the  receiving  body 

AO 

ao 

Radius  of  the  body 

ARB 

PR 

Array  of  ratios  of  body  axes 

XBLE 

Leading  edge  location  of  slender  body 

element 

XBTE 

Trailing  edge  location  of  slender  body 

el ement 

YB 

yb 

IN 

ARG 

Array  containing  the  y-coordi nates  of 

the  bodies 

ZB 

ZB 

Array  containing  the  z-coordinates  of 

the  bodies 

XS 

1 /4-chord  x-coordinate  of  slender  body 

element 

YS 

y-coordi nate  of  sending  point 

ZS 

z-coordinate  of  sending  point 

NAS 

Number  of  associated  bodies 

NASB 

Array  containing  the  associated  body 

numbers 

KR 

kr 

Reduced  frequency 

BETA2 

e2 

1-M2 

CBAR 

c 

Reference  chord  length 

Fwz 

Force  in  z  direction 

FWY 

F 

OUT 

ARG 

Force  in  y  direction 

mmm 

wy 
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MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

OUT 

ARG 

Moment  in  z  direction 

Moment  in  y  direction 

IFl 

IN 

ARG 

Flag  indicating  the  orientation  of  the 
receiving  body 

IPRINT 

Print  flag 

IBFS 

IN/OUT 

ARG 

Option  flag;  0  to  select  subroutine 

FMZY,  1  to  select  subroutine  FZY2 

Calling  Subroutines  BFM  ,  BFSMAT 

Called  Subroutines  and  Common  Blocks  FMZY,  SUBI,  FZY2 


Equations 


FWZ  =  cos y$  Fzz  -  sinys  Fzy 


+  \  [ii  F  *  +  fi  F  ^  1 

Zs  2  22  y  2y 


bl=range 


+  6  C0SYs  FzzS  +  sin^s  FzyS 


+  2|i‘zFz2I(S)  +  «yF2yI<S,| 

bl  range 


+  £jcos»s  FzzG  +  sin^s  FzyG  +  2  l°zFzzI(G)  +  pyFzyI<G) 

bi-range 

+  6c  jcosys  F22S(G)  +  6invsFzyS(G!  +X|w2F2z!(S(G,)  +  P/zyI(S(G!)| 

bi-range 

FWY  =  same  as  FWZ  replacing  F  with  F  and  replacing  F  with  F 

i.L  yz  zy  yy 


MWZ  -  same  as  FWZ  replacing  Fzz  with  Mzz  and  Fzy  with  Mzy 

MWY  -  same  as  FWZ  replacing  Fzz  with  Myz  and  Fzy  with  Myy 

The  superscripts  on  F  correspond  to 
I  -  Image  points 
S  -  Symmetry  points 
G  -  Ground  effect  points 

bl  range  refers  to  the  bodies  associated  with  the  lifting  surface 
sending  element 
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R1SQX,  XI I JR ,  XIIJI) 

Functional  Description 

Subroutine  IDF1  performs  the  integration  of  the  planar  parts  of  the 
incremental  oscillatory  kernels  according  to  Equation  (B.9)*.  The  result  of 
the  integration  is  the  complex  number  (XIIJR,  XIIJI)  which  is  returned  to 
subroutine  INCRO  via  the  argument  list  of  subroutine  IDF1. 


Input  Output  Variables 


MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

EE 

e 

Semi -width  of  sending  element 

E2 

e2 

ETA01 

y 

V°SV  +  zosini,r 

ZET01 

z 

zoCOSYr  *  yosinYr 

ARE 

Re(A1) 

AIM 

Im(A1 

Coefficients  of  the  parabola  for  the 

BRE 

Re(B-j) 

planar  part  of  kernel  integration  - 

BIM 

Im(B-|) 

IN 

ARG 

See  Eq's  (B.3)  through  (B.5) 

CRE 

Re(C1) 

CIM 

HC}) 

R1SQX 

'1 

-2  -2 
y  +  z 

XIIJR 

■bb 

Real  part  of  integral  T  planar  contri- 

XIIJI 

OUT 

Imaginary  part  of 
integral 

(  bution 

J  See  Eq.  (B.9) 

T.RM2R 

TRM2I 

Real  part 

Imaginary  part 

« 

of  the  second 
>  term  inside  the 

brackets 

TRM3R 

OUT 

I0F1 

Real  part 

of  the  third 

TRM3I 

Imaginary  part 

term  inside  the 

brackets  of 

Eq.  (B.9) 

*  Appendix  B  in  Part  I Vol.  I  of  this  report  (Reference  2) 
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MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

TESTO 

A  test  value: 

Y2  -  52  -  e2 

2  e  |  z  | 

TEST 

i 

Alternate  test  value: 

2  e  z 

__ 

y  +  z  -  e 

I 

ARGT 

OUT 

IDF1 

The  argument  to  the  arctangent  in 

Eq.  (B .6b) 

S 

The  argument  to  the  series  in  Eq.  (B.8) 

SER 

The  sum  of  the  series  in  Eq.  (B.8) 

ALPHA 

See  Eq.  (B.8) 

FUNCT 

F 

See  Eq.  (B.6b) 
and  Eq.  (B.7) 

TRM1R 

| 

Real  part 

of  the  first  term 
y  inside  the  brackets 
of  Eq.  (B.9) 

TRM1I 

Imaginary  part 

> 

Calling  Subroutine  INCRO 

Equations  -  See  Appendix  B  in  Part  I  Volume  I  of  this  report  (Reference  2) 
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S. 1.7  SUBROUTINE  IDF2  (EE,  E2,  ETAOI ,  ZETQ1 ,  A2R,  A2I,  B2R,  B2I,  C2R,  C2I, 
R1SQX,  DIIJR,  DIIOI) 

Functional  Description 

Subroutine  IDF2  performs  the  integration  of  the  nonplanar  parts  of  the  in¬ 
cremental  oscillatory  kernels  according  to  Equations  (B.15)  or  (B.16)*  The 
result  of  the  integration  is  the  complex  number  (DIIJR,  DIIJI)  which  is 
returned  to  subroutine  INCRO  via  the  argument  list  of  subroutine  IDF2. 

Input  Output  Variables 


MNEMONIC 

SOURCE 

DESCRIPTION 

EE 

E2 

ETAOI 

ZET01 

e 

e2 

y 

z 

IN 

ARG 

I  See  Subroutine  IDF1 

A2R 

A2I 

B2R 

B2I 

C2R 

C2I 

Re(A2) 

Im(A2) 

Re(B2) 

Im(B2) 

Re(C2) 

Im(C2) 

1  Coefficients  of  the  parabola  for  the  non- 
(  planar  part  of  the  kernel  integration  - 
(  See  Eq's  (B.12)  -  (B.14) 

R1SQX 

n2 

>2  -2 

y  +  z 

DIIJR 

DIIJI 

TESTO 

TRM2R 

TRM2I 

Re(D2ri 

lM(D2r, 

.) 

.) 

OUT 

IDF2 

Real  part  I  Nonplanar 

°LgiSyapart  contribution 

of  integral  \  See  Eq  (B.15)  or  (B.16) 

Test  value: 

_  2  _  2  2 
y  +.z  -  e 

2  ez 

Real  part  \  of  the  second  term 

/  inside  the  brackets 
Imaginary  part  (  of  Eq  (B.15)  or  same 
)  of  Eq  (B.16) 

*  Appendix  B  in  Part  I  Vol.  I  of  this  report  (Reference  2) 
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MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION  | 

TRM3R 

•< 

Real  part 

> 

of  the  third  term 

inside  the  brackets 

TRM3I 

Imaginary  part 

■ 

TEST 

Alternate  test  value: 

2  e  z 

-2  '  -2 
y  +  z  - 

e‘ 

l 

S 

SER 

OUT 

IDF2 

ALPHA 

a 

See  Subroutine 

IDF1 

FUNCT 

F 

ARGT 

TRM1R 

Real  part 

of  the  first  term 

‘  inside  the  brackets 

TRM1I 

Imaginary  part 

of  Eq.  (B .15) 

or  same  of  Eq. 

(B.16) 

Calling  Subroutine  INCRO 

Equations  -  See  Appendix  B  in  Part  I  Vol .  I  of  this  report  (Reference  2) 


5.1.8  SUBROUTINE  INCRO  (AX,  AY,  AZ,  AX1 ,  AY1 ,  AZ1 ,  AX2,  AY2,  AZ2, 

SGR,  CGR,  SGS ,  COS,  KR,  FL,  BETA,  SDELX,  DELY,  DELR,  DELI, 

IQ,  IR,  NBXS,  NCPNB,  LHS,  NDBLE,  IMG,  NOBI,  IMGS,  USE! ,  USE2, 

USE3,  USE4,  XUSE1 ,  XUSE2,  XUSE3,  XUSE4) 

Functional  Description 

Subroutine  INCRO  prepares  the  arguments  for  the  subroutines  TKER, 

IDF1  and  IDF2  .  It  calls  subroutine  TKER  which  computes  the  incremental 
oscillatory  part  of  the  kernel  K  for  each  receiving-sending  box  combination 

at  the  three  points  of  the  bound  vortex  segment:  at  the  center  (Kc),  at  the 

inboard  point  (K.)  and  at  the  outboard  point  (KQ).  Since  even  a  relatively 
small  case  requires  many  kernel  computations  (e.g.  an  unsteady  100-box  all 
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panel  case,  in  symmetry,  requires  2x100x3x100  =  60,000  kernel  values)  extra 
programming  effort  was  made  to  reduce  the  number  of  kernel  computations. 
Neighboring  strips  have  common  kernels  on  the  common  boundary  lines;  this 
property  is  utilized  in  subroutine  INCRO  for  all  lifting  surface  strips  us 
well  as  their  images  inside  associated  bodies,  by  saving  the  arrays  of  kernels 
for  all  common  boundary  lines  and  using  these  in  the  subsequent  calculations. 

After  the  triplet  of  kernels,  (Kc,  Kj ,  K  )  is  obtained  for  one 
receiving-sending  box  combination,  subroutine  INCRO  computes  the  coefficients 
of  the  parabolas  for  the  numerical  integrations,  done  in  subroutines  IDF1  and 
IDF2. 

Input  Output  Variables 


MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

AX 

X 

x  -  ? 

where  (x,  y,  z) 

AY 

y 

y  -  n 

define  the  receiving 

AZ 

z 

Z  -  c 

point,  and  ?,  n,  ?  are 

AX1 

x  -  5l 

the  'center'  sending 

AY1 

IN 

ARG 

y  -  n, 

point  coordinates. 

AZ1 

z  "  ci 

?j,  n,,  ?,  are  the 

AX2 

X“52 

'inboard'  sending  point 

AY2 

y  -  n2 

coordinates  and 

AZ2 

2  -  s 

?2,  n2>  s2  are  the 

'outboard'  sending 

point  coordinates 
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MNEMONIC  .  SYMBOL 


IN/OUT 

IN 


SOURCE 


DESCRIPTION 


SGR 

sinyr 

where  'r*  is  the 

CGR 

cosyr 

index  of  the 

SGS 

CGS 

sinYs 

c°sys 

receiving  point  and 
's'  is  the  index  of 

the  sending  point 

KR 

kr 

Reduced 

frequency 

FL 

c/2 

Reference  semi -chord 

BETA 

e 

ft  «  n- 

-  M2 

SDELX 

6xs 

Average  box  length 1  of  sending 

DELY 

2es 

Box  width  j  Point  's' 

DELR 

iJR' 

OUT 

ARG 

Real  part  'lof  the  oscilla¬ 

tory  contri- 

Imaginary  partj  but,on  t0  the 

nonplanar  down- 

wash  factor 

10 

An  index  running  from  1  through 
NC^,  where  NC^  is  the 

no.  of  chordwise  boxes  in 

the  sending  panel  p 

IR 

Index  of  sending  point 

NBXS 

The  no. 

of  boxes  on  the  panel 

in  which  the  sending  point 

lies  +  the  total  no.  of  boxes 

NCPNB 

IN 

_ 

of  the  preceding  panels 

The  no.  of  boxes  in  the  first 

strip  of  the  panel  in  which 
the  sending  point  lies  +  the 
total  no.  of  boxes  of  the 
preceding  panels 

MNEMONIC 


SYMBOL 


IN/OUT  SOURCE 


DESCRIPTION 


Flag  activated  by  6  (input)  for 
the  contribution  of  the  first 
symmetry  plane  effect 

Flag  activated  by  e  (input)  for 
the  contribution  of  the  second 
symmetry  plane  effect 

Flag  activated  by  the  presence 
of  image  points  inside 
associated  bodies  for  sending 
panel 

Sequence  number  of  the  body  in 
which  the  image  of  the  sending 
point  lies 

IMGS  =  IMG  whenever  kernels  are 
saved  for  image  sending  points; 
IMGS  =  0  otherwise 


of  the  steady 
contribution 
to  the  kernel 


MNEMONIC  SYMBOL  (  IN/OUT  SOURCE  DESCRIPTION 


COUNT 


Labeled 

Common 


Utility  arrays  for  use  of  the 
planar  kernel  values  which 
are  saved  for  future  use 


Utility  arrays  for  use  of  the 
nonplanar  kernel  values  which 
are  saved  for  future  use 


Mach  Number 
Reference  sami-chord 
e  =  0.00001 


See  Subroutine  IDF1 


An  internal  flag  to  select 
logic  of  subroutine 


0 


MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

XO 

Differences  of  the  x,  y  and  z 

YO 

OUT 

INCRO 

coordinates  of  the  current 

receiving  and  sending  points 

(see  AX,  AY,  AZ, 

AX1,  AY1 ,  AZ1 

and  AX2,  AY2,  AZ2  above) 

ATI 

mm 

Argument 

■Mi 

IN 

List  of 

See  Appendix  A  in 

Part  I,  Vol  I  of 

AT  2 

TKER 

this  report 

DKRC 

for  the  center 

Planar  kernel 

point  of  the 

DKIC 

bound  vortex 

segment 

XKRC 

Nonplanar  kernel 

XKIC 

DKRI 

Planar  kernel 

for  the  inboard 

KDII 

point  of  the 

IN 

INCRO 

XKRI 

Nonplanar  kernel 

bound  vortex 

XKII 

segment 

OKRO 

Planar  kernel 

for  the  outboard 

DKIO 

point  of  the 

XKRO 

Nonplanar  kerr.°l 

bound  vortex 

XKIO 

segment 

Index  for  the  selection  of  the 

JO 

proper  utility  array  in  which 

the  kernels  are  saved 

XMULT 

ax  '/8if 

J 

ETA01 

ZETOl 

R1SQX 

ARE 


OUT 


INCRO 


See  Subroutine  IDF1 


SW 

i  N / Ob  f 

SOURCE 

D!s;ription 

AIM 

,;ur 

INCRO 

BRE 

BIM 

CRE 

CIM 

See  Subroutine  !UF1 

XI I  JR 

Argument 

List  of 

XIIJI 

IDF1 

A2R 

A2I 

B2R 

OUT 

INCRO 

B2I 

See  Subroutine  1DF2 

C2R 

C2I 

DIIJR 

IN 

Argument 

List  of 

DIIJI 

IDF2 

DELR 

DELI 

— . ~J 

iTU^s) 

OUT 

INCRO 

AD  =  D,  +  Do 
rs  Irs  2-s 

See  Appendix  B  in  Part  J,  Vol 
of  this  report 

Calling  Subroutine  SUBP 


Called  Subr-.ig* ' res  and  Comiion  Blocks 

Subroutines  *KLP.  T0F1  and  I0F2,  and  the  Labeled  Common  Blocks  flLM  and 
KDS . 


Equations  :  ppendi/.  (:  n  ^art  I  Vol.  I  of  this  report  (Ref.rence  2) 


5.1.9  SUBROUTINE  READD  /Q,  N,  NTAPE) 
Functional  Description 


This  routine  reads  the  complex  array  D  which  is  N  complex  words 
long  from  the  I/O  unit  NTAPE. 

Input  Output  Variables 


MNEMONICS 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

D 

OUT 

ARG 

Complex  array  read 

N 

IN 

ARG 

Number  of  complex  words  in 

array  D 

NTAPE 

IN 

_ 

ARG 

Input  unit  number 

Calling  Subroutines  OUTCOR,  SB 

5.1.10  SUBROUTINE  RWREC  (I FLAG,  NTAP,  A,  NCWORD,  NUMBR) 

Functional  Description 

This  routine  is  used  to  read  and  write  the  array  A  which  is  NCWORD 
complex  words  long  on  the  I/O  unit  NTAP. 

Input-Output  Variables 


MNEMONIC 

SYMBOL 

. 

IN/OUT 

SOURCE 

DESCRIPTION 

I  FLAG 

IN 

Read-write  flag 
=  0  write  A 

~  1  y'pa/J  A 

=  2  write  NUMBR,  A 
=  3  read  NUMBR,  A 

'  NTAP 

IN 

ARG 

I/O  unit  number  to  use 

A 

IN/OUT 

Complex  array  to  be  read  or  written 

IN 

Length  of  array  A 

NUMBR 

IN/OUT 

The  number  to  precede  array  A  if 
desired;  zero  otherwise 

Calling  Subroutines  BFM,  BFSMAT,  MAIN,  OUTCOR,  SB 


5.1.11  SUBROUTINE  SNPDF  (SL,  CL,  TL,  5GS,  CGS,  SGR,  CGR, 

XO,  YO,  ZO,  EE,  EIJ,  BETA,  C V) 

Functional  Description 

Subroutine  SNPDF  computes  the  steady  downwash  factors  for  one 

receiving-sending  point  combination  at  a  time  according  to  Eq.  (C .35b) , 

* 

Appendix  C  .  The  result,  DIJ,  is  returned  to  the  calling  subroutine  via  the 
argument  list  of  subroutine  SNPDF. 

Input  Output  Variables 


MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

m 

sinx  where  x  is  the 

EH 

cosx  sweep  angle  of  the 

TL 

IN 

ARG 

tanx  1 /4-chord  line  of 

the  sending  box,  i.e.  the  sweep 
angle  of  the  bound  vortex 

SGS 

CGS 

SGR 

See  Subroutine 

CGR 

INCRO 

XO 

X 

AX 

YO 

y 

AY 

ZO 

i 

AZ 

EE 

e 

DIJ 

D(s) 

rs 

OUT 

SNPDF 

Steady  contribution  to  the 

downwash  factor 

See  Eq.  (C .35b) 

BETA 

0 

IN 

ARG 

See  Subroutine  INCRO 

CV 

*xs 

_ 

*  Appendix  C  in  Part  IVol.  I  of  this  report  (Reference  2) 
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MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

RIMAG 

Ri 

See  Eq.  (C .  1 4) 

ROMAG 

*0 

See  Eq.  (C .16) 

DB2 

See  Eq.  (C.30) 

VBY 

V 

VBZ 

Vbz 

IN 

SNPDF 

VIZ 

viy 

See  Eq.  (C.33b) 

VIZ 

Viz 

VOY 

Voy 

VOZ 

Voz 

■■ 

w 

See  Eq.  (C.34b) 

Calling  Subroutines  SUBP  and  TV0R 


Equations  See  Appendix  C  in  Part  I,  Vol .  I  of  this  report 

(Reference  2) 

5.1.12  SUBROUTINE  SUBI  (DA,  DZB,  DYB,  PAR,  DETA,  DZETA,  DC GAM,  DSGAM, 

DEE,  DXI,  TL,  DETAI,  DZETAI ,  DCGAMI,  DSGAMI ,  DEE I,  DTLAMI,  DMUY, 
DMUZ,  INFL,  IQTFL 


Functional  Description 

Subroutine  SUBI  has  a  dual  role:  depending  on  the  setting  of  the 
flag  INFL,  subroutine  SUBI  either  computes  the  arguments  necessary  for  the 
calculation  of  the  influence  coefficient-contribution  of  lifting  surface 
image  points  inside  associated  bodies  (INFL  =  0),  or  it  calculates  the 
uz  for  one  body  element  according  to  Eqs.(  5.1.12-21  and  -22) 
respectively  (INFL  =1). 


75 


Input  Output  Variables 


MNEMONIC 

SYMBOL 

DA 

a 

DZB 

zc 

DYB 

DAR 

AR 

DETA 

ns 

DZETA 

5s 

DCGAM 

cosys 

DSGAM 

sinys 

DEE 

es 

DXI 

TL 

tanx 

DETAI 

nl 

DZETAI 

Si 

DCGAMI 

COSyl 

DSGAMI 

sinyl 

DEE  I 

el 

DTLAMI 

tan  XI 

DMUY 

'V 

uy 

DMUZ 

'V 

IN/OUT  SOURCE 


DESCRIPTION 


Sending  point 
data  (on  panels) 


y-coordinate  of(Eq.  5.1.12-1 

image  nv.  _iq\ 

point  or  ' 

z-coordinate  of(Eq.  5.1.12-2 

or  -20) 

cosine  of  the  (Eq.  5.1.12-3 

dihedral  angle 

y  and  -4) 

of  the  image 

of  the  bound 

sine  vortex  plane 


Semi -width  of  image  of  sending 
strip  (Eq.  5.1 .12-6) 

Tangent  of  the  sweep  angle  of  the 
image  of  bound  vortex; 

Eq. (5.1. 12-5) 

Eq.  (5. 1.12-21) 

Eq.  (5.1.12-22) 
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MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

INFL 

IOUTFL 

f 

IN/ OUT 

ARG 

'In' -flag  -  0  or  1 

'Out '-flag  -  zero  if  the  image 
of  bound  vortex  lies  beyond  the 
symmetry  plane  of  the  body  - 

1  otherwise 

PSQR 

P 

Intermediate 

Eq .  (5,1 .12-10) 

YCBAR 

variables  used 

(5.1.12-14) 

ZCBAR 

in  SUBI 

(  5.1.12-15) 

ABAR 

a 

(5.1 .12-11) 

RH012 

A 

IN 

SUBI 

(5.1.12-12) 

RH022 

2 

P2 

(5.1.12-13) 

ETAI1 

nil 

(5.1.12-7) 

ETAI2 

nI2 

(5.1.12-8) 

ZETI1 

dl 

(5.1.12-9) 

ZETI2 

CI2 

(5.1.12-10) 

RH02 

2 

P 

(5.1.12-23) 

ELLIPS 

See  Eq. (5.1 .12 

-24  or  -26) 

TEST 

Test  value;(Eq 

.  5.1.12-25 
or  -27) 

Calling  Subroutines  SUPP  and  FWMW 
Equations 

Option  1.  -  INFL  =  0 


nl  =  (nil  +  nI2)/2 


d  =  (dl  +  cI2)/2 
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(5.1.12-1) 

(5.1.12-2) 


cosrl  -  -  *1^11 

(5.1.12-3) 

slnxl  =  -  d2i!;1el1 

(5.1.12-4) 

tanxl  =  d22?I-sa 

(5  ..1 .12-5) 

el  =\/(nI2  -  nil)2  +  (cI2  -  dl)2  /2.0 

(5.1.12-6) 

where 

-2 

(5.1.12-7) 

nil  =  yc  +  p-2(nl-yc) 

-2 

(5.1.12-8) 

nI2  =  yc  +  ^>(n2-yc) 

-2 

(5.1.12-9) 

?I1  =  z  +  ^-2(?1“Z  ) 

Co  C 

1 

-2 

(5.1.12-10) 

5 12  =  ZC  +  ^z(?2"zc^ 

a  =  ^sin2e  +  AR2cos2e)3/2 

(5.1.12-11) 

pi2  =  (nl-yc)2  +  (cl*5c)2 

(5.1.12-12) 

p22  =  (n2-yc)2  +  (?2-zc)2 

(5.1.12-13) 

yr  *  a(l-AR2)cos3e  +  YB 

V 

(5.1.12-14) 

zc  ‘  a  TT  S1n3e  +  ZB 

(5.1.12-15) 

cos  e  =  AR 


sine  = 
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P  =V  {n-YB)2AR2  +  (c-ZB)2 
and  where 


(5J.12-18) 


nl  =  n-ecosy 
n2  =  n+ecosy 
d  =  c-esiny 
?2  =  c+esiny 


=  K  +  h-ri^r) 


(5.1.12-19) 

(5.1.12-20) 


uy  =  -  ^H-sinY[(n-yc)2  -  (c-zc)2]  +  2cosy(n-yc)U-zc)}  (5.1.12-21) 
2 

"z  =  -  4  (-2s1nY(n-yc)(c-zc)  -  cosY[(n-yc)2  -  (5.1.12-22) 


and 


TRM8  =  a/ [a  n/ar2-!  +  (zl-z  J]2  +  (yl-yj2 

w  L 

If  AR<1 

ELLIPS  =  2a  (5.1.12-26) 

TEST  =  TRMC  +  TRMD  (5.1.12-27) 

where 

TRMC  =  >/  [a  v/T-AR2"  -  (yl-yr)]2  +  (zl-zj2 

v  V/ 

and 

TRMD  =  V/Cav/T  -AR2  +  (yl-yj]2  +  ( zl-z  )2 

V  V 

and  where  yl  and  zl  are  the  image  point  coordinates:  nl>  t,l  for  the  'center' 
sending  point,  nil.  ?I1  for  the  'inboard',  and  nI2,  c 12  for  the  'outboard' 
sending  point. 

5.1.15  SUBROUTINE  TIME  (NN) 

Functional  Description 

Subroutine  TIME  is  a  utility  subroutine  which  is  used  to  printing 
out  the  CPU-  and  I/O-time  spent,  on  tU  computer  either  since  the  beginning 
of  the  run  (setting:  NN  =  0  ),  or  since  the  last  call  to  TIME  (NN  =  1  ). 

In  program  N5KA  subroutine  TIME  is  called  from  MAIN  for  each  major  subroutine, 
so  that  a  breakdown  on  the  computing  time  per  major  computations  is  obtained 
for  the  case. 

5.1.14  SUBROUTINE  TKER  (XO,  YO,  ZO,  KR,  BR,  SGR,  CGR,  SGS,  CGS,  T1 ,  T2,  M) 
Functional  Description 

Subroutine  TKER  computes  the  incremental  oscillatory  kernel,  K, 
for  one  receiving-sending  box  combination  at  a  time.  The  result  is  obtained 
in  six  components  and  is  returned  to  the  calling  program,  subroutine  INCR  , 
via  the  labeled  common  block  OLM  (see  Sec.  3.2).  Subroutine  TKER  also  com- 

A,  A/ 

putes  the  total  planar  and  nonplanar  kernels  K-j  and  Kg  whenever  the  flag 
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IND  is  set  to  0;  this  option  is  exercised  by  the  calling  subroutine  FLLD. 
The  flag  IND,  as  well  as  the  total  kernels,  are  transmitted  via  the  labeled 
common  block  KDS. 

Input  Output  Variables 


MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

AX  for  center 

XO 

X 

i 

AX1  for  inboard 

AX2  for  outboard 

AY  for  center 

YO 

y 

AY1  for  inboard  point  of 

AY2  for  outboard  sending 

ZO 

z 

IN 

ARG 

box;  see 

AZ  for  center  subroutine 

AZ1  for  inboard  INCRO 

AZ2  for  outboard 

KR 

kr 

BR 

br,  c/2 

SGR 

See  Subroutine  INCRO 

CGR 

SGS 

CGS 

T1 

Ti 

OUT 

TKER 

T1  =  cos^r  “  Ys} 

T2 

T2 

T2  =  C(z0C0SYr  "  Vinyr) 
x(zocosYs-y0sinYs)]/(br/10)2 

M 

M 

IN 

ARG 

Mach  Number 
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MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

-2  _2 

R1 

r 

r  =  y  +  z 

T2P 

* 

T2 

See  Eq.  (A.8)* 

2 

BETA2 

$ 

BIGR 

R 

See  Eq.  (A. 14) 

K1 

kl 

Eq.  (A. 12) 

MU1 

Eq.  (A. 11) 

IOUR 

Re(I0) 

IN 

TKER 

Eq.  (A. 30) 

IOUI 

mido) 

JOUR 

Re(J0) 

Eq.  (A. 31 ) 

JOUI 

Im(J0) 

I1UR 

Red^ 

Eq.  (A. 25) 

nui 

Im(I-j) 

I2UR3 

3Re(I2) 

Eq.  (A. 27) 

I2UI3 

3Im(I2) 

K10T1 

k!5)ti 

Planar  part 

of  the 

steady 

K20T2P 

|/(S)t  * 
^2  '2 

OUT 

TKER 

Nonplanar  part 

kernel 

K1RT1 

K1IT1 

Planar  part  of 

the  unsteady 

►  kernel ; 

K2RT2P 

Nonplanar  part  of 

see  Sec.  3.2 

K2TI2P 

■ 

t  Appendix  A  in  Part  I  Vol .  I  of  this  report  (Reference  2) 
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MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

KD1R 

KD1J 

KD2R 

KD2I 

Re(K1) 

Im( ) 

Re(K2) 

Im(K2) 

OUT 

TKER 

See  Sec.  3.2 

Calling  Subroutines  INCRO  and  FLLD. 


Common  Blocks 

Common  DLM  and  Common  KDS 


Equations  See  Appendix  A  in  Part  I  Vol .  I  of  this  report  (Reference  2) 

5.1.15  SUBROUTINE  TVOR  (  SL1 ,  CL1 ,  TL1 ,  SL2,  CL2,  TL2,  SGS,  CGS,  SGR, 
XOl,  X02,  YO,  ZO,  E,  BETA,  CBAR,  FMACH,  KR,  BRE,  BIM) 

Functional  Description 

This  routine  calculates  the  normalwash  at  a  point  (x,  y,  z)  of 
surface  dihedral  Yr»  due  to  a  trapezoidal  unsteady  vortex  ring  of  unit 
strength. 


Input  Output  Variables 


MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

SL1 

CL1 

sinx-| 

COSA-j 

A-.  is  the  sweep  angle  of 
leading  edge  of  box 

TL1 

tanx-| 

SL2 

sinx2 

IN 

ARG 

CL2 

cosa2 

a2  is  the  sweep  angle  of 

TL2 

tanx2 

trailing  edge  of  box 

SGS 

sinYs 

Sine  of  dihedral  angle  at  sending 
point 
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MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

CGS 

SGR 

CGR 

XOl 

X02 

YO 

ZO 

E 

BETA 

CBAR 

FMACH 

KR 

c°Sys 

slnYr 

COSYr 

x-?i 

X+52 

y-n 

y-5 

e 

B 

c 

M 

kr 

IN 

ARG 

-  — -  -  - 

Cosine  of  dihedral  angle  at 
sending  point 

Sine  of  dihedral  angle  of 
receiving  point 

Cosine  of  dihedral  angle  of 
receiving  point 

Distance  in  x-direction  from  receiving 
point  to  leading  edge  of  box 
?2  indicates  trailing  edge  of  box 

Distance  in  y-di recti  on  from  receiving 
point  to  center  of  box 

Distance  in  z-di recti  on  from  receiving 
point  to  center  of  box 

\/l  -  M2 

Length  of  reference  chord 

Mach  Number 

Reduced  frequency 

BRE 

Re(B) 

OUT 

ARG 

Real  part  of  the 
normalwash  See 

Fn 

BIM 

Im(B) 

OUT 

ARG 

5.1.15-1) 

Imaginary  part  of 
the  normalwash 

Calling  Subroutines  DZY 


Called  Subroutines  and  common  blocks  SNPDF,  FLLD 


Equations 


B 


=  BS 
~  2eAx 


AB 

48ti 


(5.1.15-1) 


where 


BS  -  D^  ^[(^i*  XQ-j ) ,  (y^ >  Yj Yq >6 >Ax) ] 
^[(A2»Xq2) ,  >^o ’^o*e }^»A^) ] 
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)<s>  -  i 


s  calculated  in  subroutine  SNPDF 


Ax  =  ?2  -  Cj 


ri  and  r0  1  T 

a B"  _  aKj  +  4ak : .  +  ak  . 

i  c  o 


r.  or  r  <  -r 
ri  0  4T 


AB  =  6aK. 


r_  < 


c  -  ¥ 


AB"  =  3aK .  +  3aK. 

i  ao 


where 

_  AKdl  +  AKd2 

ai\j  _  — 2~  — 5“ 

r  r 

ri2  =  ^yo  +  e  C0Sys)2  +  (z0  +  e  cosys)2 

2  2  2 
r  -y  +z 
c  o 

ro2  =  (y0  "  e  cos^s^2  +  (z0  "  e  C0SYS)2 


aKj i  and  aK(j2  are  calculated  by  FLLD 


5.1.16  SUBROUTINE  WRTFMF  (IUNT,  NBE,  FZ,  FY,  EMZ,  EMY 


Functional.  Description 


This  routine  writes  the  NBE  elements  of  the  complex  arrays  FZ,  FY, 
EMZ,  and  EMY  with  a  format  on  I/O  unit  IUNT. 
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Input  Output  Variables 


SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

IUNT 

I/O  unit  on  which  the  data  is 

to  be  written 

NBE 

Number  of  complex  words  in 
each  array 

FZ 

IN 

ARG 

Array  to  be  written 

FY 

Array  to  be  written 

EMZ 

Array  to  be  written 

EMY 

Array  to  be  written 

Calling  Subroutine  BFM 

5.1.17  SUBROUTINE  WRTFNU  (IUNT,  MODE,  NBE,  FZ,  FY,  EMZ,  EMY) 

Functional  Description 

This  routine  writes  MODE  and  the  NBE  elements  of  the  complex  arrays  FZ, 
FY,  EMZ,  EMY  without  a  format  on  the  I/O  unit  IUNT. 


Input  Output  Variables 


MNEMONIC 

SYMBOL 

HP! 

SOURCE 

DESCRIPTION 

IUNT 

I/O  unit  on  which  the  data  is 

to  be  written 

MODE 

The  mode  number  to  be  written 

FZ 

IN 

ARG 

Array  to  be  written 

FY 

Array  to  be  written 

EMZ 

Array  to  be  written 

EMY 

Array  to  be  written 

---  -  ■ 

Calling  Subroutine  BFM 
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5.2  Segment  2 


5.2.1  SUBROUTINE  DATA 
Functional  Description 

Subroutine  DATA  reads  the  x-,  y-  and  z-coordinates  that  define  each 
panel,  interference  body  and  slender  body  for  the  case  considered.  It  also 
readsthe  fractional  chordwise  and  spanwise  divisions  for  each  panel  and 
calculates  arrays  of  geometry  that  define  the  sending-  and  receiving  control 
points  on  the  lifting  surfaces.  In  the  case  when  bodies  are  also  present, 
similar  geometry  arrays  are  calculated  for  all  interference  body  elements, 
and  all  slender  body  elements.  These  basic  data  arrays  are  saved  in  the 
Blank  Common  Block  for  use  in  the  subsequent  calculations. 

Input  Output  Variables 
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MNEMONIC  SYMBOL 


ZIN(IOO)  z/P 


NASB ( 200 ) 


NAS(IOO) 


NCARAY(IOO)  nc 


NSARAY(IOO)  ns 


NBARAY(IOO)  nba 


TH(50) 


TAU(50) 


GMA(IOO) 


X( 500 ) 


XII (500) 


XI2(500) 


DELX(500) 


IN/OUT 


*xPc,s 

and 


SOURCE 


DESCRIPTION 


z-coordinate  of  inboard  edge  of 
panel 

Array  of  all  associated  bodies  for 
all  panels 

Array  of  the  number  of  associated 
bodies  for  all  panels 

Array  of  the  number  of  chordwise 
boxes  for  all  panels 

Array  of  the  number  of  spanwise 
strips  for  all  panels 

nba„  =  >  nc.  ns. 

P  i=l  1  1 

where  p  is  the  panel  number 


Fractional  chordwise  divisions 
card  for  Pane'  P 

Fractional  spanwise  divisions 
for  panel  p 


Dihedral  angle  of  panel  p 

3/4-chord  x-coordinate  of  all 
boxes  and  x-coordinate  of  inter¬ 
ference  body  section  midpoints  - 
see  Eqs.(  5/2. 1-1  and-2) 

1 /4-chord  coordinates  of  inboa rd- 
and 

outboard  edge  of  panel  boxes  -  see 
Eqs . (  5. 2. 1-3  and  -4) 

Average  chordlength  of  boxes  and  ax 
of  interference  body  sections  - 
see  Eqs . (  5. 2. 1-7  and  -8) 
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MNEMONIC 


SYMBOL 


IN/OUT 


SOURCE 


DESCRIPTION 


YS(210) 
ZS(210) 
DYS (200) 

!  nzs{2oo) 

EE (200) 
CS(200) 

SG(200) 

CG( 200 ) 

XI J (200) 

GMAR(200) 

XLAM(500) 

ZC 

YX 

RAD 
AR 

NBE 
NSBE 
NZY 


yps 

y-coordinate  of  centerline  of 
panel  strips  -  Eq.  (5.2. ; .9) 

2Ps 

z-coordi nates  of  centerline  of 
panel  strips  -  Eq. (5.2.1-10) 

Ays 

((5.2.1-12)  and 

See  Eqs 

Azs 

((5.2.1-13) 

es 

Half-width  of  panel  strips  - 
Eq. (5.2. 1-11) 

s 

OUT 

DATA 

Average  chord  of  panel  strips  - 
Eq. (5.2.1-16) 

sine 

of  dihedral  angle 

and 

for  strips  - 

cosine 

Eqs.  (5.2.1-14  and  -15) 

x-coc.  dinate  of  leading  edge  of 
strip  centerlines 

Ys 

Dihedre*  angle  of  strips  in 

radians 

tanAc,s 

Tangent  o-  the  sweep  angle  of  the 

1/4-chordl i  *e  of  all  lifting 

surface 

boxe*  -  see  Eq . (5 .2 .1 -13a) 

zc 

z- coordinate  <c  body  centerline 

yc 

y-coordinate  of  uody  centerline 

a 

Average  half-width  of  body 

Cross-secfctonal  aspect  ratio  of 

IN 

card 

body 

No.  of  interference  body  elements 

No.  of  slender  body  elements 

z-y  oriectation  flag 
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MNEMONIC  SYMBOL 


IN/OUT  SOURCE 


NBEA(10,2) 


NSBEA(IO) 


XII(IOO) 


RI(IOO) 


XIS(IOO) 


RS(IOO) 


DESCRIPTION 


Interference  body  element 
radius  flag  (Rl-flag) 

Slender  body  element  radius 
flag  (RS-flag) 

Number  of  An  -  A?  pairs  for  body 
Number  of  ei's  for  body 

Number  of  62 1 s  for  body 


Array  of  z-coordi nates  of 
body  centerlines 

Array  of  y-coordi nates  of 
body  centerlines 

Array  of  average  characteristic 
half  widths. 

Array  of  aspect  ratios, 

1.  No.  of  interference 
body  elements 

2.  z-y  flags. 

No.  of  slender  body  elements  - 
for  all  bodies 


x-coordi nates  of  interference 
body  element  endpoints 

Radii  of  interference  body 
element  endpoints 

x-coordinates  of  slender  body 
element  endpoints 

Radii  of  slender  body  element 
endpoints 


MNEMONIC 
TH1 (  24  ) 


SYMBOL 


IN/OUT 


SOURCE 


DESCRIPTION 


01 


card 


Angular  orientation  of  point  y 
on  body,  first  set 


TH2  (  24) 


02 


THIA(IOO) 


Angular  orientation  of  point  u 
on  body,  second  set 

Array  of  ei^'s  for  all  bodies 


OUT 


DATA 


TH2A(100) 


card 


IFLA( 30 ,2 ) 


OUT 


DATA 


Array  of  02  's  for  all  bodies 

_ 21 _ 

First  and  last  interference 

body  elements  associated 

with  01 

v 


1.  Array  of  'first'  inter¬ 
ference  body  elements, 

2.  Array  of  'last'  interference 
body  elements. 


NFL(IO) 


NCDl 

CD2 

CD3 


Array  of  the  number  of  pairs 
of  'first-  and  last'  elements  - 
for  all  bodies 

Panel  no.  for  An,  Ac 
y-shift  of  panel 
z-shift  of  panel 


IN 


card 


NCD4 

CD5 

CD6 

NQB0DY(5O) 

NPANEL(50) 


another  set  of  the  above 
3  items 

Array  of  the  body  no.'s. 

Array  of  the  panel  no's  for 
the  An,  Ac  shift 


MNEMONIC 

SYMBOL 

DETA(50) 

An 

DZET (50) 

A? 

RIA(IOO) 

XLE(IO) 

XTE(IO) 

XISK200) 

«sit 

XIS2(200) 

IN/OUT  SOURCE 


A0(200)  a 


A0P(200)  a. 


DESCRIPTION 


Array  of  the  y-shifts  for  panels 

Array  of  the  z-shifts  for  panels 

Array  of  the  radii  of  all 
interference  body  element 
midpoints 

Eq . (  5.2.1-17)  when  NBE/0. 

If  NBE=0,  but 
NSBE/0,  the 
Eq.(  5.2.1' 18)  slender  body 
leading-  and 

tn 1  ng  x-coordi  nates  are  used 

Array  if  the  x-coordinates  of 
leadinc.-  and  trailing  edges 
(respnc -Ivety)  of  all  slender 
body  element.;;  see  Eq's 
(5.2.1-20  and  -21) 

Array  of  the  average  radii  of 
slender  i.-ody  elements;  see 
Eq .  (  5.2.1-22) 


Array  of  tht  first  derivatives 
of  a.  for  all  slender  body 

elements;  see  Eq„{  5.2.1-23) 


Calling  Subroutine:  MAIN 
Equations 

3/4-chord  x-coordinate  for  all  boxes  of  panel: 


[(  fa  +  i  '  h~  h  +  xi)+(x3 -  V3 

+  (?9c  +  f  0c+l)(>2  ‘  xl)  +  ^  (  5‘2*1-D 
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x-coordinate  of  interference-body  section  midpoint: 

T  _  cIk+l  +  ?Ik 
X1k  l 


(5.2. 1-2) 


1/4-chord  x-coordinate  of  inboard  edge  of  panel  boxes: 
Spic,s  -  ^s  BR  +  CR 


(5. 2. 1-3) 


1/4-chord  x-coordinate  of  outboard  edge  of  panel  boxes: 
{P2c,s  *  Vl  BR  +  CR 


(5.2. 1-4) 


where 

BR  -  <?  »C  +  T  ec+l)(X4  -  *3  -  \  *  Xl>  +  (X3  -  V 


■  -  *,> +  *, 


leading 

Xrcoordinate  of  and  edges  of  interference-body  sections: 
trailing 

?I1k  “  5lk 
5l2k  =  ?Ik+l 


(5.2. 1-5) 

(5.2. 1-6) 


Average  chord-length  of  panel  boxes: 


_  Vl 


£<ec+i  ’  6c){X4  '  h 


+  (0c+i  "  0c><  h  ~ 


(5. 2. 1-7) 


Ax  for  interference  body-sections: 

A XI.  =  ?I2.  -  511 u 


(5.2. 1-8) 


y-coordinate  of  centerline  of  panel  strips: 


Ts  +  Vl 


(y2-  V+Vl 


(5.2.1-9) 
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(5.2.1-10) 

(5.2.1-11) 

(5.2.1-12) 
(5.2.1-13) 
(5 .2.1 -13a) 

(5.2.1-14) 

(5.2.1-15) 

(5.2.1-16) 


Leading  and  trailing  edge  coordinates  of  bodies 

XL.E.  =  5llk  »  k  =  1  (5.2.1-17) 

XT.E.  =  5l2k  »  k  =  NBE  (5.2.1-18) 

x-coord'inate  of  slender  body  section  midpoint 

xst  =  (?Tt  +  5st+i)/2  (5.2.1-19) 

leading 

x -coordinates  of  and  edges  of  slender  body  sections 
trailing 

5Slt  =  5St  (5.2.1-20) 

?S2t  =  5St+1  (5.2.1-21) 
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(5.2.1-22) 


Average  radius  of  slender  body  sections: 


RSt  +  RS 


t+1 


1 


First  derivative  of  aQ: 


RS 


t+1 


-  RS. 


where 

axst  =  ?S2t  -  £Slt 


(5.2.1-23) 


(5.2.1-24) 
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5.3  Segment  3 

5.3.1  SUBROUTINE  GEND  (NPR1NT,  NTAPE v  WORK) 

Functional  Description 

Subroutine  GEND  generates  the  nine  submatrices  of  the  [DT]  matrix,  then 
assembles  the  [DT]  matrix  and  writes  it  on  logical  tape  unit  no.  ITP8  in  row 
order.  The  submatrices  generated  by  GEND,  and  the  major  subroutines  that 
compute  these,  are  shown  in  the  diagrams  below. 

DT  Matrix  Structure 


Major  Subroutines  Generating  DT 


N‘ 

•p 

Subroutine 

DPPS 

Subroutine 

DZPY 

Subroutine 

DYPZ 

N‘ 

N‘ 

7 

fY 

Subroutine 

DPZY 

- - NTP - - 

— NTZ - 

- — NTY - *- 
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General  Flow  Chart  -  Subroutine  GEND 


DYPZ 


DZPY 

assembles  rows 
of  [DZP], [DZZ] 
and  [DZY] 


assembles  rows 
of  [DYP] >[DYZ] 
and  [DYY] 


Input  Output  Variables 


MNEMONIC  SYMBOL  IN/OUT  SOURCE 


IFIRST 


HAST 


YFLAG 


DESCRIPTION 


Print  flag  for  DT  matrix 
Logical  tape  unit  array 
Complex  work  array 


One  row  of  the  complex  DT  matrix  - 
also  used  for  temporary  storage  of  rows 
of  all  submatrices  of  DT,  except  those 
of  DPZ  and  DPY, 

One  row  of  the  complex  DPZ  matrix 
One  row  of  the  complex  DPY  matrix 


Floating  point  variable  for  6 
the  integer  input  e 
Reference  chord 


Do-loop  delimiters  for  the  number 
of  rows  of  a  particular  submatrix  to 
be  computed  within  the  do  loop 

Do-loop  delimiters  for  the  number  of 
elements  in  one  row  of  a  particular  sub¬ 
matrix  to  be  computed  within  do-loop. 
Running  index  of  the  section  number  of 
the  receiving  body,  KB,  with  01 
distribution. 

The  number  of  sections  of  body  KB  with 
ei  distribution, 
z-y  flag  of  body  KB 

z-y  flat  of  body  preceding  the  present 
body  KB;  0  for  KB=1 

Sequence  number  of  the  first  element  in 
the  current  section  of  body  KB  with 
ei  distribution 

Sequence  no.  of  the  last  element  in  the 
current  section  of  body  KB  with 
01  distribution 

Internal  flag  selecting  correct  do-loop 
for  bodies 
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Calling  Subroutine 


MAIN 


Called  Subroutines  and  Common  Blocks 

Subroutines  DPPS,  DPZY,  DZYP  and  DYPZ,  and  the  Blank  Common  Block. 
Semi-Detailed  Block  Diagram  of  Subrouting  SEND 


Legend 


Receiving  indices  1 r 1 

I  - 

r  point 

K  - 

r  panel 

KS  - 

r  strip  (or  equivalent 
for  body) 

KB  - 

r  body  (if  any) 

KT  - 

index  of  the  array  of 
the  'first-and-last 

elements'  for  ei 

IZ  - 

r  body  element 

Sending 

indices  's' 

J  - 

s  point 

L  - 

s  panel 

LS  - 

s  strip  (or  equivalent 
for  body) 

LB  - 

s  body 

NTZ  =  0 


NZYKB=1 


I1=NB0X+1 

I2=NB0X+NTZ 

KB=1 ,  1=0 
KT=1 


Write  one 
\row  of  DPYy 
\on  tape  / 
\  I TP 4  / 


Set  IFL,  NZYKB 
I  FIRST  and  HAST 
for  body  KB 


/Done\^  MO 
with  body^— - 
\  KB/ 


Begin  00  200 
on  1=11,12 


NZYSV=NZYKB 

KB=KB+1 


Set  new  IFL, 
NZYKB  and 
IZ=0 


hZYSV=2 

and 

sNZYKB=3> 


outside  of^ 
J FIRST, I LAS! 


NO  /  ^ 

(F>KJT2 


LBO  =  KB-1 
LSO  =  NSTRIP+LBO 
OBO  =  I 


Write  I-th 
\ro w  of  DPZy 
\on  tape  / 
\  ITP1  / 


Set  new 
I  FIRST,  I  LAST*, 
KT=KT+1 
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101 


2 


NYFLA6 
=  0 


Begin  DO  800 
on  1*11,12 
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5.3.2  SUBROUTINE  DPPS  (K,  KS,  I,  01,  J2,  SGR,  CGR.  RE FC »  FMACH,  YS,  ZS, 
NBARAY,  NCARAY ,  DT,  WORK) 


Functional  Description 

Subroutine  DPPS  prepares  the  variables  necessary  for  the  computation 
of  one  row  of  the  DPP-submatrix  and  calls  subroutine  SUBP  in  a  do-loop  for 
each  element  of  this  row.  The  resulting  matrix  row,  DT(500),  is  returned  to 
subroutine  GEND  via  the  argument  list  of  subroutine  DPPS. 


Input-Output  Variables 


MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

K 

Panel  number  in  which  the 

KS 

IN 

receiving  point  'i 1  lies 

Strip  number  in  which  the 
receiving  point  1 i *  lies 

I 

i 

IN/QUT 

Receiving  point  index 

01 

Do-loop  delimiters  for  the 

J2 

IN 

number  of  elements  in  one 

row  of  the  DPP  matrix 

SGR 

sinyr 

IN/OUT 

ARG 

sine,  and 

CGR 

COSYr 

cosine  of  the  dihedral 
angle  of  receiving  strip 

REFC 

c 

Reference  chord 

FMACH 

M 

Mach  Number 

YS(2l 0) 

y 

IN 

y-array 

ZS(210) 

z 

See  Sec.  5.2.1 

z-array 

NBARAY (10 

NCARAY (10 

mM 

See  Blank  Common,  Sec.  3.1 

DT ( 500) 
WORK 

IN/OUT 

One  row  of  the  DPP  matrix 
.Complex,  work  array _ 
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MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

10 

NBXS 

See  INCRO;  Sec.  5.1 .7 

NCPNB 

IR 

Index  of  sending  point 

YREC 

y 

y-coordinate' 

of  receiving 

ZREC 

z 

OUT 

DPPS 

z-coordinate 

point  i 

L 

Panel  number" 

in  which  the 

"  sending  point 

LS 

Strip  number J 

'j*  lies 

Calling  Subroutine  GENO 
Called  Subroutine  SUBP 


5.3.3  SUBROUTINE  DPZY  (KB,  KT,  IZ,  I,  Jl,  J2,  IFIRST,  HAST,  REFC, 

FMACH ,  YB,  ZB,  AVR,  ARB,  TH1A,  TH2A,  NT12,  NBARAY, 

NCARAY ,  NZYKB,  DPZ,  OPY,  WORK) 

Functional  Description 

Subroutine  DPZY  prepares  the  variables  necessary  for  the  computation 
of  one  row  of  either  the  DPZ,  or  the  DPY  submatrix  and  calls  subroutine  SUBP 
in  a  double  do-loop  for  each  element  of  the  row,  to  perform  the  summation 
given  in  Eqs  (  5. 3. 3-1  and  -9).  The  resulting  matrix  row  is  returned  to 
the  calling  subroutine  GEND  via  the  argument  list  of  DPZY. 

Input-Output  Variables 


MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

KB 

Body  number  in  which  the 

IZ 

IN/OUT 

ARG 

receiving  point  ’ i ‘  lies 

Body  element  number  of  body 

KB  in  which  'i 1  lies 

I 

i 

Receiving  point  index 
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MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

DELTH 

a  e 

(5. 3.3-3) 

YREC 

yr 

S6e  Et>'S  (5.3.3-A) 

ZREC 

zr 

(5. 3. 3-5) 

RHO 

s 

(5. 3. 3-8) 

SGR 

sinYrj, 

(5. 3. 3-7) 

CGR 

cosY^ 

(5. 3. 3-6) 

SMULT 

OUT 

DPZY 

(5. 3. 3-2) 

CMULT 

(5.3.3-10) 

L 

Panel  number 

in  which  the 

LS 

Strip  number 

sending  point 
'j'  lies 

10 

NBXS 

See  1NCR0;  Sec.  5.1.8 

NCPNB 

IR 

Index  of  sending  point 

Calling  Subroutine  GEND 
Called  Subroutine  SUBP 


Equations 

When  the  receiving  point  ' i '  lies  on  a  z-oriented  body,  subroutine  DPZY 
computes 


Ne 


(b), 


DPZ. •  =  Z  DP ( ARGn ,ARGd ,  ARG  ) S ( 6  ,AR'U/)A6 

10  p=1  '  o*  R*  s  v  w 


(5. 3. 3-1) 


where 

S(e  ,AR(b^)  =  I  sine  J  1  +  cos2e  (AR2-!)  (5. 3. 3-2) 

y  n  U  "  y  ' 
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(5. 3. 3-3) 


46 


6 

_  U 


+  1 


-  1 


where 


°o  =  6Ne  -  2’  and  eNetl  "  9i  +  2* 
and  u  =  1,2,  ...  Ne. 

Note  that  Ne  is  either  Ne^  or  Ne^,  depending  on  the  body  section  in 
which  the  body  receiving  point  lies;  accordingly,  ey  is  either  el^  or 
02  -  see  input  data,  Sec.  1.2,  and  DP(ARGQ,  ARGR,  ARGR)  is  given  by 

Eq^ (5 .5 .3.7-2)  computed  in  Subroutine  SUBP. 

Here,  however,  the  receiving  point  arguments  depend  on  the  e-distribution 
of  the  body  section,  and  are  defined  as  follows: 

ARG  : 

K 


*1  *  XIi 

yk  =  YB^  +  a^  cose^ 

(5. 3. 3-4) 

Zr,  -  ZB^  +  a^  AR^  sine 
k  p 

(5. 3. 3-5) 

cosY  =  sine  /c 

y»  MM 

£ 

(5. 3. 3-6) 

sinyr^  =  -AR^  cose^ 

(5. 3. 3-7) 

5  =  cos2e  (AR2-!)  +  1 
p  p 

(5. 3. 3-8) 

A  complete  summary  of  the  ARGR  and  ARG^  are  given  in  the  table  at  the  end  of 
this  section,  ARGQ  denotes  the  constant  arguments  and  M  -  this  nolds 
true  for  all  nine  submatrices  of  the  total  downwash  factor  matrix  [DT] . 


For  receiving  points  on  y-oriented  bodies  subroutine  DPZY  computes 

or:  \rl  o  ad^Nao  (5.3.3 


Ne 

DPY. j  =  E  ^DP(ARGQ,ARG| 

wher 

C(e  ,AR(b))  =  1  cose 

P  7T  P 

and 

40  is  given  by  Eq.(  ! 

(5.3.3-10) 
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DP(ARGq ,  ARGR,  ARGg)  is  given  by  Eq.  (5. 3. 7-2)  and  the  arguments  ARGR  and 
ARGg  are  summarized  in  the  table  below. 


Receiving  Point  on  Body,  Sending  Point  on  Panel 


POINTS 

ARGUMENTS 

COORDINATES 

e 

ae 

AR 

n 

■a 

z 

Receiving 
point  ' i * 

ARGr 

p 

xli 

yp 

See  Eqs 

throi 

z 

p 

>.(5. 3. 3-4 
igh  -7) 

Yrk 

e 

p 

( 

ae 

p 

See  Eq. 

5.3.3 

-3) 

Sending 
point  * j 1 
strip  V 

ARGS 

S(ARGS) 

G(ARGS) 

S(G(ARG.S): 

\ 

\ 

\ 

\ 

Egg 

Box 

semi- 

width 

Aver¬ 

age 

box 

chord 

S 

S 

e£ 

C£ 

A  detailed  description  of  the  "Arguments'  above  is  given  in  Sec.  5.3.7  along 
with  their  usage  in  the  computation  of  the  DP-matrix  element  components. 

5.3.4  SUBROUTINE  DYPZ  (KB,  KS,  LS,  IZ,  1 ,  01 ,  J2,  NYFLAG,  FLND,  FLNE,  SGR, 
CGR,  REFC,  FMACH,  KR,  ARB,  NBEA,  LBO,  LSO,  JBO,  DT) 

Functional  Description 

Subroutine  DYPZ  prepares  the  variables  necessary  for  the  computation 
of  one  row  of  either  the  DYP-,  or  the  DYZ-,  or  the  DYY  submatrix,  depending 
on  the  location  of  the  receiving  point.  In  either  case  it  calls  subroutine 
SUBB  in  a  do  loop  for  each  element  of  a  row;  latter  is  returned  to  the 
calling  subroutine  GEND  via  the  argument  list  of  DYPZ. 
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Input  Output  Variables 


MNEMONIC  SYMBOL  IN/OUT  SOURCE 


DESCRIPTION 


Body  number  in  which  receiving 
point  ' i '  lies 

Index  of  receiving  point  y- 
and  z-coordinates 


See  Subroutine  DPZY 
Sec.  5.3.3 


Do-loop  delimiters  for  the 
number  of  elements  in 
one  row  of  the  submatrix 

0  for  DYP  and  DYZ  elements, 
1  for  DYY  elements 


See  Blank  Common 
Sec.  3.1 


Sequence  number  of  first  body 
with  y-orientation 
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MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

LSO 

IN 

ARG 

y-  and  z-coordinate  index  for 
first  y-oriented  body  element 

JBO 

Sending  point  index  for  first 
y-oriented  body  element 

DT ( 500) 

IN/OUT 

ARG 

One  row  of  either  of  the 

submatrices  DYP,  DYZ  or  DYY 

NDY 

Flag  used  in  subroutine  SUBB  - 
1  for  sending  points  in 
y-oriented  bodies 

LB 

Body  number  in  which  sending 
point  ' j '  is  located 

JB 

OUT 

DYPZ 

Sending  point  index  (column 
no.  of  DT-matrix) 

JZ 

Body  element  number  for 
sending  body  LB 

NZYLB 

IN 

z-y  orientation  flag  for 
body  LB 

SL 

CL 

TL 

SGS 

CGS 

OUT 

sinA=  0  Variables  used  in 

cosa=  1  the  called  subrout- 

tanx=  0  ine  SUBB:  x=0  and 

sinys=  -1  ys  "  -90°  for  send- 

cosys=  0  ing  elements  on  y- 

oriented  bodies 

Calling  Subroutine  GEND 

Called  subroutine  SUBB 


5.3.5  SUBROUTINE  DZPY  (KB,  KS,  LS,  IZ,  I,  J1 ,  J2,  NY FLAG,  FLND,  FLNE, 
SGR,  CGR,  REFC,  FMACH,  KR,  ARB,  NBEA,  DT) 


Functional  Description 

Subroutine  DZPY  prepares  the  variables  necessary  for  the  computation 


114 


I 


of  one  row  of  either  the  DZP-,  or  the  DZZ-,  or  the  DZY  matrix,  depending  on 
the  location  of  the  receiving  point.  In  either  case  it  calls  subroutine 
SUBB  in  a  do  loop  for  each  element  of  a  row;  latter  is  returned  to  the  calling 
subroutine  GEND  via  the  argument  list  of  DZPY. 

Input  Output  Variables 


MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

KB 

KS 

LS 

IZ 

I 

J1 

J2 

NYFLAG 

FLND 

FLNE 

SGR 

CGR 

REFC 

FMACH 

KR 

ARB (10) 

NBEApO,2) 

i 

6 

e 

sinyr 

cosy., 

c 

M 

kr 

IN 

ARG 

See  Subroutine  DYPZ 

Sec.  5.3.4 

DT(500) 

IN/OUT 

One  row  of  either  of  the  sub¬ 
matrices  DZP,  DZZ  or  DZY 

NDY 

OUT 

DZPY 

Flag  used  in  subroutine  SUBB; 

0  for  sending  points  in  z- 
oriented  bodies 

MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

LB 

OUT 

JB 

See  Subroutine  DYPZ 

JZ 

Sec.  5.3.4 

NZYLB 

IN 

DZPY 

SL 

sinx=0  Variables  used  in 

CL 

cosx=l  subroutine  SUBB; 

A=0,  and  ys=0  for 

TL 

taiU-°  sending  elements 

SGS 

OUT 

sinys=0  on  z-oriented 

CGS 

C0SV1  bod1es 

Calling  Subroutine  GEND 

Called  Subroutine  SUBB 


5.3.6  SUBROUTINE  SUBB  (KB,  KS,  I,  J,  OZ;  JB,  LB,  LS,  NDY,  NYFL,  FLND, 

FLNE,  PI,  EPS,  SGR,  CGR,  SGS,  CGS,  AR,  SL,  CL,  TL,  FL,  BETA,  SUM) 


Functional  Description 

Subroutine  SUBB  computes  the  downwash  factor  matrix  elements  for 
all  receiving  points  and  sending  points  on  interference  bodies,  one  element 
at  a  time,  according  to  Eqs.  (5. 3. 6-1  and  -2) or  Eq.  (5. 3. 6-3  and  -4) 
depending  on  the  orientation  of  the  interference  body  in  which  the  sending 
point  is  located.  The  actual  computation  of  the  components  to  the  matrix 
elements  is  done  in  Subroutine  DZY  which  is  called  from  SUBB.  The  final 
result  -  one  downwash  factor,  SUM  -  is  returned  to  the  calling  subroutine 
GEND  via  the  argument  list  of  SUBB. 


Input  Output  Variables 


MNEMONIC 

i 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

KB 

IN 

ARG 

Index  of  receiving  body  -  0 
when  receiving  point  is  on  panel 

116 


MNEMONIC  SYMBOL  IN/OUT  SOURCE 


DESCRIPTION 


'Strip'  number  in  which 
receiving  point  lies 

Receiving  point  index  - 
row  number  of  DT-matrix 

Sending  point  index  - 
column  number  of  DT-matrix 


NYFL, 

NYFLAG 


See  Subroutine  DYPZ 
Sec.  5.3.4 


0.00001 


sinx  =  0 
cosx  =  1 
tanx  =  0 
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MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

BETA 

0 

IN 

ARG 

0  =  V  1  -  M2 

SUM 

°Tio 

IN/OUT 

One  element  of  DT-matrix  where 
the  sending  point  * j *  is  an 
interference  body  element 

IND 

BR 

ANOT 

DXS 

c/2 

OUT 

SUBB 

Flag  for  subroutine  TKER;  0 
for  'total'  kernel,  1  otherwise 

Reference  semi -chord 

Local  characteristic  half¬ 
width  of  interference  body 
element  ' j ' 

Body  element  length 

TEST1* 

TEST2* 

IN 

|ycr  - 

r  s 
lzc  -  2cl 

same  for  se 

line 

where  y*^  ,  zT  are 
the  y-  and  z  coord¬ 
inates  of  receiving 
body  centerline  - 

and  y*  zs  are 
c  c 

nding  body  center- 

D2D 

See  Eq.  (5. 3. 6-2) 

XX 

Y 

Z 

XII 

XI2 

ETA 

ZETA 

AO 

IDZDY 

OUT 

See  argument  list  of 

Subroutine  DZY; 

Sec.  5.1 .3 

*  Used  only  if  receiving  point  is  an  interference  body  element 
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MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

IGO 

IN 

SUBB 

Internal  flag,  1  through  4, 
corresponding  to  the  'quadrant' 
in  which  the  sending  point 

lies: 

1  -  upper  right  (UR) 

2  -  upper  left  (UL) 

3  -  lower  right  (LR) 

4  -  lower  left  (LL) 

DPUR 

UR 

DPUL 

UL 

contribution  to  the 

DPLR 

LR 

downwash  factor 

DPLL 

LL 

Calling  Subroutines  DZPY  and  DYPZ 


Called  Subroutine  and  Common  Blocks  DZY,  the  Blank  Coinnon  Block  and  the 
Labeled  Common  Block  KDS 


Equations 


A.  Sending  points  * j  *  on  z-oriented  bodies  (6). 

When  the  receiving  points  'i1  are  on  panels*  subroutine  SUBB  computes 


DZPi j  =  DZUR  +  6DZUL  +  eDZLR  +  6eDZLL 


(5. 3. 6-1) 


where  the  DZ(ARGQ,  ARGR,  ARGj)  elements  are  described  in  Sec.  5.1.3 
(Subroutine  DZY). 


When  the  receiving  points  i  are  on  a  z-oriented  body  *b ’  subroutine 
SUBB  computes  DZZ,.j  as  follows: 

a)  If  the  two  bodies  b  and  b  are  not  identical  (b  f  5)  and  their 
centerlines  do  not  coincide, 

DZZij  *  DZUR  +  SDZUl  +  eDZLR  +  6eDZLL 
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b)  If  b  f  b,  but  the  centerlines  of  the  two  bodies  coincide 

DZZu  =  0 

c)  if  b  =  b,  but  i  t  j 

DZZio  -  0 

d)  if  b  =  b  and  i  =  j 

DZZij  =  D2D 


(5, 3. 6-2) 


2u(a0b)2  )1  +  ARb) 


When  the  receiving  point  * i *  is  on  a  y-oriented  body,  subroutine  SUBB 
computes  DZY- •  as  follows: 

1  J 

a)  if  b  f  b 

DZYi j  =  dzur  +  6DZUL  +  eDZLR  +  6eDZLL 

b)  if  b  =  b 

DZYij  "  0 

B.  Sending  points  'j'  on  y-oriented  bodies  (b). 

When  the  receiving  points  ' i '  are  on  panels,  subroutine  SUBB  computes 


DYPij  =  DYUR  +  6DYUL  +  cDVLR  +  6gDLL 


(5. 3. 6-3) 


where  the  DY(ARGq,  ARGr,  ARGs)  elements  are  described  in  Sec.  5.1.3. 

When  the  receiving  points  are  on  z-oriented  bodies 

a)  if  b  f  b 

°YZi  .  =  0Yur  +  6DYul  +  cOYlr  +  5eDYLL 

b)  if  b  =  b,  or  b  f  b  but  their  centerlines  coincide 

DYZlj  -  0 
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When  the  receiving  points  are  also  on  y-oriented  bodies 

a)  if  b  t  b  and  their  centerlines  do  not  coincide 

OYYi  •  =  DYur  +  6DYul  +  eDYQR  +  6cDYLL 

b)  if  b  f  b,  but  their  centerlines  coincide 

DYY .  •  =  0 

*  J 

c)  if  b  =  b,  but  i  f  j 

DYYl .  =  0 


d)  if  b  =  b,  and  i  =  j 
DYY-.  =  DY2D  =  D2D/AR 

•  J 

_  1  1 
"  m  ' 


2*(ab)2  (l+ARb) 


(5. 3.6-4) 


The  arguments  that  are  used  in  the  computations  of  the  DZ  and  DY 
elements  are  tabulated  below. 


Receiving  Point  on  Panel,  Sending  Point  on  Body 


POINTS 

ARGUMENTS 

COOF 

IDINATES 

DIHEDRAL 

ANGLE 

HALF 

WIDTH 

OF  BODY 

AVG  j 

RflnTIK  ,  AR 

Hi 

m 

z 

Receiving 
point  i, 
strip  k 

ARGr 

xi 

m 

B 

YRk 

Sending 
point  j 

ARGS 

S( ARG  ) 

S 

G(ARG$) 

S(G(ARGs) 

% 

and 

% 

\ 

\ 

“nc 

4 

-C 

Ys=0  for 
z-body 

Ysc  -tt/2 

for  y- 
body 

* 

e 

ao 

AR(b) 

*Note  that  e=aQ/2»/  | l-ARl2’  for  z-oriented,  AR>1  and  for  y-oritT.ed  AR<1  bodies 


and  e  =  aQ/372/  ) l -AR | 2  for  z-oriented,  AR<1  and  for  y-oriented,  AR>1  bodies. 
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Receiving  Point  on  Body,  Sending  Point  on  Body 


POINTS 

COORDINATES 

DIHEDRAL 

ANGLE 

ax 

m 

AR 

ARGUMENTS 

■ 

■ 

z 

OF  BODY  SECTION 

Receiving 
point  i 

ARGr 

XI  i 

-Q 

CO 

>- 

V’ 

Yr=0  for 
z-body, 

Yr  =  -n/2 

for  y- 
body 

«k 

ao 

1 

Sending 
point  j 

ARGS 

See  Table  above 

Flow  Chart  -  Subroutine  SUBB 
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5.3.7  SUBROUTINE  SUBP  (I.  L,  LS.  J.  IQ,  JR.  NBXS,  NCPNB.  SGR.  CGR.  YREC, 

ZREC,  SUM,  WORK) 


Functional  Description 

Subroutine  SUBP  computes  the  downwash  factor  matrix  elements,  DP-  -  for 

*  J 

all  receiving  points  on  panels  and  interference  bodies  and  sending  points  on 
panels,  one  element  at  a  time,  according  to  Equations  (5.3. T-1!  through  -6). 

The  computation  of  the  individual  components  that  make  up  each  downwash  factor 
element  is  done  in  Subroutines  SNPDF  and  INCRO,  which  are  called  from  SUBP. 
Additional  arguments  that  are  needed  for  the  calculation  of  the  downwash 
factor  contribution  of  image  sending  points  inside  associated  bodies  are  cal¬ 
culated  by  subroutine  SUBI,  which  is  also  called  from  SUBP.  The  resulting 
total  downwash  factor  element  'SUM'  is  returned  to  the  calling  subroutine  via 
the  argument  list  of  SUBP. 


Input  Output  Variables 


MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

I 

L 

See  DPPS,  Sec.  5.3.2 

LS 

J 

j 

Running  index  of  the  element 
number  in  DPP  matrix  row  (column 
number) 

10 

IR 

See  INCRO,  Sec.  5.1.7 

NBXS 

IN 

ARG 

NCPNB 

SGR 

sinyr 

CGR 

cosyr 

YREC 

yr 

See  DPPS,  Sec.  5.3.2 

ZREC 

zr 
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MNEMONIC  SYMBOL  IN/OUT 


SOURCE 


DESCRIPTION 


SUM 

IN/OUT 

! 

One  element  of  the  DT  matrix  for  receiving 
point  1 i 1 ,  sending  point  1 j 1 »  when  ' i 1  is 
on  a  panel.  Note  that,  if  * i '  represents 
a  receiving  interference  body  section, 

SUM  is  only  the  contribution  (to  this 
DT-element)  of  one  point  on  the  surface 

of  the  body  section  - 

-  see  Eq . ( 5 . 3 . 3-1 ) . 

TL 

SL 

CL 

XO 

YO 

I 

See  Subroutine  SNPDF, 

Sec. 5.1 .1 1 

ZO 

ES 

CV 

■ 

OUT 

SUBP 

AZ 

AX1 

AY1 

. 

See  Subroutine  INCRO, 

,  Sec.  5.1.8 

AZ1 

AZ2 

DIJS 

D 

IN 

Argument 

Steady  contribution 

of  sending  point  ‘j 

List  of 

to  the  downwash 

on  panel 

of  the  image  of  'j' 

JJIJI 

I  CD) 

SNPDF 

factor  1 

• 

inside  the  current 
associated  body 

I 


WORK 


IN/OUT 


ARG 


Complex  work  array 


MNEMONIC  SYMBOL  IN/OUT  SOURCE 


DELRS  ReUD) 


DELIS  Re  C I ( AD ) ] 

DELRI  Im(AD) 
DELII  Iiti[I(aD)] 


INCRO 


NA1 

NA2 


DA 

DZB 

DYB 

DAR 

DETA 

DZETA 

DCGAM 

DSGAM 

DEE 


DXI 

DETAI 

DZETAI 

DCGAMI 

DSGAM I 

DEEI 

DTLAMI 


DESCRIPTION 

Real  part  of  the 

• 

of  sending  point  'j' 

unsteady  contribution  < 

on  panel 

to  the  downwash 

of  the  image  of  'j' 

factor 

L 

Imaginary  part  of 

unsteady  contribution 

to  the  downwash 

factor 

See  Eqs.  (5. 3. 7-2  through  -6) 

The  number  of  associated  bodies  for  panel 

in  which  the  sending  point  'j1  lies 

Delimiters  to  the  do-loop  for  all 

associated  bodies  of  panel 

Arguments  to  Subroutine  SUBI;  see 

Sec.  5.1.12 

Image  point  arguments; 

see  Sec.  5.1 .12 
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I(ARG<«)*  =  images  of  ARG^  inside  associated  body  IR 

S(ARGg)  =  images  of  ARG<.  with  respect  to  the  y  =  0  plane 

G(ARGg )  =  images  of  ARG^  with  respect  to  the  z  =  0  plane 

The  arguments  that  constitute  ARGR  and  ARGS  are  tabulated  below. 


Receiving  Point  on  Panel,  Sending  Point  on  Panel 


Points 

Arguments 

Coordinates 

Dihedral 

Angle 

BBS 

IM 

Kr! 

Arg 

Box 

Chord 

B 

n 

z 

Receiving 
point  1 i 1 , 
strip  ' k1 

ARGr 

y . 

1 

a 

Yrk 

Sending 
point  1 j ' , 
strip  V 

ARGS 

nc* 

YS£ 

WSk 

CA 

S(ARGS) 

-nc£ 

B 

m 

G(ARG$) 

nc* 

-?c£ 

S[G(ARGS)] 

~n_ 

ci 

Xsj 

For  receiving  points  on  interference  body  elements,  the  ARGR  are  defined  in 
Section  5. 5. 3. 3  (Subroutine  DPZY). 


The  evaluation  of  D  in  Equations  (5. 3. 7-3  through  -6) is  done  by  the  follow¬ 
ing  subroutines: 

SNPDF  -  for  the  steady  case  (kr  =  0),  yields  D^s 

SNPDF,  INCRO,  TKER,  I0F1  and  IDF2  -  for  the  uisteady  case,  yielding  AD 

(aD  =  0  for  steady  cases). 

Then,  in  general 

D  =  D(s)  +  aD  (5.3. 7-7) 

Description  of  the  above  five  subroutines  is  given  in  Sections.!. 


*I(ARGS)  is  computed  by  Subroutine  SUBI  described  in  Sec.  5.1. 
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DP=CMPLX( (DIJS-DI JI ) - ( DELRS-DELRI ) ,  -(DELIS  -  DELI  I ) ) 
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See  Subroutine  SUBB 
for  notation 


where 

ARGq  represents  the  constant  arguments  kr  and  M; 

ARGp  represents  the  variable  receiving  point  arguments 
ARG<j  =  sending  point  arguments 
and 

/  n\ 

RANGE  refers  to  the  bodies  associated  with  panel  p  in  which  the 
sending  point  lies 
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5.4  Segment  4 


5.4.1  SUBROUTINE  RDMODE(IA,  NA,  NIN,  NM,  NOUT) 

Functional  Description 

This  routine  reads  the  panel  and  body  modal  data  from  cards.  The  data  for 
each  modal  coefficient  consists  of  the  coefficient  and  an  integer  containing 
the  mode  number,  the  panel  or  body  number,  the  exponent  for  the  x  coordinate, 
and  a  flag  indicating  whether  relative  or  absolute  coordinates  for  x  and  y 
are  to  be  used.  This  data  is  then  sorted  according  to  mode  number,  then 
panel  or  body  number  and  then  the  x  and  y  exponents,  respectively.  The 
routine  returns  the  input  coefficients  and  their  identifiers,  the  number  of 
coefficients  for  panels  and  the  z  and  y  oriented  bodies  and  the  number  of 
modes. 

Input-Output  Variables 


MNEMONICS 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

IA 

OUT 

ARG 

Array  containing  the  modal  coefficients 
and  coefficient  description  for  panels 
and  bodies. 

NA 

OUT 

ARG 

Array  containing  the  number  of  co¬ 
efficients  for  panels  and  z  and  y 
oriented  bodies. 

NIN 

IN 

ARG 

I/O  unit  number  containing  the  input 
modal  data. 

OUT 

ARG 

Number  of  modes 

IN 

ARG 

I/O  unit  number  for  printing  on. 

Calling  Subroutines  MAIN 

Error  Messages 

SUBROUTINE  **RDM0DE**  INVALID  DATA  CARD.  —  IGNORED  - 

The  data  card  which  has  been  read  is  invalid  and  has  been  ignored. 
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5.5  Segment  5 

5.5.1  SUBROUTINE  SB  (A,  NM,  NAY,  NAZ,  WORK,  NWORK) 


Functional  Description 

Subroutine  SB,  'Slender  Body',  calculates  the  normalwash  at  lifting  surface 


boxes  and  interference  body  elements  due  to  slender  body  elements. 
Input  Output  Variables 


MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

A 

IN 

ARG 

Array  containing  the  modal  data 

mmm 

read  by  subroutine  RDMODE. 

X 

Array  of  box  and  interference  body 

AR 

AR 

IN 

IBB 

elements  x-coordinate. 

Array  of  body  aspect  ratios. 

AO 

ao 

Block 

Array  of  body  radius. 

KR 

kr 

Reduced  frequency. 

NB 

Number  of  bodies. 

— 

IN 

ARG 

Number  of  modes. 

AOP 

mm 

IN 

Array  of  rate  of  change  of  body  radius 

HU 

Block 

with  rexpect  to  x. 

NAY 

IN 

ARG 

Number  of  modal  coefficients  for 

y-oriented  bodies. 

NAZ 

IN 

ARG 

Number  of  modal  coefficients  for 

z-oriented  bodies. 

■  ■ 

NBY 

NBY 

Number  of  y-oriented  bodies. 

NBZ 

NBZ 

Number  of  z-oriented  bodies. 

NTP 

NTP 

Number  of  lifting  surface  boxes. 

NTY 

NTY 

■ 

Number  of  y-oriented  interference 

NTZ 

NTZ 

IN 

IB 

body  elements. 

Number  of  z-oriented  interference 

Block 

body  elements. 

CBAR 

c 

Reference  chord  length. 

NSBE 

NSBE 

Array  of  number  of  slender  body 
elements  per  body. 

NTYS 

NTYS 

Number  of  y-oriented  slender 
body  elements. 
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NTZS 


FMACH 

NWORK 


SYMBOL 

IN/OUT 

NTZS 

S 

?S1 

IN 

5s2 

M 

SOURCE  DESCRIPTION 

Number  of  z-oriented  slender 
body  elements 
Complex  working  array 
ARG  Array  containing  the  leading  edge 

x-coordinate  of  the  slender  body  elements 
Array  containing  the  trailing  edge 
x-coordinate  of  the  slender  body  elements 
Mach  number 

Length  of  working  array  WORK 


Calling  Subroutines  MAIN 

Called  Subroutines  and  Conmon  Blocks 

Blank  Common,  RWREC,  MYZYC,  DUMULT,  READD,  DZYMAT,  OUTCOR 
Flow  Chart  -  Subroutine  SB 
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5,5.2  SUBROUTINE  DUMULT  (NT,  N2,  NTZS,  NTYS,  W,  DZ,  UZ,  DY,  UY) 
Functional  Description 

This  subroutine  performs  the  following  complex  matrix  operation. 

{“out}  '  [°z]  {>*4  +  [Dy]  { “y }  *  {“w} 


The  result  of  this  operation  (aw)  is  returned  to  the  calling  routine  as  w. 
Input  Output  Variables 


MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

DY 

Dy 

Input  matrix. 

DZ 

Dz 

Input  matrix 

N1 

The  first  row  of  aw  to  calculate 

N2 

The  last  row  of  aw  to  calculate 

NTYS 

IN 

ARG 

Number  of  columns  in  Dy 

NTZS 

Number  of  columns  in  Dz 

UY 

n 

Input  matrix 

UZ 

M 

Input  matrix 

M 

AW 

IN/OUT 

ARG 

Input  matrix  and  output  result 

Calling  Subroutines  SB,  OUTCOR 

5.5.3  SUBROUTINE  DZYMAT  (D,  NFB,  NLB,  NTZYS,  IDZDY,  NTAPE,  XP,  BETA) 
Functional  Description 

This  subroutine  sets  up  the  proper  argument  lists  for  the  calculation  of 
each  row  of  the  Dz  or  matrix  and  then  calls  subroutine  ROWDYZ. 

Input  Output  Variables 


MNEMONIC 

SYMBOL 

IN/OUT 

-  -  -  .  - 

SOURCE 

DESCRIPTION 

D 

Working  array  used  to  store  a  row 

IN 

ARG 

of  Dz  or  Dy 

CG 

COSy 

Array  containing  the  cosine  of  the 

lifting  surface  strip  dihedral  angle 
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MNEMONIC 


SYMBOL 


IN/OUT 


SOURCE 


DESCRIPTION 


NB 

mm 

i 

Total  number  of  bodies 

NC 

Array  containing  the  number  of  chord- 

wise  boxes  per  strip  in  a  panel 

NP 

NP 

Number  of  panels 

NS 

NS 

Array  containing  the  number  of 

strips  in  a  panel 

SG 

siny 

Array  containing  the  sine  of  the 

lifting  surface  strip  dihedral  angle 

XP 

IN 

ARG 

x-control  point  coordinate  of  lifting 

surface  boxes 

YB 

y-coordinate  of  center  of  bodies 

YP 

y-control  point  coordinate  of  lifting 

surface  strip 

ZB 

z-coordinate  of  center  of  bodies 

ZP 

z-control  point  coordinate  of 

lifting  surface  strip 

NBY 

Number  of  y  oriented  bodies 

NBZ 

Number  of  z  oriented  bodies 

NFB 

Number  of  the  first  body  with  the 

orientation  requested 

NLB 

Number  of  the  last  body  with  the 

orientation  requested 

NTP 

Total  number  of  boxes 

BETA 

6 

IN 

ARG 

V 1  -  M2 

MACH 

M 

Mach  number 

NBEA 

Array  containing  number  of  interferen 

body  elements  per  body  and  the  body 

orientation 

IDZDY 

Flag  indicating  whether  the  Dz(o) 

or  D  (1)  matrix  is  to  be  calculated 

NTAPE 

I/O  unit  number  which  the  output 

matrix  is  to  be  written  on 

NTZYS 

Number  of  z  or  y  oriented  slender 

body  elements 

Calling  Subroutines 

SB 

— 

- 1 - - — 

Called  Subroutines  and  Common  Blocks 
ROlJDYZ,  Blank  Common  Block 
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5.5.4  SUBROUTINE  MUZYC  (NMODE,  NCOEF.  K,  NTZY ,  NF60DY,  NLBODY ,  NSBE,  KR, 
IA,  A,  CBAR,  AO,  AOP,  XIS1,  XIS2,  AR,  UZY,  CPZYl 

Functional  Description 

Subroutine  MUZYC  calculates  the  axial  doublet  strengths  and  loading  for 
slender  bodies. 

Input  Output  Variables 


MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

A 

Real  array  containing  the  modal 
coefficients  and  a  key  identifying 
the  coefficients. 

K 

Code  identifying  which  type  coefficient 
to  use 

=  2  use  az 
=  3  use  a 

AR 

AR 

IN 

ARG 

y 

Array  containing  aspect  ratios  of 

the  bodies. 

AO 

ao 

Array  containing  the  radii  of  the 
bodies. 

IA 

Integer  array  equivalent  to  the 
array  "A". 

KR 

kr 

Reduced  frequency. 

AOP 

ao 

UZY 

V  y 

OUT 

The  ^  or- Py  column 

CBAR 

c 

IN 

Reference  chord  length. 

CPZY  AC  aA(1+AR)  OUT 

ARG 

A 

C  aA(1+A 

py 

R)  OUT 

NSBE 

IN 

Array  containing  the  number  of  slender 
body  elements  per  body 

NTZY 

NTZS 

IN 

ARG 

Total  number  of  z  or  y  oriented  slender 

NTYS 

body  elements 

XIS1 

?6l 

IN 

Array  containing  the  leading 
edge  x-coordinate  of  the  slender 
body  elements. 

iff 


MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

XIS2 

?S2 

Array  containing  the  trailing  edge 
x-coordinate  of  the  slender  body  elements 

NCOEF 

Number  of  modal  coefficients  inpjt  for 
this  orientation  of  body. 

nni 

IN 

ARG 

Number  of  modes  in  analysis. 

Number  of  the  first  body  with 
orientation  requested. 

NLBODY 

Number  of  the  last  body  with 
orientation  requested. 

Calling  Subroutines 


SB 


Equations 

_  — *  WS 

uz  and  ^  =  D2D 

AC  aA  AR  \ 

Z 

and  |  =  [*(ws)L  +  i (ws " )  (M)j  2a0  AR 
AC  aA  AR  1 

py  I 

where  for 

K  =  2  or  z  subscript 

D2D  =  2-rr  (1  +’aR)  a0z 

L  =  2n  (a0'  +  i  kr/c} 

M  =  u  ao/c" 
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K  =  3  or  y  subscript 


L  =  2tt  /ft  (a,/  +  i  kr/  c) 

M  =  &  it  aQ/  c 


and 

ws  = 

:  wsR  +  i  wsl 

wsR 

=  E  A  n  ( 

n=o 

9 

/  \  n 

wsl 

■  T,  A  2kr 

n=l 

w 

ws  " 

=  wsR "  + 

i  wsR(2kr) 

9 

/yc\  n_2 

wsR'”* 

=  ^  A  n(n 
n=l 

-’)  (f) 

5.5.5  SUBROUTINE  OUTCOR  (WORK,  NWORK,  NTS,  N,  NM,  NUTAP,  NDWIN, 
NDWOUT ,  NDTAP) 

Functional  Description 


This  routine  performs  the  following  matrix  operation 


[aw]  =  [aw  ]  +  [D]  [U] 

It  is  assumed  that  all  the  matrices  will  not  fit  into  core.  As  many  columns 
of  U  and  aw  '  are  read  into  core  as  possible  and  the  matrix  D  is  read  a 
row  at  a  time  performing  the  necessary  operations  to  calculate  the  columns 
of  aw  which  are  written  on  an  I/O  unit.  This  operation  is  repeated  until 
all  the  columns  of  U  have  been  read  and  hence  all  of  aw  written. 
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Input  Output  Variables 


MNEMONICS 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

BE*?® 

Complex  working  array. 

BSE&SB 

Length  of  working  array. 

NTS 

Number  of  columns  in  the  matrix  D. 

N 

Number  of  rows  in  the  matrix  D. 

NM 

Number  of  columns  in  the  matrix  U. 

NUTAP 

I/O  unit  number  containing  the 

matrix  U  in  column  sort. 

NDWIN 

IN 

ARG 

I/O  unit  number  containing  the  aw' 
matrix  in  column  sort.  If  it  is  equal 

0,  the  w  matrix  is  set  to  zero. 

NDWOUT 

I/O  unit  number  which  the  Aw  matrix 

is  written  on  in  column  sort. 

NDTAP 

I/O  unit  number  containing  the  D 

matrix  in  row  sort. 

Calling  Subroutines  SB 

Called  Subroutines  and  Coirmon  Blocks 
RWREC ,  DUMULT,  READD 
Error  Messages 

SUBROUTINE  **OUTCOR**  NUMBER  OF  MODES  IS  LESS  THAN  OR  EQUAL  TO  ZERO. 
***CALCULATIONS  SKIPPED**'4' 

5.5.6  SUBROUTINE  ROWDYZ  (NFB,  NLB,  ROW,  NTZYS,  D,  DX,  DY,  DZ,  BETA, 
IDZDY,  NTAPE,  SGR,  CGR) 

Functional  Description 


This  routine  performs  the  logic  required  to  set  up  the  argument  list  to  DZY 

for  the  purpose  of  calculating  a  row  of  the  Dz  or  D  matrix. 

Input  Output  Variables 


MNEMONICS 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

D 

OUT 

ARG 

The  output  row  of  D_  or  D..  (complex) 

mm 

m 

IN 

z  y 

Aspect  ratio  at  the  body. 

HU 

ao 

IN 

Block 

Radius  of  the  body. 

DX 

X 

IN 

ARG 

x-coordinate  of  receiving  point. 
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MNEMONICS 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

DY 

Y 

IN 

y-coordinate  of  receiving  point. 

DZ 

Z 

IN 

Bi 

z-coordinate  of  receiving  point. 

KR 

— 

mm 

Reduced  frequency. 

NO 

n 

Symmetry  flag. 

NE 

IN 

Ground  effects  flag. 

YB 

Block 

Array  containing  y-coordi nates  of  bodies 

ZB 

Array  containing  z-coordi nates  of  bodies. 

CGR 

GOSYr 

Cosine  of  dihedral  angle  of  receiving 
point. 

NFB 

Number  of  the  first  body  having  the 
desired  z  or  y  orientation. 

NLB 

IN 

ARG 

Number  of  the  last  body  having  the 
desired  z  or  y  orientation. 

ROW 

The  row  number  of  Dz  or  to  be 

calculated. 

SGR 

Sinyr 

Sine  of  dihedral  angle  of  receiving 
point. 

BETA 

3 

•y/l  -  M2 

CBAR 

c 

mm 

Length  of  reference  chord 

MACH 

M 

Mach  number. 

NSBE 

EBB 

Array  containing  the  number  of  slender 
body  elements  per  body 

XIS1 

?S1 

IN 

Array  containing  the  x-coordinate  of 
the  leading  edges  of  the  slender 
body  elements 

XIS2 

5S2 

Array  containing  the  x-coordinate  of 
the  training  edges  of  the  slender 
body  elements 

Calling  Subroutines  DZYMAT 


Called  Subroutines  and  Common  Blocks  Subroutine  DZY  and  the  Blank  Common 
Block 
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See  Subroutine  SUBB 
for  notation 


5.6  Segment  6 


5.6.1  SUBROUTINE  WANDWT  (A,  IA,  NSARRY,  NBARAY>  X,  YP,  Y1 ,  ZP,  Zl , 

NR,  CBAR,  KR,  NP,  NTP,  NM,  COEF,  NDW,  NWT,  NOUT,  W,  DM,  NB,  I  PRINT) 

Functional  Description 

This  subroutine  calculates  the  complex  boundary  conditions  (w)  on  the 
lifting  surfaces  due  to  lifting  surface  motion  and  adds  to  it  the  incremental 
normalwash,  aw,  due  to  the  slender  bodies.  Also  included  in  aw  is  a  normalwash 
induced  at  the  interference  body  elements. 

Input  Output  Variables 


MNEMONICS 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

A 

a 

IN 

ARG 

Array  containing  modal  coefficients 

W 

w 

OUT 

ARG 

A  column  of  the  matrix  w 

X 

xp 

IN 

ARG 

3/4  chord  x-coordinate  t  lifting 

surface  box 

DW 

AW 

IN 

irzilM 

Input  aw  matrix  column 

WT 

OUT 

i/o(nwt; 

The  output  matrix  column  WT 

IA 

Coded  array  describing  the  modal 

coefficients 

KR 

kr 

Reduced  frequency 

NM 

NM 

Number  of  modes 

NP 

NP 

Number  of  lifting  surface  panels 

NR 

N 

Total  number  rf  rows  in  the  WT  matrix 

YP 

yps 

IN 

ARG 

y-coordinate  of  the  lifting  surface 
strip  3/4-chord  point 

Y1 

Y1 

y-coordinate  of  the  inboard  edge  of  panel 

ZP 

ZPS 

z-coordinate  of  3/4-chord  point  of 
the  lifting  surface  strip 

Zl 

Zl 

z-coordinate  of  the  inboard  edge  of  panel 

NDW 

I/O  unit  number  containing  the  aw  matrix 

NTP 

NTP 

Number  of  lifting  surface  boxes 

NWT 

I/O  unit  number  on  which  the  WT 

matrix  is  written 

CBAR 

c 

Length  of  the  reference  chord 
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MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

_ _ 1 

NBARAY 

IN 

ARG 

NSARRY 

NS 

NB 

IPRINT 

DESCRIPTION 

I/O  unit  number  for  printing  output  on. 
Number  of  modal  coefficients  in  the 
array  A. 

Array  containing  the  number  of  the 
last  box  in  each  panel. 

Array  containing  the  number  of  strips 
in  a  panel . 

Number  of  bodies. 

Pr.nt  flag. 


Calling  Subroutines 


MAIN 


The  array  describing  the  modal  coefficients  is  an  integer  at  least  7  to  8 


the  mode  number 
the  panel  number 
m  in  the  equation 

n  in  the  equation 

[>18  in  the  ^  equation 
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5.7  Segment  7 

5.7.1  SUBROUTINE  SOLVIT  (A,  RA,  ND,  MD,  KD,  NI,  NM,  NO,  NW,  NPR1 ) 


Functional  Description 

Subroutine  SOLVIT  solves  the  system  of  simultaneous  linear  equations  repre¬ 
sented  by  the  augmented  rectangular  (n  x  m)*  matrix  [DT'jWT],  which  is 
written  on  logical  tape  unit  NI,  row  by  row,  in  the  MAIN  program.  All  other 
information  necessary  for  the  operation  of  SOLVIT,  is  entered  via  its  argument 
list.  The  solutions  obtained  by  SOLVIT  are  saved  on  logical  tape  unit  NW  in 
column  order,  i.e.,  one  set  of  solutions  per  record;  the  data,  entered  from 
input  tape  NI  is  not  preserved. 


A  detailed  description  of  Subroutine  SOLVIT  can  be  found  in  Reference  3. 
Here  we  give  only  a  brief  description  of  the  variables  in  the  argument  list. 

Input  Output  Variables 


MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

A(lOOOO) 

Complex  work  array 

RA(2,100i 

0) 

Equivalenced  real  work  array  in  which  all 
the  computations  are  done 

NO 

n 

IN 

ARG 

Total  number  of  unknowns,  i.e.,  size  of  the 
square  DT -matrix 

MD 

Total  number  of  righthand  sides  in  the 
system  of  simultaneous  linear  equations 

solved  in  SOLVIT 

KD 

Work  array  size  (total  real -variable 
dimension)  -  present  value  is  8000 

NI 

Tape  number  assigned  to  logical  tape  unit 
containing  all  rows  of  the  augmented 
matrix  [DT | WT] 

*where  n  =  number  of  unknowns  and  m  =  n  +  the  number  of  right  hand  sides 
for  which  the  solutions  are  obtained. 
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For  the  present  assignment  of  tape  units,  see  Section  4.1. 


Calling  Subroutine  MAIN 


5.8  Sepent  8 


5.8.1  SUBROUTINE  BFSMAT  (ND,  NE,  NB,  NP,  NTP,  N TOTAL,  IQ,  IPRNT,  NAS, 

FMACH,  YB,  ZB,  YS,  ZS,  X,  DELX,  EE,  XIC,  SG,  CG,  AR,  RIA,  NBEA1 ,  NBEA2, 
NASB,  NSARAY,  NCARAY,  BFS,  AVR,  CBAR,  AO,  XIS1 ,  XIS2,  KR,  NSBEA, 

IBFS) 


Functional  Description 

This  subroutine  is  the  basic  calling  routine  that  forms  the  [FZ]^ 
and  [FY]^  matrices.  (See  equation  (2.6-35  of  Vol .  1.  These  matrices 

are  written  on  tape  a  row  at  a  time,  alternating  first  a  row  of  FZ  then 

one  of  FY,  etc.  This  set  of  matrices  is  formed  for  each  body.  The  element 

FZ. .  gives  the  force  in  the  z-direction  due  to  either  (1)  a  lifting  surface 

*  J 

box  or  (2)  an  axial  doublet.  The  element  FY^.  gives  the  force  in  the 
y-direction.  The  formulas  differ  depending  on  whether  a  lifting  surface 
box  or  doublet  is  considered.  The  point  pressure  doublet  is  ACp6A  for  the 
case  of  the  lifting  surface  boxes.  The  point  pressure  doublet  for  an  axial 
doublet  involves  a  derivative  with  respect  to  x. 


Input  Output  Variables 


MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

ND 

<5 

Symmetry  flag 

NE 

e 

Ground  effect  flag 

NB 

Number  of  bodies 

NP 

NB 

Number  of  panels 

NTOTAL 

2 1  NSBEA. 

2x  (total  number  of  slender  body 

i=l 

IN 

ARG 

elements) 

10 

Logical  tape  unit  on  which  rows  of 
the  BFS  matrix  are  written 

IPRNT 

Print  flag 

NAS 

Array  containing  the  number  of 
associated  bodies  for  each  panel 

i 
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MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

FMACH 

Mach  Number 

YB 

Array  of  y-coordi nates 

of  the  bodies 

ZB 

ZB 

Array  of  2-coordinates 

YS 

y 

Array  of  y-coordi nates 

of  strips  and 
bodies 

ZS 

2 

Array  of  2-coordinates 

X 

X 

Array  of  the  3/4-chord  locations  of 
boxes  and  midpoints  of  interference 
body  elements 

DELX 

Ax 

Array  of  lengths  of  boxes  and 
interference  body  elements 

EE 

e 

IN 

ARG 

Array  of  the  semi -widths  of  strips 

XIC 

?c 

Array  of  1 /4-chord  locations  of 
all  boxes 

SG 

siriYs 

Array  of  the  sines  of  the  dihedral 
angles  of  strips 

CG 

COSYs 

Array  of  the  cosines  of  the  dihedral 
angles  of  strips 

AR 

Array  of  the  cross  sectional  aspect 
ratios  of  the  bodies 

RIA 

Array  of  the  radii  of  interference 
body  elements 

NBEA1 

Array  of  the  number  of  interference 
body  elements  per  body 

NBEA2 

Z-y  orientation  flag  array  per  body 

153 


MNEMONIC  SYMBOL  IN/OUT  SOURCE 


DESCRIPTION 


Array  of  the  bodies  associated  with 
panels 

Array  of  the  number  of  strips  per 
panel 

Array  of  the  number  of  chordwise 
boxes  per  panel 

One  row  of  the  z-  and  y-forces  in 
the  BFS  matrix 

Array  of  the  average  radii  of  all 
bodies 

Reference  chord 

Array  of  the  radii  of  slender 
body  elements 

Array  of  the  slender  body  element 
leading  edge  coordinates 

Array  of  the  slender  body  element 
trailing  edge  coordinates 

Reduced  frequency 

Array  of  the  number  of  slender  body 
elements  per  body 

Option  flag  =  1  to  select  subroutine 
FZY2  for  the  computation  of  the 
individual  force  element  contri¬ 
butions 
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Calling  Subroutine 

MAIN 

Called  Subroutine 

FWMW 

Equations 


The  forces  in  the  z-  and  y-directions  on  the  slender  body  elements 
are  related  to  the  lifting  surface  box  pressures,  ACp,  and  axial  doublet 
strengths  through  the  [FZ]  and  [FY]  matrices. 

|Fz}'b)  =CFZ](b)  jP* 

\ 

(Fy|(b)  =[FV](b)  Jp[ 
where  {p}  =  (  Acp 

*  ~(y) 

The  elment  FZ^  =  FWMW^  aA^  when  the  sending  element  is  a  lifting 
surface  box. 


FZ 


ij 


=  FWMw!^ 

•  J 


*ik— •  J 


3+i  k~-j  +1 
c 


where  the  superscript  (z)  on  FWMW^ .  indicates  that  the  force  in  the 
z-direction  is  desired,  i.e.,  IF1  =  1.  If  the  y-direction  is  desired 
set  IF1  =  3.  The  subtraction  indicated  above  numerically  performs  the 
differentiation  required  of  and  p^. 


5.8.2  SUBROUTINE  FZY2  (XIJ,  XI,  X2,  ETA,  ZETA,  YB,  ZB,  A,  BETA2,  CBAR, 
K,  FZZR,  FZZI ,  FZYR ,  FZYI,  FZYR,  FZYI,  FYYR,  FYYI,  MFLG) 


Functional  Description 

This  subroutine  calculates  the  forces  in  the  z-  and  y-directions 
on  all  slender  body  elements  of  circular  cross  sections  due  to  a  unit 
pressure  doublet  located  either  inside  or  outside  of  the  cross  section. 
Subroutine  FZY2  is  called  by  subroutine  FWMW,  which  in  turn  is  called  by 
BFSMAT  and  main  only  if  the  option  flag  IBFS  =  1  (card  input);  it  is 
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bypassed  for  all  other  cases  (IBFS  =  0). 


Input  Output  Variables 


MNEMONIC  SYMBOL  IN/OUT  SOURCE 


DESCRIPTION 


1 /4-chord  x-coordinate  of  slender 
body  element 


Leading 


Trailing 


edge  of  slender  body 
element 


y-coordinate 


z-coordinate 


y  coordinate 


z  coordinate 


of  sending  point 


of  body  centerline 


Radius  of  slender  body  element 


BETA2 


1-M  ,  where  M  =  Mach  Number 


Reference  chord 


K,  KR 


Reduced  frequency 


part  of  the  F' 


Imaginary  I 

Real  part  of  the  F* 

Imaginary  Eq*  (5. 8. 2-2) 
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MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

FYYR 

Real 

part  of  the  F  ^ 

OUT 

ARG 

Eq.  (5. 8. 2-3) 

FYYI 

Imaginary 

MFLG 

IN 

ARG 

Detail  print  flag 

TEST! 

The  number  1/7 

TEST  2 

IN 

FZY2 

The  number  1/2 

KBAR 

E 

Eq.  (5. 8. 2-9) 

DX 

AX 

X2  "  X1 

RAA 

RaA 

Eq.  (5. 8. 2-6) 

DELTA 

6 

“/RaA 

XA 

xa 

(x-j  +  X£)/2 

RWIG 

'V 

r 

Eq.  (5.8. 2-7) 

RAIJ 

r3ij 

Eq.  (5. 8.2-5) 

QR 

Eq.  (5. 8. 2-8) 

QI 

CAPA 

A 

Eq.  5.8.2-15) 

CAPDR 

Re(A) 

Eq. (5.8.2-12)  or  Eq.  (5.8.2-16) 

CAPO  I 

Im(A) 

FTHR 

Re(f  J6)) 

0 

Eq.  (5.8.2-10) or  Eq. (5.8.2-13) 

FTHI 

FRR 

Re(f  (r)) 

,  ,,r(r) 

Eq.  (5.8.2-11)  or  Eq.  (5.8.2-14) 

FRI 


MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

CTH 

STH 

cose 

sine 

IN 

FZY2 

*3* 

t 

CM 

CO 

LT> 

cr 

LU 

11 

through 

Ill 

h 

through 

'ii 

The  integrals  I-j  through  1^ 

Calling  Subroutine  FWMW 


Equations 

Subroutine  FZY2  computes  the  unsteady  body  forces  for  circular  cross 
sections  according  to  the  following  equations. 


FlZu  =  cos2ef0^0^  +  sin2efr^ 

FZYij  =  FYZij  =  cose  sine  (fr(r)  "  f0(9^ 
,  =  sin20f^6^  +  cos20  f  ^ 


FYYi j  =  si"  0fe 


where 


nj  -  YB. 

cos©  =  - -  ,  sin© 


(5.8. 2-1) 

(5.8. 2-2) 
(5.8. 2-3) 


aij 


a  -  ZBi 

r  •  • 
aij 


(5. 8. 2-4) 


and 


raij  “  <"j  -  ¥Bi)2  +  (tJ  -  ZBf>2 


(5.8 .2-5) 


where  i  is  the  index  of  the  receiving  body  element 

and  j  is  the  index  of  the  sending  element. 

Note  that  we  set 

cose  =  1.0  and  sine  =  0.0  whenever  r  •.=  =  0.0. 
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The  formulae  for  f0^  and  f  ^  are  calculated  depending  on  the  location 
of  the  sending  element,  as  shown  below. 

Define  a  number  6  as  follows 


if  r 


2 

aij  >  a 


2 


otherwise 


(5.8. 2-6) 

(5. 8. 2-7) 


AX  ■  X2j  -  X„ 

..  _  xli  +  x2i 
XA - 2 - 

Also,  define 

Q  =  m  e  1?a  CM(XA  '  '  W  (5. 8. 2-8) 

where 

k  =  2  k  M  ai/c  (5. 8. 2-9) 

Then,  if  6  _<  1/7 

f0(0)  =  Qa(62al]  +  ifc  I4)  (5.8.2-10) 

and 

fr^r)  =  fe(9)  +  a  (5.8.2-11) 


where 


QraijC-3e2a(e2aI6  +  ik  I9)  +  k2I-| 3  if  rj.  >  a2 


otherwise 


(5.8.2-12) 


I4  ■  </Ra 

>6  ■  ‘*1 
h  '  5'Ra 


If  6  >  1/7 


f0(e)  =  Qa[(B2aI1  -  V-AIg)  +  ik(AI2  +  I4  -  - 
0  1  62a  5  2  4  2R 

fr(r)  =  fQ(e)  +  L 


(5.8.2-13) 


(5.8.2-14) 


where 


A  -  M-(xfl  -  £j)/Rafl 


(5.8.2-15) 


Qra1j  +  +  3wlo> 


+  i302ak(-AI?  + 


i8b2?2 


-  y +  <- 


k"3  AL 


otherwise 


2  2 
if  r  . .  >  a 
aij 


(5.8.2-16) 


and  where  the  integrals  I-j  through  1^  are  calculated  two  different  ways, 
depending  on  the  value  of  the  number  5: 


If  6  <1/2 
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I5  =  -t63/6 

Vft-  [' +  zr  ♦  7'2) 

aA  L 

17  =  fy-  (  >5t/12) 

aA 

18  =  «3/02  Rj> 
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In  =  In 


x2  -  5  +  Ra2 


xn  -  5  +  R 


al 


*3  =  hi  '2(xa  "  ?)I2  ■  RaA  h 


I  =  ^ 
a4  * 
Br 


Atan  x2  "  5  -  Atan  X1  “  5 


er 


-•X/ 

er 


(&"■£■)'  <XA'5)'6 


T  -  ■>/  1 

I? 

a2  •'al 


!8  =  h  "2(xA  "  "  RaA  *6 

(x2  -  5)  (x1  -  e) 


l9‘S¥ 


a2 


2 — +  U 


al 


l10 


2e2r2 


<XA  ~  {)ax/2  *  RaA  ,  <»A  ~  {)tx/2-RaA 


+  (> 


a2 


al 


where 


Ral  =  ^(x]  -  z)Z  +  e2r2  and  Rfl2  =  ^ (x2  -  c)2  +  e2r2 


'A  -  «>I 
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5.9  Segment  9 


5.9.1  SUBROUTINE  BFM  (I WORK,  RWORK,  WORK,  NWORK,  NPTAP,  NPSTAP, 

NBFM,  NPOT,  NM,  IPRNT,  IERR0R,  IBFS) 

Functional  Description 

Each  and  every  singularity  in  the  flow  field,  whether  it  be  inside  or 
outside  of  a  body,  contributes  to  the  force  distribution  on  a  body. 

Subroutine  BFM  determines  this  loading,  force  and  moment,  on  bodies  that 
occurs  due  to  these  singularities.  The  following  are  the  flow  singularities 
considered. 

(a)  Slender  body  singularities. 

(b)  Interference  element  singularities. 

(c)  Lifting  surface  boxes  external  to  the  body. 

(d)  Lifting  surface  box  images  both  inside  and  outside  of  the  body. 

(e)  The  additional  contributions  due  to  considerations  of  symmetry 
and  ground  effects  for  the  above  singularities. 

Subroutine  BFM  is  essentially  a  colling  routine.  Generally  speaking,  FWMW 
determines  the  loading  on  the  body  for  a  pressure  singularity  of  unit 
strength,  including  images,  symmetry,  and  ground  effect.  Subroutine  ORGAN 
gives  the  region,  on  the  b^dy,  over  which  the  loading  acts.  Subroutine 
SBLOAD  determines  the  pressure  singularity  strength  from  the  various  flow 
singularity  strengths  and  gives  the  detailed  loads  on  each  slender  body 
element.  BFM  then  uses  the  unit  loading  obtained  from  FWMW  and  the  results 
from  SBLOAD  in  conjunction  with  ORGAN  to  find  the  load  (both  forces  and 
moments)  on  each  of  the  slender  body  elements.  The  slender  body  elements 
have  been  chosen  as  a  convenient  set  of  elements  over  which  all  body  forces 
are  distributed. 
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MNEMONIC 


SYMBOL 


IN/OUT  S' 


DESCRIPTION 


AR 

At 

■1 

Aspect  ratio  of  the  bodies. 

CG 

cos(y) 

Cosine  dihedral  angle. 

EE 

e 

IN 

NB 

ND 

NB 

5 

Number  of  bodies. 

Symmetry  flag 

NE 

c 

Ground  effects  flog 

— 

HK88M 

IN 

ARG 

Number  of  modes. 

NP 

NP 

Number  of  panels. 

SG 

sin(y) 

Sine  dihedral  angle. 

YB 

yc 

y  coordinate  of  body  center. 

YP 

y 

y  coordinate  of  centerline  of  strip 

ZB 

2c 

IN 

z  coordinate  of  body  center. 

ZP 

z 

z  coordinate  of  centerline  of  strip 

AVR 

Average  radius  of  bodies 

NAS 

Number  of  associated  bodies  per  panel 

1  NBY 

Number  of  y-oriented  bodies. 

j  NO7 

Number  of  z-oriented  bodies. 

IN 

Blank 

Common 

Total  number  of  boxes  and  interference 

body  elements. 

NTP 

Block 

Total  number  of  boxes 

NTY 

Number  of  y-oriented  interference 
body  elements 

NTZ 

Number  of  z-o. iented  reference 

body  elements 

XLE 

\E 

IN 

Leading  edge  of  bodies 

XTE 

*LE 

Trailing  edge  of  bodies 

DELX 

Length  of  boxes  and  interference 
body  elements 

MACH 

M 

Mach  number 

NASB 

Number  of  the  associated  bodies  per  panel 

NBFA 

Number  of  interference  body  elements 
per  body  and  body  orientation. 

IPRNT 

Print  flag 

IERR0R 

IN 

ARG 

Error  flag 

IBFS 

Body  force  calculation  method  flag 

MNEMONICS 

SYMBOL 

IQS 

SOURCE 

DESCRIPTION 

NBFM 

IN 

ARG 

I/O  unit  which  output  body  and  forces 
are  to  be  written  on  (unformatted) 

Print  output  unit 

NSBE 

Number  of  slender  body  elements  per  body 

NTYS 

NTYS 

Number  of  y-oriented  slender  body  elements 

NTZS 

NTZS 

Number  of  z-oriented  slender  body  elements 

WORK 

Complex  array  for  temporary  storage  of 
results.  The  length  in  complex  words 
required  is  4  (No.  Strips  +  NB2 
+  NBY)  +  4 (NTZS  +  NTYS) 

XIS1 

Sc 

si 

?s2 

Leading  edge  of  slender  body  elements 

XIS2 

IN 

Training  edge  of  slender  body  elements 

IWORK 

ARG 

Integer  array  for  temporary  storage  of 
results  its  length  is 

2(NTP  +  NTZ  +  NTY  +  NTZS  +  NTYS)  words 

NPTAP 

1/0  unit  containing  the  P  matrix 

NWORK 

Size  of  the  array  WORK 

RWORK 

Real  array  for-  temporary  storage.  The 
length  is  =  NTP  +1+2  {maximum 

C (NTP  +  NTZ  +  NTY)  or  (NTZS  +  NTYS)]} 

NCARRY 

Array  containing  the  number  of  chordwise 
boxes  per  strip  in  a  panel 

NPSTAP 

IN 

ARG 

I/O  unit  number  containing  the  PS  matrix 

NSARRY 

Array  containing  the  number  of  strips 
in  a  panel. 

Calling  Subroutines  MAIN 

Called  Subroutines  and  Comnon  Blocks 


RWREC,  WRTFMU,  FWMW,  ORGAN,  SBLOAD,  WRTFMF,  Blank  Common 
Error  Messages 

SUBROUTINE  **BFM**  NWORD(NNNNNNNN)  IS  GREATER  THAN  NWORK(NNNNNNNN) 

***J0B  TERMINATED*** 

The  available  working  array  is  too  small  to  execute  the  subroutine.  Either 
decrease  the  problem  requirements  or  increase  the  size  of  the  work  area. 
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5.9.2  SUBROUTINE  FMZY  (DYB,  PZB,  DA,  PAR,  DY,  DZ,  DKR,  DM,  DCBAR, 

DRFZZ,  DIFZZ,  DRFZY ,  DIFZY,  DRFYZ,  DIFYZ,  DRFYY,  DIFYY,  DRMZZ, 
DIMZZ,  DRMZY ,  DIMZY,  DRMYZ,  DIMYZ,  DRMYY,  DIMYY,  I  FI ,  IF2) 

Functional  Description 

This  subroutine  calculates  the  force  and  moment  on  an  elliptic  cross 
section  due  to  a  unit  pressure  doublet  located  either  inside  or  outside  of 
the  cross  section. 


Input  Output  Variables 


MNEMONICS 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

DYB 

DZB 

DA 

DAR 

DY 

DZ 

DKR 

DM 

DCBAR 

YB 

ZB 

a 

AR=b/a 

y  and  y 
z  and  z 

k 

M 

c 

IN 

ARG 

Origin  of  ellipse 

Width  of  ellipse  in  y-di recti  on 

Ratio  of  semi  height  to  semi  width 
of  the  ellipse 

Coordinates  of  pressure  doublet.  For 
bared  quantity  y=y-yB,  z=z-zB 

Reduced  frequency,  wc/ZU^ 

Mach  number 

Reference  chord  length 

DRFZZ 

LU’  ‘ 

Real  and  imaginary  parts  of  z-force  due 

DIFZZ 

to  doublet  oriented  in  z-di recti on 

DRFZY 

F  (y) 

Real  and  imaginary  parts  of  z-force 

DIFZY 

due  to  y-doublet 

DRFYZ 

C  (z) 

‘  w 

Real  and  imaginary  parts  of  y-force  due 

DIFYZ 

y 

to  z-doublet 

DRFYY 

F  (y) 

OUT 

ARG 

Real  and  imaginary  parts  of  y-force 

DIFYY 

y 

due  to  y-doublet 

DRMZZ 

My(2) 

Real  and  imaginary  parts  of  z-moment 

DIMZZ 

due  to  z-douulet 

DRMZY 

„  (y) 

etc 

DIMZY 

z 

DRMYZ 

Mw(z) 

etc 

DIMYZ 

a 

DRMYY 

M  <y) 

etc 

DIMYY 

y 
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MNEMONIC 

IF! 


SYMBOL 


IN/OUT  SOURCE 
IN  ARG 


IF2 


DTHETO 


0 


IN 


ARG 


o 


_ DESCRIPTION _ 

1  if  body  z-oriented,  2  if  body  z  and  y 
oriented  and  3  if  y-oriented 
0  if  doublet  outside  of  ellipse 
1  if  doublet  inside  of  ell i pse 
Elliptic  coordinate  at  which 


FMUR 

FMUI 

PMUR 

PMUI 

DKJO 

DKYO 

DKJ1 

DKY1 


f 


kJ_(k) 

kY?(lc) 

kJ?(k) 

kYjdc) 


EF1 

EF2 

DKK 

DRR 


cosA 

sinA 

k 

R 


RCURV  a 

DEE  e 


IN 


FMZY 


integration  starts 

Real  and  imaginary  parts  of  f 

Real  and  imaginary  parts  of  p 


J  ,  Y  are  Bessel  function  of  order  v 
Y  y  ' 


A  is  the  surface  slope  of  the  elliptic 
cross  section  at  the  point 

Distance  from  field  point  to  surface 
point 

Surface  curvature  of  the  ellipse 
Element  half  width 


Calling  Subroutines  FWMW 

Equations 


If  km  =  0  and  /R  =  1  and  r/a  >  I 
Fz(z)  =  (y2  -  *2)  y  -P5- 


Fz^  =  -  (z  y)  a2/r'1 

F  (z)  =  F  (y) 

y  z 

F  (y)  -  F  (z) 

y  z 
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0 


M  <2>  .  M  <»)  =  M  (2)  =  H  <*>  = 

z  z  y  y 

z  =  z  -  zB 

y  =  y  -  yB 

If  kM  =  0  and  r/a  <  1 

Fz(z)  =  1/  (1  +  PR) 

Fy(y)  =  PR/  (1  +  PR) 

F  (y)  =  p  (z)  =  m  (2)  =  M  {y)  =  M  (z)  =  M  {y) 

z  y  z  z  y  y 

Otherwise  numerical  integration  is  performed 

©o  =  0j  if  z  and  y  >  0 

0q  =  -  0!  if  y  >0  and  z  <  0 

0  =  0i+Trify<O  and  z  <  0 

0 

0  =  tt-01  if  y  <0  and  z  >  0 

0 

cos2©!  =  1  -  sin2©! 

sin20!  = 

A  =  (y  -  yB)2  +  (z  -  zB)2  +  a2(PR2  -  1) 

If  |  1  -  PR2|  <  .001 

sin20i  =  (z  -  zB)2  /  r2 
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0 


0,  +  A  0 


A0  =  2tt  /N 
c  U)  _ 


(Currently  N  =  (NUMB) 


j=l 

N 


f  (k.)  EFZZ. 


Fz(y)  =  J  f  (kj)  EFZY. 

jN] 

FyU)  =  f  (Eo>  EFyzd 

K! 

c  (y) 


>  f  (kJ  efyy. 

fa  3  3 

N 

M_V4'  =  y  P  (k.)  EF1.EIM2, 

0=1  3  3  3 
N 

2  p  <u 

o=i 
N 

J  -  P  (kj)  EF2j  EIM2j 

j=l 

N 

m >>  .  y  - 


M 


M, 


(2) 


(y)  ... 


(z) 


EF1.EIM1  • 
J  J  J 


2  - p  <V  EF2j  EIM1j 

0=1 


f ■  -  Hi  v‘rkd> +i  j‘(ko>l 


p<ko>  ■  TEJ  °o  <k0>  -  1  V0  (kj) 


=  30) 


169 


Y  and  are  Bessel  functions  of  order  v  and  are  approximated  by  polynomials 
taken  from  Reference  4 


k 

2  k  M  R/c 

EF1 

COSA 

EF2 

sinA 

EIM1 

2e  (y  -  n)/R 

El  M2  = 

2e  (z  -  c)/R 

R2 

(y  -  n)2  + 

e  =  "  t  /c2  +  ar4  n2 


n  =  a  coso 
e  =  b  si  no 


COSA 

y?2  +  PR4  n2 

sinA 

-  PR2  n 

y?2  +  PR4  n2 

If 

R/e  >  5.0 

EF 

»> 

k_2e 

K  -RZ- 

IF 

1.5  >  R/e  <  5.0 

EF 

_  £  2e  ,2  ,eu  fT  R 
"  K"F”  +  3  (R}  [l  a  " 

BJ  ,  K 
R"  +  R 

<1 


) 
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!  (|)5[I  {-<n  +  ^  ■  3(fw  -  D/a)  ♦  <fr} 

+  j{2|  (1  -  D/a)  -  (B/R)3 1  +  l|{(B/R)2  -  1  +  D/aJj 


IF  R/e 

<  1.5 

EF  = 

{- 

OB  ,  77  3  DJB  ,  DK 
2  K  “  4  a  2a 


I  }  l  *ta„  (lh3) 


+ 

2el  . 

T“  + 

1 

+  R2 

|  eBj  (1  +  D/a) 

+ 

t 

! 1  = 

3/?  sin2x  -  1/2 

EFZZ 

=  EF 

with 

j  = 

sinx  cosx  -  inx  (■ 

K  = 

-  cosx  (z  -  c) 

r 

;  i  ■ 

3/2  COSySiny 

EFYZ 

=  EF 

with  ' 

o  = 

cos2x  -  sinx 

IL. 

K  = 

-cosx  (y  -  n) 

' 

I  = 

3/2  cosxsinx 

EFYZ 

=  EF 

with 

!  o . 

-sin2x  -  cosx 

IL. 

K  = 

sinx  (z  -  c) 

with  ! 

i 

i  I  = 

3/2  cos2x  -  1/2 

EFYY 

=  EF 

0  = 

-  cosxsinx  - 

cosx 

R  = 

sinx  (y  -  n) 

t* 


ilz  a), 

a 


B 

D 

a 


-2  (cosx  (y  -  n)  +  sinx  (z  -  c) 
sinx  (y  -  n)  -  cosx  (z  -  s) 

2/3 


|  sin2©  +  /R2cos2© 
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5.9. 


SUBROUTINE  ORGAN  (I START,  ISTOP>  NLBE,  NTP,  X,  DELX,  PERCNT , 
XLE,  XTE,  XIS1,  XIS2,  ITYPE) 


Functional  Description 

It  is  assumed  that  the  effect  of  a  flow  singularity  in  a  body  acts  exactly 
at  the  same  x-location  as  the  flow  singularity  itself.  This  routine  there¬ 
fore  generates  two  output  arrays,  one  identifying  the  slender  body  elements 
which  the  leading  edges  of  the  sending  elements  lie  within  and  the  other 
the  trailing  edges. 


Input  Output  Variables 


MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

X 

X 

IN 

AR6 

ITYPE^l 

Array  of  sending  box  or  body  element 

location 

ITYPE  =  1 

Array  of  sending  slender  body  elements 
leading  edge  x-location. 

NTP 

NTP 

IN 

AR6 

Number  of  boxes  or  elements  sending 

loads 

XLE 

\.E. 

Leading  edge  of  receiving  slender  body 

XTE 

*r.E. 

Trailing  edge  of  receiving  slender  body 

DELX 

IN 

ARG 

ITYPE^l 

Array  of  sending  box  or  body  element 
length 

ITYPE  =  1 

Array  of  sending  slender  body  elements 
trailing  edge  x-location 

NLBE 

IN 

ARG 

Number  of  elements  in  the  receiving 
slender  body. 

XIS1 

Si 

Array  of  leading  edge  x-coordi nates 

of  the  receiving  slender  body  elements 

XIS2 

?s2 

Array  of  the  x-coordi nates  of  the 

trailing  edges  of  the  elements  of  the 
receiving  slender  body 

2 


ITYPE 


ISTART 


I  STOP 


PERCNT 


IN 


ARG 


OUT 


ARG 


OUT 


ARG 


IN 


ARG 


Flag  indicating  the  type  of  the 
sending  elements 

tl 

s ending  elements  are  either  boxes  or 
interference  bodies. 

=  1 

sending  elements  are  slender  bodies. 
An  array  containing  the  first  slender 
body  element  receiving  a  load  con¬ 
tribution  from  the  sending  element. 
The  length  of  this  array  is  (NTP) 

An  array  containing  the  last  slender 
body  element  receiving  a  load  con¬ 
tribution  from  the  sending  element. 
The  length  of  this  array  is  (NTP). 

The  fractional  location  of  the  input 
x  with  respect  to  the  length  of  the 
sending  element.  This  is  used  to 
calculate  the  leading  edge  of  the 
sending  element. 


Calling  Subroutine  BFM 

5.9.5  SUBROUTINE  SBLOAD  (COEF,  IFIRST,  HAST,  KSt  XIS1,  XIS2,  DELTA,  K, 
PERCNT,  XLE,  XTE,  NBE,  FWZ,  FWV,  MMZ,  MWY,  FZ,  FY,  MZ,  MY 

Functional  Description 

It  is  assumed  that  the  x-distribution  of  the  flow  singularity  is  constant 
over  the  flow  of  the  flow  singularity  element.  It  is  further  assumed  that 
the  effect  of  a  flow  singularity  on  a  bcdy  acts  exactly  at  the  same 
x-location  as  the  flow  singularity  Itself.  This  routine  distributes  the 
forces  and  moments  on  a  body  due  to  internal  and  external  flow  singularities. 
The  center  of  force  in  each  of  the  slender  body  elements  is  also  calculated. 
The  output  from  subroutine  ORGAN  determines  which  slender  body  elements  are 
affected. 
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Input  Output  Variables 


MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

— 

IN 

ARG 

x-coordinate  of  sending  element  or  box 

FY 

wm 

OUT 

y-force  array 

FZ 

WM 

OUT 

2-force  array 

KS 

mm 

IN 

ARG 

reserved 

MY 

My 

OUT 

y-moment  (yawing)  array 

MZ 

M- 

2 

OUT 

m 

z-moment  (pitching)  array 

FWY 

Fwy 

unit  y-force  on  sending  body  or  strip 

FWZ 

Fwz 

IN 

ARG 

unit  z-force  on  sending  body  or  strip 

MWY 

Mwy 

unit  y-moment  on  sending  body  or  strip 

MWZ 

Mwz 

unit  z-moment  on  sending  body  or  strip 

NBE 

N 

Number  of  body  elements  or  receiving 

body. 

XLE 

IN 

ARG 

Reserved 

XTE 

Reserved 

COEF 

XIS1 

«si 

Array  of  leading  edge  x-coordinate 

XIS2 

5s2 

IN 

ARG 

Array  of  trailing  edge  x-coordi nates 

of  the  slender  body  elements. 

DELTA 

AX 

Length  of  sending  box  or  element. 

I  LAST 

Last  slender  body  element  to  which  the 

loads  are  to  be  applied. 

PERCNT 

PERCNT  =  (x  -  xle)/DELTA 

where  is  the  x-coordinate  of 

the  leading  edge  of  the  sending 

element  or  box. 

Calling  Subroutines  BFM 
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5.10  Segment  10 

5.10.1  SUBROUTINE  AERO  (NMODE.  NSTRIP,  NW,  NBFM  ,  NEWBFM,  IBFS) 

Functional  Description 

Subroutine  AERO  computes  the  aerodynamic  parameters  for  each  mode.  These 
include  the  sectional  lift  and  moment  coefficients  for  all  strips  of  all 
lifting  surface  panels,  and  the  center  of  pressure  locations  (Eqs.  5.10.1-1 
through  5.10.1-3):  lift  and  moment  coefficients  (z-  and  y-components)  for 
all  slender  body  elements  (Eqs.  5.10.1-4  through  5.10.1-9);  total  lift  and 
moment  coefficients  for  each  slender  body  (Eqs.  5.10.1-10  through  5.10.1-13); 
and  the  total  lift  and  moment  coefficients  including  body  effect  (Eqs. 

5.10.1-14  through  5.10.1-18).  The  above  aerodynamic  parameters  are  printed 
along  with  the  strip  number,  or  body  element  number  in  case  of  a  slender  body, 
and  the  y-  and  z-coordi nates  of  the  strips  and  slender  bodies. 

Input  Output  Variables 


MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

NMODE 

Total  number  of  modes 

NSTRIP 

Total  number  of  strips  on  all  lifting 

surface  panels 

NW 

IN 

Argument 

Tape  number  for  logical  tape  unit  con- 

List  of 

taining  the  solutions  (aC- )  for  all 

AERO 

unknowns*  ' 

n1  and  all  modes 

NBFM 

Tape  number  for  logical  tape  unit  con- 

taining  all 

slender  body  element  forces 

and  moments 

for  all  modes  when  IBFS  =  0 

NEWBFM 

Tape  number  for  logical  tape  unit 

containing  all  slender  body  element 

forces  when 

IBFS  =  1 

IBFS 

Body  force  calculations  method  flag 

FZ(200) 

3Vq 

z-forces  -> 

FY(200) 

v 

Tape 

y-forces 

for  all  slender  body 

NBFM  or 

•  elements  and  all 

MZ( 200) 

NEWBFM 

z-moments 

modes 

MY (200) 

IN 

y-moments  - 

*n  =  total  number  of  boxes  +  total  number  of  interference  body  elements. 
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MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

DCP(5C0) 

ACp 

Tape  MW 

Pressure  coefficients  (solutions)  for 
all  unknowns  and  all  modes 

CN(200) 

C"qj 

See  Equation  (5.10.1-1) 

CM  (200) 

See  Equation  (5.10.1-2) 

CPR(200) 

CPRqj 

See  Equation  (5.10.1 -3a ) 

CPI (200) 

CPIqj 

See  Equation  (5.10.1 -3b) 

FZLB 

c(b), 

Lz 

See  Equation  (5.10.1-4) 

FYLB 

r(b), 

cy 

See  Equation  (5.10.1-5) 

MZLB 

c(b)- 

mz 

See  Equation  (5.10.1-6) 

MYLB 

c(b)i 

my 

OUT 

AERO 

See  Equation  (5.10.1-7) 

DFZX 

3fzt/3x 

mm 

aee  tquanon 

ins 

DFYX 

See  Equation  (5.10.1-9) 

Q 

CZB(IO) 

c(b) 

CZq 

See  Equation  (5.10.1-10) 

CYB(IO) 

c(b) 

See  Equation  (5.10.1-11) 

CMG(IO) 

f(b) 

lM 

Se:  Equation  (5.10.1-12) 

CNB(IO) 

c(b) 

See  Equation  (5.10.1-13) 

2e  cos 

Y 

■ 

2e  sin 

Y 

GUCJ 

GJ 

IN 

1  or  0.5;  Equation  (5.10.1-19) 

GLCB 

9o» 

1  or  0.5;  Equation  (5.10.1-20) 

SYMA 

AERO 

1  +  6 

where  6  is  the  symmetry 

SYMB 

1  -  5 

flag  (input) 

CZT 

°Zq 

Bee  Equation  (5.10.1-14) 

CYT 

CYq 

OUT 

See  Equation  (5.10.1-15) 

CMT 

CM 

See  Equation  (5.10.1-16) 

CNT 

CN 

See  Equation  (5.10.1-17) 

CLT 

ii _ 

See  Equation  (5,10.1-18) 

CLT 


Calling  Subroutine  MAIN 

Common  Blocks  The  Blank  Common  Block 


Equations 

The  following  aerodynamic  parameters  are  calculated  for  each  mode: 


A.  Lifting  Surface  Strips 


AX. 


(5.10.1-1) 


*2 

0  •  *  • 

J  1=1] 

Re<c,»  >  , 

CPRqj=WdVT  (5-10-1-3a) 

qj; 

lB<v>  i 

CPIqj ■  isrtfV + 1  (5-10-1-3b) 

where 

i*1  and  ig  are  the  indices  of  the  first  and  last  boxes  in  strip  j, 
where  j  =  1,  NSTRIP; 

Cj  is  the  chord! ength  of  strip  j,  and 
-  514.  is  the  x-coordfnate  of  the  1 /4-chord  point  on  the  centerline  of 

J 

strip  j. 

B.  Slender  Body  Elements 

f(b) 

fc(b)»  =  ifl*  (5.10.1-4) 

Z  AxSB[b) 


177 


(5.10.1-5) 


and 


where 


c(b). 

y 


*(b) 

yqt 

~Tb r 


AX  SB 


m 


(b) 

Vt 


AXSB 


m 


m. 


(b) 

qt 

IbT 


AxSB 


(5.10.1-6) 


(5.10.1-7) 


AxSbI^  =  §S2ib^  -  ?Slib^  =  length  of  slender  body  element  1 1 1  in 
u  '  body  1 b 1 ;  see  Blank  Common,  Sec.  3.1. 

In  addition  to  the  above  aerodynamic  parameters,  subroutine  AERO  also  computes 
modified  lift  coefficients  (z-  and  y-components)  for  all  slender  body 
elements  according  to  the  following  equations. 


and 


afzt/3x 


afyt/ax 


=  AR 


gz 


tt 


8fzt 


GY 


tt 


(5.10.1-8) 


(5.10.1-9) 


where  the  elements  of  the  NSBE^  x  NSBE^  matrices  GZ^  and  GY^  are 
defined  as  follows: 

GZ^  =  G(a^b\  k,  M,  c ,  AR,  R  =  RZ^,  x  -  ?  =  xSBt  -  xSB^) 

GYt^  =  G(a^,  k,  M,  c,  PR  R  =  RYtp  x  -  %  =  xSBt  -  xSB^) 


where 


G(a,k,M,c,AR,R,x-c)  =  ^ 


+  .  2kM\J<2kMW) 
Tt 


f  =  C(x  -  ?)M  -  R]/c 
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and 


RZt^  =  ^(xt  -5|-)2  +  ^a2m(l  +  /R  )/2 


RYt^=  ^(xt  -?r)2  +  Ml  +  AO/2 

The  indices  1 t '  and  't'  denote  the  'row'  and  'column'  subscripts  for  the 
matrices  [GZ^]  and  [GY^t-]  corresponding  to  the  slender  body  elements 
' t 1  and  't' ,  and  NSBE(b)  =  number  of  slender  body  elements  in  body  ' b ' . 
The  total  lift  and  moment  coefficients  (z-  and  y-directions)  for  all  slender 
bodies  are  defined  by  the  following  equations. 


■  (b) 


NSBE 


(b) 


cr°jr£  f^b’ 

q  t=1  qt 


'(b)  _  1 
'Y_  "  S' 


NSBE 


(b) 


t=l 


,(b) 

yqt 


NSBE 


(b) 


t«l 
NSBE^y 


eft)  =  1-  V 

N  aF  s-j 

AC  t=1 


-fy  (5(b)  _x(b))  +  m(b> 

yqt^t  XLE  >  +  my_4. 


qt  J 


(5.10.1-10) 


(5.10.1-11) 


(5.10.1-12) 


(5.10.1-13) 


where 

b  is  the  index  of  the  slender  body, 

NSBE^)  is  the  number  of  elements  in  slender  body  ' b ' 

5^)  is  the  x-coordinate  of  element  ' t '  in  slender  body  ' b ' ,  and 

x^}  is  the  leading  edge  x-coordinate  of  slender  body  ' b ' . 
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The  total  lift  and  moment  coefficients,  including  body  effect  are  defined  as 
follows: 


NSTRIP 


qi 


MB 

cz  .<1  ♦,){!■  E  Veicicn  cosYi  ♦  V  g<b>C<»>  1  (5.10.1-14) 

1  p\  J  J  3  U  ’  U  1  ' 


NSTRIP  NB 

cy  =  (1  -  5) \l  E  GJ2eicJcn  sin  Yj  +E  9<b)cYb)  1  I5-'0-1-'5) 

1  jTf  J  J  J  qj  b^i  q  1 


NSTRIP 


j=l 


"1 


C  \i  _J2eJ  cos 


NB  , 

+  E  g(b)[4b)  ~4b)(x[b)  -  XM)/c]  |  (5.10.1-16) 

b=l 


NSTRIP 


CN  =  (1  -6)<^~  V  “CCn  (5l4i  “XM)32ei  sin  Y4 

N  |Ac  TTj'  J  mqj  "qj  J  J  J 

NB 

+  E  9(b)CCj$b)  ~4b)(xLE}  "XM)/5]  }  (5.10.1-17) 

b=l  q 


and 


NSTRIP 


C*=(1 


6)  7s  Z  VCn  >j  cos  +  zj  sin  V2eJ 
j=l  qJ 


b=l 


NB  NB 

V  g<b>C<bUb>  +E  *(b)4b>4b> 

mrnmm  ^  mrntamm  1  n 

b=l 


(5.10.1-18) 


where 


A  =  reference  area 
XM  =  moment  axis 
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<$  =  symmetry  flag 


=  y-  and  z-coordinates  of  the  centerline  of  strip  ' j 1 


=  y-  and  z-coordinates  of  the  centerline  of  slender  body  V 


1/2  if  y.  =  0  and  cosy^  =  0  and  6  /  0 
6.  =  J  J 
J  1  otherwise 


(5.10.1-19) 


,(b)  J 


1/2  if  yc(b)  =  0 
1  otherwise 


(5.10.1-20) 


5.10  Segment  11 

5.11.1  SUBROUTINE  GENF  (NMODE,  NSTRIP,  NW,  NEWBFM,  IMODE,  AA,  NAI,  NPR2.IBFS) 
Functional  Description 

Subroutine  GENF  corputes  generalized  forces  for  all  pressure  and  deflection 
modes  according  to  either  of  the  two  definitions  given  in  Equations  (5.11.1-1) 
(AGARD  definition)  and (5.11 . 1 -2 )  (conventional  generalized  forces)  depending 
on  the  setting  of  the  input  flag  NPR2  -  1  for  AGARD  forces,  2  otherwise.  It 
also  prints  the  pressure  coefficients  ' ACp'  for  all  boxes  of  all  lifting 
surface  panels  along  with  the  panel,  strip  and  box  number,  and  the  fractional 
chordwise  locations  (x/c)  for  the  'ACp' .  Note  that  Subroutine  GENF  can  be 
bypassed  through  MAIN  by  specifying  NPR2  =  0  (input)  whenever  generalized 
forces  and  pressures  are  not  desired. 


Input  Output  Variables 


MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

NMODE 

NSTRIP 

NW 

NEWBFM 

IN 

ARG 

See  AERO,  Sec.  5,10.1 

IMODE 

(2,150,3) 

AA(2, 

150,3) 

NAI (3), 
NA(3) 

NPR2 

Modal  input  array,  IA(2 ,150,3)  in  Sub¬ 
routine  RDM0DE;  see  Sec.  5.4.1 

Floating  point  equivalent  of  the  array 

IMODE 

See  RDM0DE,  Sec.  5.4.1 

Control  flag;  1  for  AGARD  generalized 
forces,  2  for  conventional  generalized 

forces 

FZ( 200) 
FY(200) 
MZ(  200) 
MY( 200) 

■ 

IN 

K 

See  Sec.  3.2 

DCP(500) 

ACp 

IN 

Tape  NW 

See  AERO,  Sec.  5.10.1 
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MNEMONIC 

SYMBOL 

IN/OUT 

SOURCE 

DESCRIPTION 

N8 

N8 

N 

n 

One  set  of  modal  input  for  panel  — 

M 

m 

see  Sec.  5.4.1 

AMODE 

ai 

nm 

X40C 

x4./c 

K 

See  Equation  (5.11.1-4) 

y4^c 

DELA 

AAk 

Area  of  box  ‘k1 

IN 

GENF 

See  Equations (5.11 .1-11  through  -13) 

See  Equations(5.11 .1-14  and  -15) 

az. 

One  modal  coefficient  for  body;  see 

< 

Sec.  5.4.1 

QW ( 50) 

One  row  of  each  of  the  generalized  force 

QZ(50) 

fvj 

O' 

OUT 

GENF 

components  for  all  pressure  modes  and  one 
deflection  mode  -  see  Equations  (5.11.1-3 

QY ( 50) 

-5  and  -8) 

QIJ(50) 

Qij 

See  Equation  (5.11.1-1  and  -2) 

Calling  Subroutine  MAIN 

Common  Blocks  Blank  Common  Block 


Equations 
AGARD  definition 

-  fjj  (QWiJ  +  QZiJ  +  qV 

Conventional  definition 

Qij  ’  (,“ij  +  QZij  +  QYi.i>  (5.11.1-2) 

The  Q^.  are  computed  for  all  deflection  modes  ‘i1  and  all  pressure  modes 
*  j 1 ,  where  i  and  j  run  from  1  through  NMODE. 
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The  three  components  of  the  Q. .  are  computed  by  Subroutine  GENF  as 

*  J 

follows: 


A.  Lifting  Surface  Contribution 


NBOX 

■  E  GlCpkjhkA 


(5.11.1-3) 


where 


=  1  y-coordinate  °f  centerline  of  strip  is  0  and  cos  Ystri-p  =  0 
2  otherwise 

aC  is  the  pressure  coefficient  for  box  1 k '  for  pressure  mode  'j' 
Pkj 

h^  is  the  modal  deflection  of  box  ’ k '  at  the  1/4-chord  x-coordinate 


of  the  panel  box  centerline,  ‘x4[<1,  defined  by 


5  5 


n=0  m=0 


X4,\n/y4,\m 


nm  \  c 


(5.11.1-4) 


where 


I -  2  2" 

y4k  =  y(yz  -  (N8)yi(p))  +  (zz  -  (N8)Zl(b)) 

yz  and  z^  are  the  y-  and  z-coordi nates  of  the  centerline  of  strip  V, 
and  Yl^  and  Zl^  are  y-  and  z-coordinates  of  the  inboard  edge  of 
panel  ' p ' ,  and  aA^  =  area  of  box  1 k ' . 

The  variables  ainm  and  N8  are  modal  input  data;  see  Section  1.2. 

B.  Slender  Body  Contribution 

The  QZ^j  and  QY^-  in  Equations  (5.11.1-1  and  -2)  denote  the  z-  and 
y-components  of  the  slender  contribution  to  the  total  generalized  forces. 


MB  NSBE^/  dh^  \ 


(5.11.1-5) 
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where 


X>i 


,<b>\n 


(5.11.1-6) 


%  s  /*sb<« 

T  =  Z^nazi  \-j- 
n=0  n\  c 


(5.11.1-7) 


(5.11.1-8) 


where 


h<b>  =  c 


E/  XOD,. 

M— ; 

n=0  *  \  c  / 


,W\" 


yt!  _  ...  /XSEt 


TFi=E^1, 


(b)\n  ' 


n  \  c 


(5.11.1-9 


(5.11.1-10) 


and  where  the  variables  aj  ,  az-j  ,  ay-}  and  N8  are  modal  input  data  (see 

j  |v\  *  ""  *  *  *  * 

Section  1.2);  xSB^  '  is  the  x-coord'  nete  of  the  “'ender  body  element 
midpoint  for  element  ' t ‘  of  slender  body  V,  and 

.(b)  J 1  if  f  - 0 

2  if  y£b)  f  0 

Note  that,  to  facilitate  programming,  one  array,  K1A(500)  is  generated  in 
Subroutine  GENF  for  all  boxes  and  all  slender  bod,,  elements,  defined  as 


HQA(K)  =  GhkAAk,  k  =  1 ,  NBOX 

and  hk  is  defined  in  Equation  5.11.1-4  as 
HQA(KZ)  =  g^hkz,  kz=kz],kz2 


(5.11.1-11) 


(5.11.1-12) 
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where 


NBY 

ky-,  =  kz2  +  1,  ky2  =  kz2  +  ^  NSBE^ 

b=l 

and  h^y  is  given  in  Equation  (5.11.1-9). 

Also,  one  array,  DHQ(IOO),  is  generated  to  contain  all  the  dhz/dx  and 
dhy/dx  in  Equations (5.11 .1-7)  and  (5. 11. -1-10)  as  follows: 

dhz 

DHQUz)  =  g^  ,  iz  =  izr  IZ2  (5.11.1-14) 


where 

and 


and 


DHQUy) 


ZZj  =  1 

NBZ 

£z2  =  ^  NSBE^ 
b=l 


*y  =  ^yp  *y2> 


(5.11.1-15) 


where 

and 


jiy-|  =  iz2  +  1 

NBY 

zy2  =  £Z2  +  ]T  NSBE^ 

b=l 
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6.0  PROGRAM  LISTING 
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